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PREFACE 


The Joint University Program for Air Transportation Research is a 
coordinated set of three grants sponsored by NASA Langley Research Center, 
one each with Massachusetts Institute of Technology (NGL-22-009-640) , Ohio 
University (NGR-36-009-017) , and Princeton University (NGL-31-001-252), to 
support the training of students for the air transportation system. These 
grants, initiated in 1971, encourage the development of innovative curricu- 
lums and support the establishment of graduate and undergraduate research 
assistantships and internships. 

An important feature of this program is the quarterly review, one held 
at each of the schools and the fourth at a NASA facility. This latter review 
for 1981 was conducted at NASA Headquarters in Washington, D.C., December 11, 
1981. At these reviews the program participants, both graduate and under- 
graduate, have an opportunity to present their research activities to their 
peers, professors, and invited guests from government and industry. 

This conference publication represents the second in a series of yearly 
summaries of the activities of the program. (The 1980 summary appears in 
NASA CP-2176.) The majority of the material is the efforts of the stu- 
dents supported by the grants. Because of the ongoing nature of some of the 
work, certain graphics (notably photographs) are not of high quality; however, 
it was decided to publish the best material available, and if the reader's 
interest is sufficient, the appropriate advisor can be contacted for more 
recent and complete data. 

Three types of contributions are included. Completed works are repre- 
sented by the full technical papers. Research previously published in the open 
literature, for example, theses or journal articles, is presented in an anno- 
tated bibliography. Status reports of ongoing research are represented by 
copies of viewgraphs augmented with a brief descriptive text. 

Use of trade names or names of manufacturers in this report does not 
constitute an official endorsement of such products or manufacturers, either 
expressed or implied, by the National Aeronautics and Space Administration. 
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INVESTIGATION OF AIR TRANSPORTATION TECHNOLOGY AT MASSACHUSETTS 


INSTITUTE OF TECHNOLOGY, 1981 


Professor Robert W. Simpson 
Director, Flight Transportation Laboratory 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139 

SUMMARY OF RESEARCH 


There have been three areas of research sponsored by the Joint 
University Program at MIT during the past year: I) Dynamic Scheduling of 
Runway Operations at a Major Airport; 2) General Aviation Electronic Flight 
Displays - (P-POD) ; 3) Prevention of Airframe Icing with Microwaves. There 
are three faculty responsible for supervising student research in the program; 
Robert Simpson, Antonio Elias, and Walter Hollister. 

1. DYNAMIC SCHEDULING OF RUNWAY OPERATIONS 

In past years, an ATC simulation facility called TASIM (Terminal 
Area Simulation) was developed to provide a tool for investigating human 
factors problems associated with automated decision making in ATC. The 
current research goals are associated with demonstrating in real time the 
potential reductions in delay from the introduction of dynamic scheduling 
of landing and takeoff operations for a system of runways at a major airport. 
This potential was first indicated by Dear's research (ref. 1). 


In the past year two decision-making modules have been designed. 

The "Traffic Scheduler" module has been successfully coded and tested in the 
Boston scenario for a single runway with landings and takeoffs. A schedule 
of runway operations is displayed on an auxiliary CTID (Controller Traffic 
Information Display), and can be seen to change dynamically as real time 
arrivals and departures enter the system. On the controller's "radar" display, 
the current FAA ARTS III display is augmented to show "boxes" for landing 
arrivals along the extended runway centerline. These landing boxes show the 
desired position of landing aircraft for its landing schedule. A second module 
"Flight Path Generation" has been designed but is not yet tested successfully. 

It generates a sequence of commands of altitude, heading, and speed for landing 
aircraft such that they will be vectored into their box. Further modules 
associated with conformance and hazard monitoring have yet to be designed, and 
the problems of man-machine interactions with this automated "decision support" 
software which arise from incorporating controller "override" capability (so 
that he remains in command of the computer system) and from operational problems 
(such as changing runways) have not yet been addressed. There is an annotated 
slide presentation by John Pararas in this publication. 
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2 . 


GENERAL AVIATION ELECTRONIC FLIGHT DISPLAYS (P-POD) 


This research is focused on the interface between the pilot and 
the display systems of his aircraft. We are interested in the information 
handling processes and consequent pilot workloads encountered in flying 
I FR in today's and future ATC systems. In the past year, a low cost digital 
electronic flight display has been assent) led using microprocessor components 
commercially available to the home computer hobbyist. This display is called 
P-POD (Programmable Pilot-Oriented Display) and is described in a conference 
paper by Professor Antonio Elias. There are two research projects currently 
underway using P-POD. The first is described by James Littlefield in an 
annotated slide presentation given at the conference. Its goals are to create 
an operational "long range" RNAV display by interfacing P-POD to a commercially 
available Loran-C receiver. This will be tested on the bench, and in flight 
test. The second project has interfaced P-POD to a table-top general aviation 
cockpit simulator to create a novel electronic flight director display which 
provides rate of crosstrack deviation to the general aviation pilot. The 
display has other proposed projects concerned with LNAV approaches to a runway 
and the development of an "Electronic Kneepad" for flight management activities 
of the I FR pilot. These are described briefly by Elias in his paper. 

3. PREVENTION OF AIRFRAME ICING WITH MICROWAVES 

John Hansman describes his research into the possibilities for using 
microwave energy to preheat supercooled water droplets before they reach an 
airfoil. This consists of some theoretical work on droplet movement and 
flattening due to accelerations due to the flow around the airfoil. The droplet 
shape is important in determining the optimal frequency for the microwave and 
whether there is sufficient time to heat the droplets at speeds typical of 
general aviation aircraft. There is also an experimental portion of the research 
in a rain-fog section of a 1 foot by 1 foot wind tunnel at MIT. This is aimed 
at a proof of principal experiment with a prototype system mounted on an airfoil 
section in the tunnel where measurements can be made of microwave power re- 
quirements and droplet temperatures. 


REFERENCE 


1 . 


Dea r , Roger D . : 
Area . Flight 
Astronaut ics , 
Sept. 1976. 


The Dynamic Scheduling of Aircraft in the Near Terminal 
Transportation Report R76-9, Dept, of Aeronautics and 
Massachusetts Institute of Technology, Cambridge, Mass., 
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MASSACHUSETTS INSTITUTE OF TECHNOLOGY 


ANNOTATED BIBLIOGRAPHY, 1 980-81 


Pararas, John: Man Vehicle Systems Research Facility: Functional 

Specification of the ATC Subsystem. Flight Transportation Laboratory 
Report R80-13, Department of Aeronautics and Astronautics, MIT, 

December 1 9 80 . 

This document establishes the functional and operational specifications 
for the ATC subsystem of the Manned Vehicle Systems Research Facility to be built 
at NASA Ames Research Center. The ATC subsystem will provide a realistic ATC en- 
vironment for two transport cockpit simulators. Fuel mission, high fidelity simu- 
lation of extended flights in both today's and future ATC systems is required. The 
ATC subsystem design recommended consists of three audio/visual controller 
stations and a corresponding three pseudo-pilot stations with a voice disguiser 
capable of representing several different pilot voices. The ATC subsystem makes 
use of software developed for TASIM under the Joint University Program. 


Natarajan, Krishman: Use of Loran-C for General Aviation Aircraft 

Navigation. Flight Transportation Laboratory Report R81-2, 

Feb ruary 1 9 8 1 . 

This report describes an extensive evaluation of Loran-C for use by 
general aviation. Flight, ground, and antenna tests were done. Flight tests 
measured the accuracy and the ability to make approaches. Receiver reliability 
and susceptibility to atmospheric noise were also studied. Ground tests looked 
into grid stability and grid warpage. Antenna tests were done to evaluate three 
antenna configurations -- ADF, vertical whip, and trailing wire antennas. 

The measured accuracy met FAA AC 90-^5A requi rements for al 1 phases of 
flight. Loran-C was found to be satisfactory for approaches with AC 90-45A 
specifications. Reliability was 99-7%, and the receiver was insensitive to 
atmospheric noise. The time difference grid was stable in the long run. 

Antenna tests showed the ADF and vertical whip antennas to be suitable for 
ai rborne use . 

It is concluded that Loran-C is suitable for navigation as an alternative 
to VHF RNAV. This navigation system is suitable for use in general aviation 
ai rcraft . 
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DYNAMIC SCHEDULING 
OF RUNWAY OPERATIONS 


John Pararas 

MIT, Flight Transportation Laboratory 



RESEARCH OBJECTIVES 


—Demonstrate the feasibility of automated ATM/C 
decision making in the terminal area 

— Develop prototype algorithms for automation 
functions 

- Runway Scheduling 

- Flight Plan Generator 
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Terminal Area Geometry 


Due to the airway structure of the enroute airspace, landing aircraft 
tend to enter the terminal area airspace at specific points called entry 
fixes. Nominal approach routes (the dashed lines) define typical paths 
leading from the entry fixes to the runway threshold(s) . Similarly, 
nominal departure routes (dotted lines) define typical paths for departing 
ai rcraft . 

Whenever a new arrival enters the terminal area, its preferred time of 
arrival at the runway is estimated based on the entry fix and the associated 
nominal approach path. For departures the preferred time of arrival at 

the runway is determined using an estimate of the taxiing time. The 
nominal sequence of operations at the runway is determined based on the 
preferred times of arrival at the runway of all aircraft currently in 
the system. This sequence is perturbed to optimize the runway utilization 
for the given traffic, thus producing the optimal schedule of operations. 
Given the latter, ^-dimensional flight plans are generated for all aircraft. 
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Schematic Diagram of the Terminal Area ATM/C System 


The air-to-ground data acquisition (surveillance) system provides 
position information on all aircraft in the system. Nominal flight plans 
are generated for new entries and provided to the runway scheduler which 
updates the optimal schedule of operations at the runway(s). The traffic 
flight plan generator updates the flight plans according to the new 
schedule. The decision process is now complete until a new aircraft enters 
the system. 

The command processor is responsible for the timely dispatch of commands 
to the air traffic controller. This is done through the traffic display. 

The commands can also be transferred directly to the pilot via an air-to- 
ground digital data link if one is available. The commands are such that 
will keep the ai rcraft in conformance with their flight plans. The traffic 
display provides the controller with a situation display of all the traffic 
in the system and allows him to communicate with the remaining system 
through "display requests". 

The conformance monitor compares actual and desired aircraft positions 
and generates an alert when the discrepancy detected exceeds prespecified 
limits for some aircraft. This in general will cause a new flight plan 
to be generated for this (and possibly other) aircraft. The hazard 
monitor compares separations of short term projections of aircraft positions 
and generates an alert when violations of ATC separations are imminent. 

Again such an alert will require new flight plans to be generated but first 
avoidance paths are generated for the aircraft in hazard. 
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Final Vector Controller Display 


This display is designed to facilitate precise delivery of landing 
aircraft at the outer marker without automated flight plan and command 
generation. Each box corresponds to the landing aircraft having the 
same, number in the landing sequence. The boxes are moving along the runway 
centerline extension at the landing speed of the corresponding aircraft. 

The distance from the threshold is such that the aircraft reaches the runway 
at its scheduled landing time. Given the current and the landing speeds 
of the aircraft we can determine the relative bearing between the aircraft 
and its box such that if the aircraft then turns to intercept the I LS it 
will intercept the box before it reaches the outer marker. The wands shown 
are the pictorial representation of the appropriate relative bearing. 

The controller needs to approximately time the aircraft so that it will be 
within the final vectoring area when it is touched by the wand of its box. 

At that time the turn to intercept the i LS should be given. 


.ARRIVAL ON . RHS . .. 
WITH SEQUENCE NO. 


FINAL VECTORING AREA 


VECTOR WAND 
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VECTOR MARKER 








Tri-University Conference 
Accomplishments to Date 


A major portion of the work has been consumed in creating TASIM 
(Terminal Area Simulation). It is a flexible research simulator for ATC 
problems involving automated decision making for ATC, but can be useful 
to many researchers interested in ATC problems. It forms the basis for 
the ATC subsystem to be built for the MVRSF (Manned Vehicle Research 
Simulation Facility) at NASA Ames Research Center, and is of interest to 
the FAA Technical Center at Atlantic City. 

At present the simulation can be operated by a human controller through 
the keyboard in a scenario of the Boston Logan terminal area. The automated 
runway scheduler based on the CPS (Constrained Position Shifting) method 
of Dear has been coded and operates successfully with the display scheme 
shown for the Final Vectoring. The automated Flight Plan generator has 
been designed but is not yet operational. The Command Processor is 
operational and passes controller-generated commands directly to each target 
aircraft. There are two controller displays. The ATC or PVD display 
is a reproduction of the current ARTS III display (with the Final Vectoring 
display augmentat ion) . There is also an auxiliary CTID (Controller Traffic 
Information Display) which displays the dynamically changing schedule for 
runway operations, some simulation statistics, and has message areas, 
command preview areas, etc., for controller interaction. 


1. TASIM INTERACTIVE SIMULATION 

2, SOFTWARE 

RUNWAY SCHEDULER 

CPS 

TRAFFIC PLAN GENERATOR 

TERMINAL AREA NETWORK 

r- design of the algorithm 
COMMAND PROCESSOR 

CONTROLLER INPUTS 

TRAFFIC PLAN GENERATOR INPUTS 

DISPLAYS 

ATC DISPLAY (PVD) 

CTID 
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SCHEDULE INFORMATION AREA 


SIMULATION STATISTICS AREA 


AIC-JNF-QRMATION & MESSAGE / 
COMMAND PREVIEW AREA 







PPOD PROGRAMMABLE PILOT-ORIENTED DISPLAY 


Antonio L. Elias 
Assistant Professor 
Flight Transportation Laboratory 
Massachusetts Institute of Technology 
Cambridge, MA 02139 


SUMMARY 


A general-purpose low cost research microprocessor system for 
General Aviation has been developed at MIT under the sponsorship 
of the Joint University Program. This system is intended to be 
the vehicle for individual research efforts in low cost airborne 
hardware and software as well as advanced microprocessor based 
navigation systems and techniques. During 1981, two such 
research projects were undertaken, yielding results in the areas 
of micro hardware/software design, cost and performance, and 
pilot/computer interface. Five new projects are being developed 
for 1982, including a new approach to low-cost flight software 
reliability and a time-difference based Loran approach procedure 
that eliminates the need for propagation corrections and 
latitude/longitude transformations . 


INTRODUCTION 


Before the advent of microprocessors, the design of navigation 
systems was strongly influenced by the lack of convenient 
computational capabilities. This lack, combined with the need to 
minimize pilot workload, led to "direct readout" systems such as 
VOR/DME, where the signal scheme itself directly yields naviga- 
tion information in easy to use format, e.g. rho-theta. 
Similarly, systems which required any significant amount of 
"post-processing", such as Loran, were deemed unsuitable for 
aeronautical use, unless a crewmember was dedicated to the 
navigation task. 


Microprocessors make available so much computational capability 
that this traditional limitation no longer exists, thus opening 
up important new possibilities in the design and use of avionic 
systems. However, this technology has evolved in such a way that 
the principal costs involved in the use of microprocessors in 
avionic systems are not the cost of the processors itself, but of 
the environment required to support the digital processing 
functions, namely: 

1) The development (non-recurring) cost of flight 

software 

2) The certif ication of the digital hardware and 
software (also non-recurring) 

3) Repetitive design of the discrete (i.e. non-LSI) 
hardware, e.g. busses, cards, interfaces, etc. 
(mainly a recurrent cost) 
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These high non-recurrent and recurrent costs appear to have 
driven the avionics industry to a very specific technological and 
marketing niche, namely the high cost, high performance devices 
used on executive aircraft. Thus we find units such as the 
Garrett AirResearch 2000 navigator, with both high performance and 
high price. Ironically, most of these aircraft operate with two 
crewmembers so that, in spite of the increased complexity of the 
aircraft, the resulting workload may be lower than the typical 
single-pilot IFR situation on a simpler but less equipped 
aircraft. It is precisely this situation that could benefit the 
most from workload- reducing microprocessor based avionics. 


If this segment of the General Aviation community is to receive 
the benefits of microprocessor-based systems, ways must be found 
to make the low-cost market attractive to industry. The first 
step is to demonstrate that microprocessor avionics can be 
produced at a low unit cost without the risk of committing to 
very large production runs. Only after this has been 
demonstrated is industry likely to embark on larger production 
run projects with larger non-recurrent costs. 


We believe that there may be major cost-reduction demonstration 
opportunities on each of the cost areas mentioned before, in 
particular -. 


1) Improved software techniques, which may reduce both 
the development and the certification costs of 
embedded software 

2) The use of standard ("off the shelf") modules for 
commonly-used functions, such as power supply, data 
busses, interfaces, memory, processors systems, dis- 
play drivers, etc. 

3) Multifunctionality; that is, the use of a single 
digital processor to perform a number of avionic 
functions, rather than embedding individual 
microprocessors on each avionics box (1) 


(1) 


Distributed processing, that is, 
processors for each function, may 
other factors, such as performance 
over cost. 


the use of individual 
be the best design when 
and reliability, dominate 
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OVERVIEW OF PPOD PROJECT 


The motivation for developing PPOD was to make available to 
researchers a low-cost general-purpose tool with which to inves- 
tigate the areas mentioned above. The result is a 
three-processor system (main, I/O, display) using standard hard- 
ware and software, which can be configured (both from the 
hardware and the software point of view) for a wide range of 
operations and uses. The intent is that individual research 
projects be undertaken using this system, in order to try out 
novel digital avionics ideas while at the same time accumulating 
the experience of using and developing low-cost hardware and 
software. 


The PPOD hardware consists ^ of an S-100 box (including the 
motherboard and ground power supply), a single-board Z-80 
processor (including serial I/O, disk controller and EPROM 
burner), 64K of dynamic memory, 32K of ROM space, a 
high-resolution monochrome raster display (using its own Z-80), 
and a slave processor (also a Z-80) with DMA capability, 24 bits 
of parallel I/O and small amounts of RAM and ROM. 


The emphasis of all of these sub-projects will be the reduction 
of single-pilot IFR workload, with the secondary goal of reducing 
the physical complexity of cockpit instrumentation. 


SPECIFIC RESEARCH ISSUES: SOFTWARE 


Although ultimately it will be up to industry to develop the 
actual designs for low-cost microprocessor based avionics, some 
of the research that may demonstrate the viability of doing so 
may be best carried out from the detachment of a non-commercial 
research organization, such as NASA or the Universities; in 
commercial avionics development, marketing and timing reasons may 
prevent the kind of general-application research that PPOD is 
intended to support. 


In the software areas, three specific issues are addressed: 


1) The use of high-order languages 

2) The use of resident-based (as opposed to host-based) 
development systems 

3) The development of standard software functions 
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Use of high - order languages 


Large-scale aerospace processors, such as the IBM AP101, have 
modern high-order language support (e.g. HAL), whereas 
microprocessor flight software is commonly done in assembler 
code. One of the reasons for this situation is the large 
investment required to develop a high-order language environment 
(up to $1M for the HAL system) . 


On the other hand, the explosion of the personal computer market 
has resulted in the development of a number of modern high-order 
language systems (Pascal, PL1, Forth, C, etc.) for 
microprocessors. Naturally but unfortunately, most of them are 
geared towards non-real time applications, and require disk 
drives for program operation. We would like to identify exactly 
the shortcomings of these systems, and what has to be added or 
modified to make them suitable for flight software use. 


Finally, we would like to quantify the cost benefits of 
high-order languages in the GA avionics environment, in particu- 
lar the tradeoffs between the reduction in software cost, the 
increase in memory needs, and the reduction in processing 
efficiency for significantly different languages and systems. 


Alternative development systems 


There are three possible approaches to developing microprocessor 
software, listed here in order of increasing procurement cost: 


1) The use of an existing time-sharing system to edit, 
compile, link and, under certain circumstances, 
execute (by simulation) the flight code, or, 
alternatively, download it to the target machine for 
execution and checkout 

2) The use of the target microprocessor itself both for 
development and checkout 

3) The use of a specially-configured micro or 

minicomputer system geared towards software develop- 
ment (a "Development System") 


Traditionally, embedded microprocessor software has been 
developed by the first and last methods, sometimes in combina- 
tion. Use of the target machine for software development is 
another technique that has come of age with the explosion of the 
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personal computer industry. Each of the methods has advantages 
and disadvantages; with PPOD, we can experiment with the first 
two methods (the use of a Development System being beyond the 
reasonable reach of the project's budgetary expectations). 


Standard software libraries 


The last software research issue addressed is the use of standard 
functions in flight software, that is, the identification of 
functions that are general enough, and common enough so as to 
warrant the development of "library" modules to reduce subsequent 
software development/certif ication costs. Of particular interest 
is the human interface/display area, where a standardization of 
protocols and/or formats could also result in uniform user 
procedures across different designs. 


SPECIFIC RESEARCH ISSUES: HARDWARE 


In the hardware side, we will identify the hardware functions 
that are required, and relate these functions with board-level 
modules that are either currently available, or should be made 
available as "off the shelf" components for GA avionics use. 


It should be noted that the purpose is not to determine the 
suitability or unsuitability of exi sting components for GA 
avionics use, but rather to identify the functions required, and 
to quantify the benefits that would accrue from the availability 
of such standard components. As a consequence of this research, 
we hope to identify opportunities for GA-oriented industry 
standards for busses, hardware interfaces, etc. that could be an 
incentive for production of standard boards. 


Finally, in the systems design area, we hope to identify and 
quantify the acquisition, maintenance, and operating cost reduc- 
tions that would accrue from integrating, rather than 
distributing, the processing functions. 


Status of project 


As of December 1981, the following milestones had been achieved: 
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1) A hardware configuration was selected, procured and 
assembled. 

2 ) The hardware has been interfaced with a ground-based 
simulator for testing and development purposes. 

3) An operating system environment for the various 
high-order languages has been selected and procured, 
including elements of a resident software develop- 
ment system (e.g. compilers, assembler, linker, 
full-screen editor, etc.) 

4) The first high-order language system (PL1) has been 
procured. 

5) Two pilot projects have been started, and are 
yielding initial research results. 


Hardware configuration 


The "S-100" family of hardware components was selected for P-POD; 
the main basis for this selection was the large number of 
board-level components available for this family, including 
processors, memory, device controllers and I/O, ROM, video 
displays, voice I/O, etc. The bus used in this family (and from 
which the name derives) is not without its faults, particularly 
in the area of noise immunity, bandwidth, and signal path 
architecture. On the other hand, it is a flexible design whose 
noise performance has been improved with active termination and 
heavy-duty construction, and whose signal definition faults have 
been at least partially remedied by an IEEE standard definition 
(IEEE 696) of the bus. 


The initial configuration consisted of a Teletek FDC-1 
single-board processor, a Measurement Systems and Controls (MSC) 
64K dynamic memory board, a Digital Research Computers 32K 
selectable-address ROM board, and a Scion Corp. Microangelo 
display board. 


The FDC-1 board includes, in addition to the Z-80 chip, two 
serial ports, floppy controller, 2K of ROM, and an EPROM burner. 
This board is used both for software development (using the CP/M 
floppy-based operating system) and for actual execution of the 
flight code. The Microangelo display board is a 
medium-resolution (512x480) monochrome bit-mapped raster display 
with a dedicated Z-80 processor, 8K of ROM and 2K of RAM (in 
addition to the screen RAM) . Communication between the main 
processor board and the display board is by buffered interrupt 
driven physical I/O. 
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Earlier into the first project (Loran-C pilot interface) it was 
determined that the amount of I/O that could be accomplished by 
the main processor board was limited, and therefore a Teletek 12 
board was procured. This board, which is produced commercially 
as an Intelligent Interface for devices such as Winchester hard 
disks (and therefore geared towards the OEM market), consists of 
a Z-80 processor, 2K of RAM, 2K of ROM, and the bus support 
circuitry required to interface with the bus owner (the Z-80 on 
the FDC-1) by a combination of vectored interrupts and Direct 
Memo r y Access. 

These boards are contained in a California Computer Systems 
motherboard and card cage, which includes an unregulated 
heavy-duty 115 volts power supply to feed the series regulators 
on each of the boards. For flight use, the power supply is fed 
AC from a 300 watt static inverter. 


PROJECT #1: VOR/RNAV-LIKE INTERFACE FOR LORAN-C 


The first project initiated under P-POD was a display system that 
would allow a pilot to use a Loran-C receiver as a rho-theta 
navigation device similar to a VORDME/RNAV system. The motiva- 
tion for this project was the observation that pilots prefer the 
rho-theta presentation of information to almost any other type of 
presentation, including maps. Whether this preference is due to 
the considerable experience that pilots have with the VOR system, 
or whether the popularity of the VOR system is due to its 
rho-theta presentation, is subject for debate. 


In any case, the advantages of Loran (low equipment cost, wide 
area coverage) could be combined with the ease of use of VORDME 
by a device that would store names and geographical locations of 
VORTAC's. The aircraft latitude and longitude, as determined by 
Loran, is then converted to great-circle bearing and distance to 
a VORTAC selected by the pilot, and presented in a form vaguely 
resembling the familiar CDI and DME range window. In addition, 
the system allows the pilot to define (in-flight or before the 
flight) waypoints using either latitude and longitude or offsets 
from existing waypoints, and to assign three-letter identifiers 
to these new waypoints. 


To test this concept, P-POD was interfaced to an existing 
commercial receiver (a Digital Marine Northstar 6000). However, 
one of the factors that was considered in selecting this as the 
first P-POD project was the expected availability of low-cost 
Loran-C receiver boards, such as the Ohio University unit. The 
display produced, the definitions of the navigation terms used. 
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and the commands available to the user are summarized in Figures 
1 , 2 and 3 . 


The system has been tested both on the bench and on a moving 
vehicle. Flight tests during late summer of 1981 were hampered 
by power supply and antenna problems. Flight tests with a new DC 
power supply, retractable downwards-pointing whip antenna, and 
improved keyboard will begin this spring. A preliminary result 
of some interest is the figures for utilization of high-order 
language (PL1) versus assembler code, summarized in Figure 4. In 
essence, it is estimated that the use of PL1 reduced programming 
time by a factor of 7 (1000 lines of code vs. 7000 lines of 
code), while the program memory required increased by 40% (14,000 
bytes vs. 10,000 bytes). 


Of some concern is the total amount of memory used by the 
applications code (some 7700 bytes), considering that the total 
address space of these machines is only 64K bytes, and that the 
operating system takes some 6500 bytes out of them. On the other 
hand, one-third of this code is made up of pilot input handling 
and display routines that could be shared by other functions. 
Also, the system incorporates a large number of ""features" not 
all of which may be considered "indispensable". 


A major obstacle encountered during development of this code was 
that the PL1 system subroutines do not initialize storage memory, 
because they expect its contents to be initialized from diskette 
as part of the program load. When the code is resident in ROM, 
however, there is no program load, and thus no initialization of 
storage memory. A system was developed to simulate this function 
at a cost of 25% of the storage memory required by systems 
programs (application programs were coded so that they perform 
their own initialization). In spite of this fix, it was 
concluded that any high-order language system that is to be 
useful for GA avionics use must either explicitly initialize 
memory storage itself, or at least separate storage memory into 
two areas, one of which requires initialization and thus must be 
"replicated" in ROM. 


Code for this project was developed both with the resident system 
and the time-sharing system. On the basis of this experience we 
conclude that: 


1) The resident system is superior to the time-sharing 
system for developing I/O (including user interface 
and display) code, and for final code integration, 
due to its "closeness" to the target machine. 
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2) The time-sharing system can be extremely useful for 
development and coding of algorithm-intensive func- 
tions (i.e. math-intensive), due to its faster 
compilation turnaround and its more powerful mathe- 
matical debugging capabilities (debugging can be 
performed at the PL1 level). 

3) The time-sharing system is superior for maintaining 

and preprocessing data bases (e.g. VORTAC 

identifiers and positions). 


PROJECT #2: RATE -ESTIMATING ILS DISPLAY 


Flight Directors are devices that combine navigation and attitude 
information (including rate information) and generate an attitude 
command that is displayed on the same instrument that displays 
attitude. By suitable generation of the command attitude, the 
pilot workload required to fly a high-performance aircraft 
(requiring substantial command lead compensation ) can be greatly 
reduced. 


On the other hand. Flight Directors are expensive devices, 
because they require attitude and attitude rate (and, sometimes, 
control surfaces position) information, and because the computa- 
tion and display devices are electromechanical in nature. 


General aviation aircraft usually exhibit benign dynamics(l) 
which make full Flight Directors somewhat of an unnecessary 
luxury. The dynamics of lateral control during ILS precision 
approaches, however, do involve a significant amount of pilot 
workload due to the piloting strategy required by the low-gain, 
low rate information localizer display. The data processing and 
display capabilities of a P-POD like device in the cockpit could 
be used to generate a rate-aided display, maybe not as sophisti- 
cated as a Flight Director, but at essentially zero marginal cost 
given the ILS receiver and the P-POD display. 

For that cost to be indeed "zero", the display must not require 
the attitude and attitude rate data (and the corresponding 


(1) Under the sponsorship of this same NASA Joint University 
Program, Princeton is examining the dynamic and handling 
features of GA aircraft that make them more or less suitable 
for single-pilot IFR. 
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electro-mechanical instruments) that the true Flight Director 
requires . ( 1 ) 


This project attempts to produce a simple display, using ILS 
localizer information and digital filtering techniques to esti- 
mate lateral deviation rate. This rate information is then 
presented to the pilot in the form of a rotating CDI needle, very 
much like the one used in the Loran-C display. 


The specific research issues addressed by this project are: 


1) Adequacy of 8-bit processor throughput for fast 
repetition rate processes 

2) Adequacy of bit-mapped displays for fast update rate 
control oriented (as opposed to navigation oriented) 
displays 

3) Adaptability of interface subroutines developed for 
one project (Loran-C) for other projects 


This project was developed and tested using a desktop simulator 
(ATC-610) of a Piper Cherokee Arrow single-engine retractable. A 
PDP-11/10 was used to extract simulator X-Y position and synthe- 
size ILS localizer information from it. This localizer "signal" 
was then fed to P-POD, which then generated the display. Flight 
tests will begin as soon as a navigation receiver is procured for 
this purpose. 


One of the results obtained so far is that bit-mapped displays, 
such as the one used by P-POD, complicate the display design 
process because they "erase" the background over which a symbol 
is written, unless special precautions are taken to "save" the 
background information. In essence, three approaches are possi- 
ble : 


1) Design the display itself so that symbols do not 
overlap; this was the approach taken on the Loran-C 
display. In the case of the Localizer, this led to 
awkward displays. 

2) Save the pixels that are over-written by a symbol 
and restore them at the time the symbol is erased 


(1) At most, a low cost single axis fluidic rate gyro may be 
required. 


25 



for re-positioning. This requires complex graphics 
primitives and additional memory. 

3) Write the moving symbols by complementing, rather 
than 1-ing the pixels; erase the symbol by 
re-complementing, thus restoring the background sym- 
bol. Does not require complex code or memory, but 
two overlapping symbols have blank pixels at their 
common locations. 


All three techniques will be explored during the coming year. 


PROPOSED NEW PROJECTS 


A number of projects, all using the P-POD hardware and software, 
and all contributing towards the ultimate goal of identifying 
cost-reducing techniques, have been identified and will be 
started as resources become available. Some of these projects 
are .- 


1) The possible use of redundancy to improve the 
cost/reliability ratio of flight software 
("Layered-redundant coding") 

2) A Loran-C non-precision approach procedure based on 
Time Differences, not on latitude and longitude of 
runways 

3) An "electronic kneepad" that would make available 
flight information to the pilot in a dynamic fash- 
ion, as well as perform flight-plan following compu- 
tations (using pilot inputs, not direct navigation 
information) 

4) An R/Nav device using a single, digitally- tunable , 

fast-locking DME and a large VORTAC 

position/frequency catalog to achieve country wide 
"hands-off" navigation 

5) The use of voice I/O to augment/supplement the above 
functions 


The first' two projects require some elaboration: 
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There are two approaches to reliability in avionics hardware.- 
producing hardware that is essentially very reliable, and using 



less reliable hardware in redundant numbers, with the appropriate 
error detection and correction protocols. Two excellent examples 
of these approaches are the Apollo Guidance Computer, a single 
unit with MTBF's in excess of 20,000 hrs, and the Shuttle Orbiter 
Main Computer Complex, using four identical units with MTBF's of 
less than 1000 hours. Economic considerations dictate which 
approach is the most cost-effective for a particular situation. 


Software reliability can only be achieved by the first (or "brute 
force") method, since two copies of the same software are 
guaranteed to have the same identical bugs; as opposed to 
hardware, there is no distinction between the design and the 
fabrication stages, and thus there are no advantages to redundant 
copies of the same code. 


If, however, identical functional s pecifications for software are 
given to different software designers (i.e. coders), then there 
will be no correlation between the bugs in one version of the 
software and another. Such "truly redundant" software then has 
to be executed in an environment that recognizes differences in 
their outputs, and takes appropriate action in case of disagree- 
ment, much in the same way as redundant hardware. 


The cost of producing two independent codes for the same function 
is, indeed, twice that of producing a single version. (1) 
However, that cost may very well be much less than that of 
testing and inspecting the single version to achieve the same 
degree of confidence that would result from the redundant 
execution. Also, any deficiencies in the specifications them- 
selves will show up on both versions, requiring either duplicate 
specifications (hence the term "layered software redundancy"), or 
the use of conventional brute force reliability techniques for 
all higher level software and specifications. 


To our knowledge, this approach has never been tried. In order 
to evaluate its potential, an experience base must be developed, 
including cost and performance statistics. We propose to use the 
P-POD project, and the sub-projects developed under it, as a test 
bench for the redundant coding concept. Selected functions will 
be double-coded by different individuals (using different lan- 
guages, if possible) from the same specifications. The Operating 
System environment will be modified to handle the double-calling 
of the subroutines, comparing their outputs, and taking appropri- 
ate action when disagreement is detected. Finally, selected 
equivalent functions will be subject to the conventional testing 


(1) The apparent duplication of execution cost may be reduced by 
alternating codes during successive "passes", where this is 
possible . 
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and inspection approach, and records kept of the performance and 
cost of both approaches. 


The Time-difference based Loran approach procedure tries to cope 
with a major objection to the use of Loran-C for non-precision 
approaches: the use of corrections for signal propagation 
anomalies . 


Due to variations in the propagation characteristics of the 
terrain between a Loran station and a receiver position, the 
"electric" grid of hyperbolas of constant time-differences is 
"distorted" with respect to the "geometric" grid of constant 
distance-difference hyperbolas. This difference can be of the 
order of half a microsecond in time difference, or 150 meters in 
position. On the other hand, it is contended that these 
differences are repeatable within 0.1 microseconds, or some 15 
meters, leading to the suggestion that correction terms be used 
to achieve 15 meter type performance. This suggestion is viewed 
with concern by certification authorities because it may lead to 
blunders (if entered manually) or complex secular changes (simi- 
lar to magnetic variation, but infinitely more complex) . To make 
matters worse, there is disagreement on the procedures that must 
be used to convert time differences to latitude and longitude, 
and differences in the navigated position of different types of 
Loran receivers have been traced to differences in the mathemati- 
cal generation of the geometric grid. 


The proposed technique is to base the navigation solution for the 
approach entirely on time-differences, i.e. navigate on the 
electric, rather than the geometric grid. In this procedure, the 
"Loran approach" is based on five parameters: two time differ- 
ences at each end of the runway, and the distance between these 
points, namely the runway length. (1) The benefits of this 
technique are twofold: a reduction of the complications (and 
potential errors) resulting from the use of propagation correc- 
tions and various lat-lon conversion schemes, and a reduction of 
the complexity of the navigation computations required during 
approach, which can lead to faster display update rates. 


(1) Additionally, some approximate information about the location 
of the airport with respect to the stations may be used to 
improve the linearity of the approximations involved, partic- 
ularly at airports situated close to a Loran station; but in 
any case the navigation solution guides to time-differences , 
not to a point in lat and Ion. 
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Figure 1. Loran-C navigator commands. 
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Figure 2. Enroute Loran navigation - geometry and definitions. 






E — enter new waypoint in waypoint list 
I — insert new waypoint in waypoint list 
D — delete a waypoint 

A — advance the "current waypoint" to the next 
waypoint in the list 

B — back up the "current waypoint" to the previous 
waypoint in the list 

S — show a waypoint (from the waypoint list) in 
the waypoint display area 

T — set CDI tracking to the "TO" position 

F — set CDI tracking to the "FROM" position 

R — set desired tracking radial to/from a waypoint 


Figure 3- Summary of Loran-C navigator commands. 


PROGRAM MEMORY 


ITEM 

SOURCE LINES 

CODE (BYTES) 


PL1 CODE 

706 

7177 

(10 BYTES/LINE) 

ASSEMBLER CODE 

365 

511 

(1.4 BYTES/LINE) 

SYSTEM CODE: 

_ 

6522 



TOTAL PROGRAM MEMORY: 14210 

ESTIMATED SAVINGS AND PENALTIES 
( PL1 VS. ASSEMBLER) 

- MEMORY PENALTY: 4000 BYTES (10000 VS 14000) 

- SOURCE CODE SAVINGS: 6000 LINES (1000 VS 7000) 

Figure 4. Core and memory usage summary - enroute Loran navigation 

pilot interface software. 
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THE P-POD PROJECT: PROGRESS REPORT 


James A. Littlefield 
MIT, Flight Transportation Laboratory 



P-POD PROJECT: PROGRESS REPORT 


The programmable pilot oriented display (P-POD) is a multiprocessor 
based flight information processing and display unit. Host of the work 
to date has been directed toward optimizing communication protocols 
between the three Z-80 processors inside P-POD. As described in a previous 
report, one processor (PO) is dedicated to 10 control, data buffering, 
and formatting; a second processor (Pi) is reserved for coordinate trans- 
formations and supervision of the other two Z-80's; the third unit (P2) 
is dedicated solely to display generation. 

At the last quarterly meeting of the Tri -Uni vers i ty program a 
description of the interface between the video processor, P2 , and the 
main processor, PI, was given. Subsequent reliability testing led to 
significant revisions in both hardware and software. Figure 1 illustrates 
the cycle of data flow between PI and P2 . Should a failure (either 
hardware or software timing) occur at any point in the data transmission 
cycle, the process will hang indefinitely, resulting in cessation of 
display generation and an overflow of the P1/P2 command buffer. While 
running diagnostic routines to exercise the P1/P2 interface at maximum 
speed an intermittent failure was observed. A careful check of possible 
noise sources or race conditions in the hardware diagramed in figure 2 
did not reveal the source of this failure. Possible software timing 

problems were also fully investigated. Due to the sporadic nature of 
this failure it was not possible to isolate its source; however, all 
test results point to a portion of the status generation circuit on the 
video board as the error source. 

To eliminate the effects of this type of failure a performance 
monitor routine is now being used to detect these occasional failures 
and co r re ct fo r th em . Th is s of twa re co rrection for har dwa re failure was 

deemed prudent for two reasons; detailed debugging of the video board 
would be very time-consuming with the available documentation; these failures 
are serious and might be induced by external noise sources. At present. 


32 



the real time clock on board PI is used to trigger execution of a routine 
which observes the history of interface operation over the previous second. 
Based on the setting of a flag and contents of a status byte, the data 
transmission cycle between P1/P2 is either allowed to continue normal 

operation or is re- i ni ti al i zed and restarted when an error condition is 

detected . 

One of the advantages of a software reliability monitor is that it 

will correct for any temporary operational failure in the interface (i.e., 

power line spike, missed interrupt, software timing error, EMl). 

In addition the assembly coded monitor routine occupies less than 41 bytes. 

The performance check routine runs every second in the current implementation. 
Since all commands to the graphics processor, P2, are buffered, no bytes are 
lost due to momentary hardware failures of the type observed. Command 
buffer contents are preserved until the transmission cycle is restarted by 
error detection code. An additional function of the performance monitor 
is accumulation of a total error count indicating the number of failures 
detected and corrected since system startup. Presently failure rates of 
approximately 6/(3 x 10^) bytes are being observed. 

Loran data processing has been the second major area of work during 
the past quarter. Two 8-bit latches have been connected to the Loran 
receiver data bus. Whenever new data is latched into these buffers a 
mode 1 Z-80 interrupt is sent to P0. P0 multiplexes the 16 bits of data 

onto its own 8 bit 10 port by cycling the tri-state enable lines of the 
octal latches. Once Loran data has entered the P0 memory map it is 
processed by the control program outlined in figure 3. 

The Northstar 6000 receiver is capable of transmitting digital data 
to four external devices (D0-D3) • Data for each of these devices is 
multiplexed onto the same data bus. IUSER3, the P0 control program, 
decodes 2 of the four available devices. D1 data consists of TD's for 
the first and second slaves in the selected chain. D2 data is composed 
of TD's, SNR's, ECD's, mode numbers, and latitude/longitude for all 
stations in the chain. D1 data is available once every 2 seconds. 
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D2 data is available about every 2.7 seconds. iUSER3 will observe both 
D1 , D2 data and transmit either or both to PI via a DMA operation. 

IUSER3 can also be configured to select the best data source; if the 
first and second slaves have high enough SNR's, then D1 data is preferred 
since it is available more often. SNR's are determined by observing 
the D2 data stream at all times. Should first or second slave SNR's 

fall below a programmed threshold value, the D2 data stream will be sent 
to PI. Whenever data is transferred to PI memory, the first byte of 
the transmit buffer contains a tag indicating the contents of that 
buffer ( D 1 , D2 data). 

Use of the additional 10 processor frees the PI processor from 
several tasks. The P0 control program formats each data frame and does 
length checking to ensure that each frame is complete. Thus PI is assured 
of having a valid position fix before attempting coordinate conversion. 

Objectives for the next quarter are completion of the interface 
between the P0 control program and the existing PL/1 Loran data processing 
code which was demonstrated at Ohio University. Successful flight 
testing of this configuration will be followed by some execution speed 
tests to more fully document the performance of the interprocessor 
communications package. Final documentation will complete my thesis 
work. The interprocessor communications package will provide subsequent 
users of PPOD with a modular operating environment which will ease adaptation 
to a wide range of experimental goals. 
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Figure 2. P1/P2 communication data flow 


















Figure 3. (Concluded) 












MICROWAVE ICE PREVENTION 


R. John Hansman, Jr. 
and 
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MIT, Flight Transportation Laboratory 



RESEARCH ON ANTI-ICING BY MICROWAVE HEATING OF SUPERCOOLED WATER DROPLETS 


The study of the novel concept of using microwave energy to 
provide aircraft ice protection, specifically an anti -icing system, 
and the feasibility of such a system is the aim of this research. 

In a microwave anti-icing system impinging supercooled water droplets 
are heated to above freezing by the resonant absorption of microwave 
energy located upstream of the airfoil. This process is inherently 
more efficient than existing anti-icing devices due to the saving of 
the latent heat of fusion (a substantial 334 joules/gm (80 cal/gm)) 
and the fact that only the droplets are heated, thereby reducing 
convective losses to the air. 

Efficiency improvements of factors of 10 or 20 are, in principle, 
possible. This would make anti-icing feasible in applications where 
less desirable de-icing devices are now used. The light weight and 
low aerodynamic drag of microwave systems are also desirable characteristi 

In studying the microwave anti-ice concept three major issues 
became obvious: 

1. Is there sufficient time to heat the supercooled water 
droplets? (i.e., will it work?) 

2. What efficiencies are realistically possible? 

3. What are the applications problems? (i.e., cost, high 
frequency, i nterference , runback, etc.) 

The work to date has been focused on the first two questions. 

The question of droplet heating time has been approached on 
both a theoretical and experimental basis. The theoretical work 
has resulted in a model for the heating of a droplet as it approaches 
an airfoil. The model has undergone several levels of refinement 
and is manifest in a computer code which calculates droplet trajectories 
and temperatures. At present, the model includes the following 



effects or assumptions: 

Rayleigh-Gans absorption 
Inviscid airflow around body 

Electromagnetic field (can be varied — normally exponential radial , 

2 

cos angular behavior) 

Droplet trajectories -- iterative solution of equations of motion 

from reference 1 
Increased heating due to: 

droplets slowing before impact 

increased absorption from flattening due to pressure gradients 
Polarization and near field effects 
Evaporative and convective losses 
Max. electric field limited by gaseous breakdown 

The code is 75% complete and when completed should provide 
detailed design and efficiency data. Hand calculations of droplet 
heating including the aboVe effects indicate that there is clearly 
time to heat the droplets if the microwave frequency is sufficiently 
high. Pending completion of the code calculations the operating 
frequency is expected to be between 10 GHz and 30 GHz. 

One effect which is considered important but was not included 
in the above computer simulation due to a lack of a quantative model 
is the increase in droplet absorption as it passes through the ice-water 
phase transition. An increase in absorption is expected at the phase 
transition due to the nonspherical shapes of freezing and melting 
droplets . 

The theoretical work described above is being supprted by an 
experimental program. The predictions of droplet trajectories and 
droplet flattening are being compared with wind tunnel investigations 
of droplet flow around cylindrical and airfoil shapes. Velocity and 
flattening information are obtained from fast strobe photography and 
fixed "time of flight" photography in a rain-fog section of the MIT 



1 ft by 1 ft tunnel. In addition an experiment is being set up to 
quantify droplet absorption through the ice-water phase transition. 

This sensitive measurement is made by observing the change in Q 
(Q = power stored/power dissipated) of a resonant cavity which 
contains melting or freezing water droplets. For more details of 
the experimental or theoretical work to date, see figures I through 
17, which were presented at the NASA Joint University Air Trans- 
portation Research Conference. 

With the completion of the computer simulation of droplet 
heating and its experimental verification the question of whether 
there is sufficient heating time will be answered, and we are therefore 
prepared to begin the next phase of the work. This phase consists 
of two parallel tasks. 

(1) Proof of principal experiment, consisting of the design of 
a microwave anti-icing system using the above analysis. The system 
will be mounted on an airfoil section and tested in the wind tunnel. 

The system can be tested both at room temperature and at icing temperatures 
(the winter operation of MIT 1 x 1 tunnel is adequate for this purpose). 
Measurements can be made of microwave power drain and of droplet 
temperature. The temperature measurement can be made in situ by 
using advanced fiberoptic thermometry which allows an unperturbing 
measurement of temperature even in the presence of microwave fields. 

(2) Completion of droplet heating work by measuring absorption 
during the ice-water phase transition described earlier. This work 
should provide some insight into both the "mixed phase" icing conditions 
about which there is only limited information and the phenomenon of 

the "bright band" . 
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CLOUD DROPLETS (FAR 25) 


FREEZING RAIN 


T 0 -¥ -20°C 0 -5°C 

LWC 0 -*■ 1 gm/m 3 0 -► 3 gm/m 3 

D 10 40 fim 1 -► 5 mm 

Generally below 10,000 ft in localized layers and zones 
PROBLEM BECAUSE OF: 

weight ice-ingestion 

drag antenna 

flutter (anti-ice) 

CURRENT TECHNIQUES 

boots hot air 

electrothermal liquid-glycol or alcohol, TKS 

Figure 1. The general parameters of the icing problem and 
the current techniques used to deal with it. 

1. Ice prevention — use resonant absorption properties 
of water 

A. Anti-Ice (Preheat Droplets) 

B, De-Ice (Melt Ice) 

2. The possibility of using radar to remotely detect 
and avoid icing zones was first proposed by 
David Atlas in 1954 J. Met. Vol. 11, p. 309. While 
his concept was based on a ground weather radar, 
advances in "on board" radar capability along 

with certain scattering characteristics of droplets 
at the ice-water phase transition (to be discussed) 
may make remote detection feasible. Requires 
vertical scan; 

Figure 2. Several concepts of using radiation for ice protection or remote 
detection. The idea of remotely detecting icing zones by air or ground 
radar is intriguing and was spawned by meteorological radar observations 
of a highly-reflecting "bright-band" at the melting layer in precipitation. 
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This scheme operates in two modes: anti-ice or de-ice 


ANTI-ICE 

DroDlets are preheated to above freezing before impact 
Power requirements estimated (V - 200 mph, 20 ;jm droplets) 

Wing 40 watts/linear foot 

Prop or Rotor 25 watts/llnear foot (Goodrich "Hot Prop", 

200 watts/ linear foot) 

DE-ICE (Hagenheim) 

Ice on surface waveguide causes the effective thickness of the 
dielectric to increase which causes the wave to be more tightly 
bound. The high fields which result cause heating at the ice- 
dielectric interface, thus unbinding the ice from the surface. 

Figure 3. One possible scheme for using radiation for 
ice prevention, primarily anti-icing. First-order 
conceptual scheme employing surface waveguides. 


(1) Only water heated — low convective losses 

(2) Saves latent heat 

L iw “ gm" L wv = 800 

(3) Only draws power when liquid water present 
( i . e . , serves as its own detector) 

(4) Aerodynamical ly clean 

(5) Low maintenance 

(6) High power microwave technology on o strong 
growth trend 

Figure 4. Several potential advantages 
of a microwave anti-icing system. 
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RATIO OF ENERGY REQUIRED 
TO HEAT WATER COMPARED T: 

THE ENERGY REQUIRED TO HEAT 
AND KLT ICE 

joules/gn. (60 cal/grr.) 
C^-4. 19 joules/cr °C (1 cal/gir. c C) 




TB 2o S Jo 35 5fT 

41 (°C) 


Figure 5. Efficiency improvement gained by avoiding the loss of the latent heat 
of fusion at the ice-water phase transition (80 ml/gm). Note that a similar but 
greater effect occurs with the latent heat of evaporation (600 cal/gm). 

1. Is there sufficient time to heat the droplets? i.e., what are 
the absorption properties of the impinging droplets? 



efficiency goes as h^/h r j implies shorter wavelength desirable 


h^ -w 1 cm 


3. Runback refreeze problems 


4 . Applications problems (i.e., cost, HFI etc) 


While all the above are important, the most germaine is 11. 

Figure 6. Major questions which must be addressed in order to determine the 
feasibility and operating frequency of a microwave anti-ice system. 
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initial position x Q and velocity v fl 
droplet size D 





l 

( y es ) 


I 

output trajectories 
and temperatures 


Figure 7. Flow chart for a computer simulation of droplet trajectories and heating 
This code has been written to help address the questions of droplet heating time 
and efficiency. The absorption cross section includes Rayleigh or Mie absorption 
with pressure-induced shape effects, polarization, and phase effects included. 
Equation of motion is from reference 1. 
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Figure 8. Computer-calculated trajectories are checked by computing the terminal 
fall velocity of drops (m/sec) versus droplet diameter (mm). With the drag 
coefficients modified to include droplet flattening the code is in good agreement 
with established experimental data. 
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Y AXIS 


Figure 9. Droplet trajectories near half body. Droplets of four 
different sizes, 5, 10, 20, and 40 pm, approach a "half-body" 
airfoil. The smaller droplets are accelerated around the body 
while the larger droplets impact. The effect of the increasing 
catch efficiency with droplet size is discussed in figures 10 
and 1 1 . 


Fai r 

Weather 

Cumulus 


Cumulus / 
Conjestus^ 

1/ 


Stratus 


/. 

/ I 

// 

/ 1 
/ 

/ 

/ 

/ 


d (microns) 


Figure 10. Fraction of liquid water with diameter smaller than d, 
showing typical droplet size distributions in several different 
types of clouds. 


KHRGIAN-MAZAN DISTRIBUTION (Arb. Units) 


1. Drop size in cloud average diameter 



Figure 11. Illustration of the effect of the larger droplets being 
caught more efficiently by the airfoil. The larger droplets 
contribute appreciably more to the icing mass than the more 
prevalent smaller droplets. 


Current anti-ice systems evaporate all water 

L wv = 600 cal/gm C w = 1 col/gm K 


Preliminary heat transfer calculations over wet wing show factor 
of 10 to 20 improvement over evaporative systems 
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Figure 12. Some preliminary calculations of heat loss during runback. It is clear 
that the coolest feasible surface temperature is the most desirable for an elec- 
trical system, as opposed to the 20° to 30°C temperatures used in hot air systems. 
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Figure 13. Photographic setup used to attempt to validate the 
computer simulation of droplet flattening and trajectories. 



Figure 14. Photographs taken on the stagnation streamline of a 4. 
diameter cylinder with a free-stream velocity of 60 mph. 








Figure 15. Droplet velocity measurement by double strobe photography. 



d (mm. ) 

Figure 16. Plot of velocity along tunnel axis (1.25 m downstream of injection site) 
versus droplet diameter for droplets injected into an unobstructed tunnel. These 
velocity measurements were made by the technique described previously and are in 
good agreement with the computer simulation values (dashed lines). 
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Figure 17. Schematic of phase-change absorption experiment 
being done to quantify absorption of droplets as they 
pass through ice-water phase transition. 
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INVESTIGATION OF AIR TRANSPORTATION TECHNOLOGY 
AT OHIO UNIVERSITY, 1981 


Professor Richard H. McFarland 
Avionics Engineering Center 
Department of Electrical Engineering 
Ohio University 
Athens, Ohio 45701 


INTRODUCTORY REMARKS 

Loran-C has had much of its development motivated by military and 
marine applications. The increased availability of Loran signals In the 
United States has encouraged consideration of Loran for airborne applica- 
tions and has motivated a somewhat similar evaluation and approach as was 
performed by Ohio University several years ago in this NASA-supported joint 
university program then using Omega navigation. 

The low-frequency band because of its relatively high noise contami- 
nation, requires high-quality signal processing to obtain effective signal- 
to-nolse ratios which permit good reliability in position determination 
consistent with airborne applications. To this end Ohio University has 
investigated techniques for deriving air-navigation quality information 
from Loran-C. A part of this objective has been to produce the navigation 
data in as efficient and cost-effective a manner as possible. Specifically, 
the volume, weight and power requirements of the equipment are to be com- 
patible with light airplane flown in the general aviation environment. 

The Lilley— McCall paper describes the effort which led to a prototype 
Loran-C receiver which was designed, built and flight tested as a part of 
this year's work effort. Two undergraduate students, two graduate students 
and two staff supervisors worked on this project for over a year. The 
final results have been to produce a receiver which, for demonstrations, 
operates with a master-dependent, single chain tracking three stations. 

The Burhans paper which follows reveals the development of the RF cir- 
cuitry. During the evolution of this receiver prototype, students in par- 
ticular have also obtained rich experiences in microprocessor applications, 
software development, fabrication techniques and final proofing of concepts 
by measurements of navigation performance in flight on board a Piper 
Cherokee and a Douglas DC-3. 

The use of a navigation receiver during flight leads to consideration of 
several ancillary aspects of the navigation equipment itself. One con- 
sideration is the development of a more intuitively satisfying and useful 
reference coordinate frame. The hyperbolic geometries inherent in the 
Loran time-difference system are not easily used by pilots; consequently, a 
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coordinate conversion is highly desirable* One approach to a coordinate 
conversion uses an on-board microprocessor for producing the convenient 
and common polar coordinate geometries available with contemporary VORTAC 
navigation. 

Finally, and with work in progress, the concern is with the human 
interface. Display of the navigation information in a manner which is con- 
sistent with the general aviation cockpit environment and in a form which 
is convenient and conducive for interpretation by the pilot is the goal. 
Completion of the display work described by Novacki will mean that the pro- 
ject has been involved with a totality of aspects of the Loran-C airborne 
equipment ranging from conceptualizing and design to flight testing. 

Further movement along this line will take place with colleagues at MIT 
obtaining data on performance and pilot usage of the equipment. 

Following also are a bibliography of reports and presentations of 
reports completed during 1981 on the tasks integral to the Loran airborne 
receiver developed at Ohio University. 

The reports referred to as OU NASA TM-XX in the papers and in the annotated 
bibliography are the work of Ohio University and should not be confused with 
the NASA TM series of reports, which are represented by five-digit numbers 
(NASA TM-XXXXX) . Research done at Ohio University under the Joint University 
Program is supported by NASA and reported in the OU NASA TM-XX series. These 
reports may be obtained from the University directly. 



FIRST TECHNICAL AREA - OMEGA NAVIGATION 


Develop Technology to Make Omega Suitable for General Aviation 
® RF Circuit Design 

• Signal Processing 

• Receiver Design 

• Flight Evaluation 

IDENTIFY FUNDAMENTAL PROBLEMS AFFECTING GENERAL AVIATION 

• Noise Factors 

• Geometries 

® Propagation factors 

• Accuracy requirements 

SECOND TECHNICAL AREA - LORAN-C 

Develop technology making Loran-C applicable to general aviation 

• RF circuits 

• AGC 

• Envelope detection 

• Prototype Receiver 

Investigations of noise effects 

Phase-lock tracking loops 

Applications of microprocessor techniques to Loran-C 
Coordinate Conversions 
Display Interface 
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BENEFITS 

Technology product - 82 Technical Memorandums 

Student product - Over 31 graduates 
(BSEE, MS and PhD levels) 

Student improvement in the lab and classroom 

Industry benefits - technical information provided to more than 
100 contacts- Examples ITT, Litton, Teledyne, 

Magnavox, Rockwell 
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OU NASA TM Number 


ANNOTATED BIBLIOGRAPHY 


1 SIMULTANEOUS PAIR OMEGA RECEIVER 
Ralph W. Burhans, August 4, 1972. 

A new concept of OMEGA receiver operation is presented. The 
simultaneous comparison of 13.6 and 10.2 KHz in a single time 
slot results in 3.4 KHz, 10.2 KHz, and phantom 40.8 KHz dif- 
ference pair lanes. Direct lane count without intermediate 
storage or sequential feedback loops is possible. Small boat 
receivers using low cost recorders such as a Rustrak might be 
possible at $300 market prices. Completely digital readout 
devices for general aviation use could become $500 instru- 
ments. Even with restricted lane pair use tills could provide 
low cost navigation aids for over half of the U.S.A. and all 
of the coastal areas including most of the Great Lakes. 

2 SIMPLE BAND PASS FILTERS 

Ralph W. Burhans, December 28, 1972. 

A single power source operational amplifier provides tunable- 
high Q band pass resonators of possible utility for OMEGA VLF 
receiver input processors and other audio frequency applica- 
tions. 

3 LOW BIT SINE WAVE APPROXIMATIONS FOR AUDIO SIGNAL SOURCES 
Ralph W. Burhans, January 1973. 

Special purpose sine wave sources are easily obtained using 
digital counting-decoder-filter methods with second harmonic 
distortion less than 1%. Good frequency stability and reason- 
able variety in choice of output frequency result when the 
input clock signal is derived from a crystal oscillator 
through programmable divider chains. The resulting sine wave 
output would be suitable for applications such as the 90 and 
150 Hz modulation frequencies of ILS transmitters, as well as 
other limited range audio frequency use. 

4 SIMPLIFIED OMEGA RECEIVERS 
Ralph W. Burhans, March 1974. 

Circuit details are presented for a low cost OMEGA receiver 
being developed for general aviation use. Some novel pro- 
cessing methods, not used in commercial systems, have been 
demonstrated in experimental bench processors. An airborne 
model is being designed. 

5 BINARY PHASE LOCK LOOPS FOR SIMPLIFIED OMEGA RECEIVERS 
Ralph W. Burhans, March 1974. 

A sample binary phase lock loop is proposed for periodically 
correcting OMEGA receiver internal clocks. The circuit is 
particularly simple to implement and provides a means of 
generating long range 3.4 KHz difference frequency lanes from 
simultaneous pair measurements. 
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PHASE LOCK LOOP SYNTHESIZER FOR OMEGA REFERENCE FREQUENCIES 
Kent A. Chamberlin, April 1974. 

An OMEGA reference frequency of 4080 KHz is provided by a 
single loop VCXO circuit driven from an atomic clock or stable 
crystal standard. The circuit is used to provide OMEGA 
frequencies for direct ranging, differential correction moni- 
toring, or as a laboratory source for calibrating OMEGA 
receivers. 

SIMULTANEOUS MASTER-SLAVE OMEGA PAIRS 
Ralph W. Burhans, April 1974. 

Master-Slave sequence ordering of the OMEGA system Is 
suggested as a method of improving the pair geometry for low 
cost receiver user benefit. The sequence change will not 
affect present sophisticated processor users other than to 
require new labels for some pair combinations, but may require 
worldwide transmitter operators to slightly alter their long 
range synchronizing techniques. 

SELECTED BIBLIOGRAPHY OF OMEGA, VLF AND LF TECHNIQUES APPLIED 
TO AIRCRAFT NAVIGATION SYSTEMS 
NASA Project Staff, August 1974. 

A bibliography of references collected during the first three 
years of the NASA Tri-University Program in Air Transportation 
Systems. 

LOW-COST OMEGA NAVIGATION RECEIVER 
Robert W. Lilley, October 1974. 

The status of Ohio University's efforts towards specifying a 
low-cost Omega receiver is reviewed at the onset of the 
fourth-year program under the NASA Tri-University Program in 
Air Transportation Systems. 

BINARY PROCESSING CONCEPTS FOR OMEGA RECEIVERS 
Robert W. Lilley, November 1974. 

Preprint of paper presented at the Second Omega Symposium, 
sponsored by the Institute of Navigation, Washington, D.C., 
November 7, 1974. 

THE MEMORY-AIDED DIGITAL PHASE-LOCKED LOOP 
Kent A. Chamberlin, November 1974. 

Preprint of paper presented at the Second Omega Symposium, 
sponsored by the Institute of Navigation, Washington, D.C., 
November 7, 1974. 



OMEGA FLIGHT-TEST DATA REDUCTION SEQUENCE 
Robert W. Lilley, November 1974. 

A series of FORTRAN computer programs for preparation and sum- 
mary of flight-test data obtained from the Ohio University 
Omega Receiver. 

FLIGHT EVALUATION: OHIO UNIVERSITY OMEGA RECEIVER BASE 
Kent A. Chamberlin, R. W. Lilley, and Richard J. Salter, 
November 1974. 

A description is given of the data-collection flight, round- 
trip from Athens, Ohio to Langley Field, Virginia, during 
which Omega data was collected on machine-readable media for 
use in the Tri-University Program in Air Transportation. 

COMPUTER PROGRAM CORDET 

R. A. Palkovic, November 1974. 

A simulation tool is described for use in the design and ana- 
lysis of digital phase-locked loops, with specific application 
to the DPLL in the Ohio University Omega receiver base. 

A SIMULATION ANALYSIS OF PHASE PROCESSING CIRCUITRY IN THE 
OHIO UNIVERSITY OMEGA RECEIVER PROTOTYPE (Master's Thesis) 

R. A. Palkovic, June 1975 (CR-132707). 

A first-order digital phase-lock loop is modeled on the com- 
puter. Loop response to signal phase in noise is evaluated. 
Optimum integration time is determined. Phase jitter in a 
frequency synthesizer used as the local oscillator is quan- 
tified, and design is optimized. Design rules for use of 
synchronous rate multipliers are presented. Overall system 
response is discussed. 

GANGED SERIES POTENTIOMETER MIXER NETWORKS 
Ralph W. Burhans, December 1974. 

A ganged potentiometer with a single linear section and two 
opposite log tapered sections is rediscovered for providing a 
simple series resistor control element for mixing of audio 
frequency signals. An application is for bench evaluation of 
detector signal-to-noise ratios with Omega receivers. 

COMMON ANTENNA PREAMPLIFIER-ISOLATOR FOR VLF-LF RECEIVERS 
Ralph W. Burhans, July 1975. 

An improved high impedance preamplifier circuit provides out- 
puts to drive an Omega-VLF receiver and an ADF-LF receiver 
from a common antenna such as the ADF sense antenna on general 
aviation aircraft. The preamplifier has been evaluated with 
fixed ground station receivers and is anticipated for use in 
the second generation prototype Ohio University Omega 
receiver design. 


LOW COST, HIGH-PERFORMANCE, VLF RECEIVER FRONT-END 
Ralph W. Burhans, September 1975. 

A VLF receiver front-end has been designed using standard 
linear integrated circuits. The basic methods have been 
evaluated extensively on the Omega 10.2 KHz channel but are 
readily adaptable to any other VLF frequency in the 10.2 KHz 
to 20 KHz region. Applications for the modules exist in 
position location, time-frequency measurements, and signal 
propagation. The set provides control gates, zero crossing 
signals, and analog outputs to interface with any type of 
digital logic, microprocessor, or analog signal processor. 

DIGITAL CORRELATION DETECTOR FOR LOW-COST OMEGA NAVIGATION 
Kent A. Chamberlin, February 1976 (CR-144956) (Master's Thesis) 

This report describes the background information on the 
research that led to the development of the memory— aided 
phase-locked loop (MAPLL) which is an all-digital correlation 
device that is capable of determining the phase of extremely 
noisy fixed-frequency signals. This design is of special 
interest for Omega or other phase sensitive VLF navigation 
purposes since it is relatively inexpensive, maintenance-free, 
and can operate in a time-multiplexed fashion. 

THE MINI-0, A DIGITAL SUPERHET, OR A TRULY LOW-COST OMEGA 
NAVIGATION RECEIVER 

Ralph W. Burhans, November 1975. (CR-144923) 

A quartz tuning fork filter circuit and some unique CMOS clock 
logic methods provide a very simple OMEGA-VLF receiver with 
true hyperbolic station pair phase difference outputs. An 
experimental system has been Implemented on a single battery- 
operated circuit board requiring only an external antenna 
preamplifier and LOP output recorder. A bench evaluation and 
preliminary navigation tests indicate the technique is viable 
and can provide very low-cost OMEGA measurement systems. The 
method is promising for marine use with small boats in the 
present form, but might be implemented in conjunction with 
digital microprocessors for airborne navigation aids. 

FLIGHT TEST OF 4-HZ AND 30-HZ OMEGA RECEIVER FRONT-END 
Lee Wright, February 1976. 

A test flight in DC-3 aircraft was conducted to evaluate the 
performance of a 4-Hz ultra-narrowband, Omega receiver front- 
end compared to a more conventional 30-Hz bandwidth receiver. 
Results indicate that the 4-Hz front-end has superior signal- 
to-noise performance. Other interesting results obtained 
during the test flight were recordings of the sunset noise 
effects on amplitude, and the attenuation of signal levels 
when flying through clouds. 



22 POSSIBLE METHODS FOR USSR-VLF NAVIGATION RECEIVERS 
Ralph W. Burhans, March 1976. 

A brief study of the USSR-VLF navigation system indicates that 
very low-cost digital techniques might be applied to receiver 
systems. The transmitted signal format is of interest for 
application to other VLF systems in the future. Some possible 
circuits for simplified receiver processors are presented. 

23 A PROPOSED MICROCOMPUTER IMPLEMENTATION OF AN OMEGA NAVIGATION 
PROCESSOR 

John D. Abel, March 1976. 

Documentation of current status of research pertaining to a 
microprocessor-based Omega navigation processor to be used in 
conjunction with the Ohio University Avionics Engineering 
Center Omega sensor processor is presented. 

24 IMPROVEMENTS FOR OMEGA RF PREAMPLIFIERS 
Lee Wright, April 1976 

An Omega preamplifier with no phase shift over the ADF band 
but with bandpass filtering and gain at the Omega-VLF band has 
been designed, built, and tested. This is expected to be use- 
ful principally in planned work at MIT and Princeton Involving 
the use of the Ohio University Omega Sensor Processor 
Receiver. 

25 NARROW BAND BINARY PHASE LOCKED LOOPS 
Ralph W. Burhans, April 1976. 

Very high Q digital filtering circuits for audio frequencies 
in the range of 1 Hz to 15 KHz are implemented in simple CMOS 
hardware using a binary local reference clock frequency. The 
circuits have application to VLF navigation receivers and 
other narrow band audio range tracking problems. 

26 SIMULATION ANALYSIS OF A MICROCOMPUTER-BASED LOW-COST OMEGA 
NAVIGATION SYSTEM 

Robert W. Lilley and Richard J. Salter, Jr., May 1976. 

Preprint of paper presented at the Bicentennial National 
Symposium of the Institute of Navigation, Warminster, 
Pennsylvania, April 28, 1976. 

27 AUTOMATIC NOISE LIMITER-BLANKER 
Ralph. W. Burhans, May 1976. 

Modifications of an audio noise limiter circuit, used in WW II 
era radio communications receivers, provides a noise limiter- 
blanker for narrow bandwidth low-level audio signals. The 
method has been evaluated for noise blanking with OMEGA-VLF 
navigation receivers but is adaptable to more general audio 
frequency processing systems. 
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SMALL AIRCRAFT FLIGHT EVALUATION OF RUSTRAK CHART RECORDER 
Richard J. Salter, Jr. and Robert W. Lilley, May 1976. 

In support of the NASA Omega Prototype Receiver project, three 
short flight evaluations of the RUSTRAK chart recorder were 
flown. 

DIGITAL TIME SLOT DISPLAY FOR OMEGA RECEIVER 
Ralph E. Smith, July 1976. 

Variations of circuits to display the Omega sequence letters 
A through H have been designed, breadboarded, and tested. 

One of these is suggested as an alternative station display 
method for the Ohio University Omega sensor processor systems. 

LOW-COST MECHANICAL FILTERS FOR OMEGA RECEIVERS 
R. W. Burhans, June 1976. 

A set of mechanical filter assemblies has been obtained for 
possible use in the RF front-end of an OMEGA navigation 
receiver. The resonators provide very narrow bandwidth per- 
formance with good skirt selectivity in a simple two-stage 
circuit. It is recommended that these filters be used in a 
complete receiver system for long-term evaluation on reception 
of low-level OMEGA signals. 

A MEMORY -MAPPED OUTPUT INTERFACE: OMEGA NAVIGATION OUTPUT 
DATA FROM THE JOLT (TM) MICROCOMPUTER 
R. W. Lilley, August 1976. 

A hardware interface which allows both digital and analog data 
output from the JOLT microcomputer is described in context with 
the Ohio University software-based Omega Navigation Receiver. 

A MICROPROCESSOR INTERFACE FOR THE OHIO UNIVERSITY PROTOTYPE 
OMEGA NAVIGATION RECEIVER 
R. W. Lilley, August 1976. 

A hardware interface is described which allows a microcomputer 
to obtain data and interrupt signals from the Ohio University 
Omega Receiver Prototype. 

TEST PROGRAM FOR 4-K MEMORY CARD, JOLT MICROPROCESSOR 
R. W. Lilley, August 1976. 

A memory test program is described for use with the JOLT 
micro-computer memory board used in development of the Ohio 
University Omega navigation receiver. 

A MICROCOMPUTER-BASED LOW-COST OMEGA NAVIGATION SYSTEM 
R. W. Lilley and R. J. Salter, Jr., August 1976. 

Preprint of paper presented at the First Annual Meeting, 
International Omega Association, Arlington, Virginia, 

July 27-29, 1976 



35 MINI-0, SIMPLE OMEGA RECEIVER HARDWARE FOR USER EDUCATION 
R. W. Burhans, August 1976. 

Preprint of paper presented at the First Annual Meeting, 
International Omega Association, Arlington, Virginia, 

July 27-29, 1976. 

36 OPERATING INSTRUCTIONS: KENNEDY TEST FIXTURE 

Donald P. Seyler, September 1976. 

A unit for testing the integrity of data and control circuits 
of the Kennedy 1600/360 Incremental Tape Recorders is 
described. 

37 OPERATING INSTRUCTIONS: MEMODYNE/KENNEDY INTERFACE UNIT 

Donald P. Seyler, January 1977. 

A unit for transcribing incremental data from magnetic 
cassette tapes to magnetic reel-to-reel tapes is described. 

This includes integral testing for validity of data. 

38 KIM-1 INTERFACE ADAPTER TO 3-WIRE TELETYPE SYSTEMS 
R. W. Burhans, August 1976. 

This brief technical note has been submitted to the KIM-1 
microcomputer group publication, KIM User Notes. It is of 
interest to others who have 3-wire ASR-33 teletype systems in 
using microcomputer hardware with the Ohio University 
Prototype Omega Sensor Receivers. 

39 OMEGA DISTRIBUTION AMPLIFIER WITH FOUR CHANNEL INDEPENDENT 
LEVEL CONTROL 

Donald P. Seyler, September 1976. 

A portable unit for distributing a single small signal source 
to a maximum of four loads, with independent source level 
control for each load, is described. 

40 MEASURING CLOCK OFFSETS FOR MINI-0 WITH KIM-1 
R. W. Burhans, September 1976. 

The previous MINI— 0 Omega receiver system required adjustment 
of the local clock to a low offset for proper operation. 

Single station tracking loop software provides an easy way of 
determining the offset prior to experimental navigation tests. 

A second order software tracking loop is suggested to elimi- 
nate the local clock error problem. 

41-Mod 1 IMPROVED ANALOG OUTPUT CIRCUITS FOR OHIO UNIVERSITY PROTOTYPE 
OMEGA NAVIGATION RECEIVERS 
Lee Wright, October 1976. 

A minor hardware change is described which provides a more 
accurate analog output from the Ohio University Omega Prototype 
Receivers. 



DIURNAL MEASUREMENTS WITH PROTOTYPE CMOS OMEGA RECEIVERS 
R. W. Burhans, November 1976. 

The Ohio University Prototype CMOS Omega Sensor Processor is 
capable of receiving all eight Omega channels on 10.2 KHz. 
Diurnal recordings of selected station pairs made during the 
period October-November 1976 demonstrate the receiver perfor- 
mance and illustrate limitations for navigation using diurnal 
corrections. 

DEMONSTRATION PROGRAM FOR OMEGA RECEIVER PROTOTYPE MICRO- 
COMPUTER DATA PROCESSING 
R. W. Lilley, November 1976. 

Using the prototype Omega receiver developed for the NASA 
Joint University Program plus a digital interface to a 
commercial microcomputer, a software routine to demonstrate 
receiver operation is described. 

AN ASSEMBLER FOR THE MOS TECHNOLOGY 6502 MICROPROCESSOR AS 
IMPLEMENTED IN JOLT (TM) AND KIM-1 (TM) 

Robert W. Lilley, November 1976. 

The 6502 Assembler implemented at Ohio University for support 
of microprocessor program development in the Tri— University 
Program is described. 

INTERACTIVE OMEGA PROPAGATION CORRECTIONS 
Robert W. Lilley, January 1977. 

An implementation of Coast Guard computer programs for Omega 
propagation corrections is described. 

ANALYSIS OF A FIRST ORDER PHASE LOCKED LOOP IN THE PRESENCE OF 

GAUSSIAN NOISE 

Paul R. Blasche, March 1977. 

A first-order digital phase-locked loop Is analyzed by appli- 
cation of a Markov chain model. Steady-state loop error pro- 
abilities, phase standard deviation and mean loop transient 
times are determined for various input signal-to-noise ratios. 
In addition, results for direct loop simulation are presented 
for comparison. 

A MICROCOMPUTER-BASED LOW-COST OMEGA SENSOR PROCESSOR 
Richard J. Salter, Jr., February 1977, Master's Thesis. 

MINI-L LORAN-C RECEIVER 
R. W. Burhans, March 1977. 

A low-cost prototype Loran-C receiver front-end has been 
designed and bench-tested. This receiver concept provides 
outputs to interface with a microcomputer system. The devel- 
opment of sensor and navigation software for use with the Mini- 
L system is underway. 



SIMULATION OF DIGITAL PHASE-LOCKED LOOPS 
Paul R. Blasche, April 1977. 

This technical memorandum deals with development of simulation 
equations for first- and second-order digital phase-locked 
loops. In addition, examples of loop simulation are given to 
determine loop performance with respect to several loop para- 
meters. 

A KEYBOARD INTERFACE FOR THE JOLT MICROPROCESSOR 
Lee Wright, May 1977. 

The Ohio University Microprocessor Navigation Receiver Base 
utilizes the JOLT(TM) microcomputer. This keyboard interface 
is designed to allow data input without use of a teleprinter. 

A FOUR DIGIT MEMORY-MAPPED DISPLAY 
Ralph E. Smith, May 1977. 

An interface board has been fabricated for the Ohio University 
Microcomputer-Based Navigation Receiver to display data from 
the microprocessor. 

INTERACTIVE LORAN-C-TO-GEOGRAPHIC AND GEOGRAPHIC-TO-LORAN-C 
COMPUTATION 

Lynn M. Piecuch and Robert W. Lilley, August 1977. 

An implementation of Naval Oceanographic Office computer soft- 
ware for Loran-C is presented. 

CIRCUIT METHODS FOR VLF ANTENNA COUPLERS 
R. W. Burhans , September 1977. 

A summary of E-field antenna preamplifiers developed during 
the course of the NASA Tri-University Program studies on VLF 
methods for general aviation is presented. The circuit tech- 
niques provide useful alternative methods for Loran-Omega 
receiver system designers. 

LORAN-C DIGITAL WORD GENERATOR FOR USE WITH A KIM-1 MICRO- 
PROCESSOR SYSTEM 
James D. Nickum, December 1977. 

The digital word generator used with Mini-L front end to 
develop a Loran sensor processor at Ohio University is 
described. 

MICROPROCESSOR— TO-SYSTEM/3 70 INTERFACE 
Robert W. Lilley, February 1978. 

A hardware interface is described which allows direct memory 
load of a microprocessor from the host System/370 computer, 
eliminating paper tape handling. 


STAND-ALONE DEVELOPMENT SYSTEM USING A KIM-1 MICROCOMPUTER 
MODULE 

James Nickum, March 1978. 

Documentation of the stand-alone microprocessor development 
system used in the navigation sensor processor research at 
Ohio University is described. 

A LOW-COST LORAN-C ENVELOPE PROCESSOR (The Mini-L Loran-C 
Receiver) 

R. W. Burhans, April 1978. 

A reprint of published article on the Mini-L Loran-C receiver 
front-end is presented. Complete circuit details are given 
with a basic Introduction to Loran-C navigation and time- 
frequency standard uses. The Mini-L concept is of particular 
Interest to the low-budget experimenter as an RF signal Inter- 
face for more sophisticated end use. 

A VIDEO DISPLAY INTERFACE FOR THE LORAN-C NAVIGATION RECEIVER 
DEVELOPMENT SYSTEM 

Joseph P. Fischer and Robert W. Lilley, May 1978. 

A character-mode video unit is described which allows micro- 
processor-controlled display of program and navigation data 
with a small investment in logic. 

COMPUTING LORAN TIME DIFFERENCES WITH AN HP- 2 5 HAND 
CALCULATOR 

Edwin D. Jones, August 1978. 

Accurate Loran-C time differences can be calculated from 
known transmitter and receiver positions using the program 
described. 

PHASE-LOCKED TRACKING LOOPS FOR LORAN-C 
R. W. Burhans, August 1978. 

Two, portable, battery-operated Loran-C receivers have been 
fabricated to evaluate simple envelope detector methods with 
hybrid analog-digital phase-locked loop sensor processors. 

The receivers are being used to evaluate Loran-C in general 
aviation applications. Complete circuit details are given 
for the experimental sensor and readout system. 

LORAN-C FLIGHT TEST SOFTWARE 
James D. Nickum, August 1978. 

Described is the software package developed for the KIM-1 
Micro-System and the Mini-L PLL receiver to simplify taking 
flight test data at Ohio University. 



PREAMPLIFIER NOISE IN VLF RECEIVERS 
R. W. Burhans, September 1978. 

Rapid methods of estimating antenna preamplifier noise contri- 
bution to receiver performance are presented for JFET or CMOS 
transistors. An improved CMOS preamplifier circuit is 
suggested. 

LORAN-C TIME DIFFERENCE CALCULATIONS 
Joseph P. Fischer, October 1978. 

A simplified approach to calculate Loran-C time differences 
from a given geographic location is presented. 

INITIAL FLIGHT TEST OF A LORAN-C RECEIVER/DATA COLLECTION 
SYSTEM 

Joseph P. Fischer and James D. Nickum, November 1978. 

Described are the flight test results of a Loran-C navigation 
receiver/data collection system designed at Ohio University. 

ACTIVE ANTENNA FOR THE VLF TO RF OBSERVER 
R. W. Burhans, February 1979. 

This report is a prepublication manuscript submitted to one of 
the contemporary electronics magazines as part of a series on 
VLF-LF signal reception problems. The report presents a sim- 
ple and low-cost method of fabricating an active antenna pre- 
amplifier system covering the range of 10 KHz to 10 MHz, for 
use with tunable communications receivers. The same type of 
preamplifier system can be used with airborne VLF navigation 
receivers. 

ANALYSIS AND DESIGN OF A SECOND-ORDER DIGITAL PHASE-LOCKED 
LOOP 

Paul R. Blasche, March 1979. 

A second-order digital phase-locked loop is analyzed by 
application of a Markov chain model with alternatives. 
Steady-state loop error statistics and mean transient time 
are determined for various loop parameters. In addition, a 
hardware digital phase-locked loop was constructed and tested 
to demonstrate the applicability of the Markov chain mode. 

LORAN-C FLIGHT DATA BASE 
Robert W. Lilley, February 1979. 

A large file of Loran-C data to be used in receiver design 
and testing is documented. 

RESULTS OF THE SECOND FLIGHT TEST OF THE LORAN-C RECEIVER/ 

DATA COLLECTION SYSTEM 
Joseph P. Fischer, March 1979. 

Reported are the results of a second flight test of the 


Loran-C system under development at Ohio University using a 
variation of the techniques used for the first flight test. 

DIGITAL PHASE-LOCKED LOOP DEVELOPMENT AND APPLICATION TO 
LORAN-C 

Daryl L. McCall, September 1979. 

A digital phase-locked loop has been developed and implemented 
for use in a low-cost Loran-C receiver. This paper documents 
the DPLL design and application to Loran-C. 

ACTIVE ANTENNA COUPLER FOR VLF 
R. W. Burhans, November 1979. 

A reprint of a paper published in the "Ham Radio Magazine", 
Volume 12, Number 10, October 1979, is presented. The circuit 
designs are applicable to a variety of VLF-1IF active antenna 
receiving systems including Omega and Loran— C for airborne and 
marine users. 

EXPERIMENTAL LOOP ANTENNAS FOR 60 KHz to 200 KHz 
R. W. Burhans, December 1979. 

A series of loop antennas have been fabricated and evaluated 
for possible use with Loran-C and other VLF to LF band 
receivers. A companion low noise and very high gain pream- 
plifier circuit has been devised to operate the loop antennas 
remote from the receiver. Further work is suggested on the 
multiple loop antenna systems to provide omni-directional 
coverage and reduce E-field noise pickup in navigation or 
communications systems. 

DATA REDUCTION SOFTWARE FOR LORAN-C FLIGHT TEST EVALUATION 
Joseph P. Fischer, December 1979. 

This paper describes a set of programs written for use on 
Ohio University's 370 computer for reducing and analyzing 
flight test data. 

LORAN DIGITAL PHASE-LOCKED LOOP AND RF FRONT-END SYSTEM ERROR 
ANALYSIS 

Daryl L. McCall, December 1979. 

Various experiments have been performed to determine the 
system error of the DPLLs and RF front-end currently being 
used in a Loran receiver prototype. This paper documents 
those experiments and their results. 

RESULTS OF A LORAN-C FLIGHT TEST USING AN ABSOLUTE DATA 
REFERENCE 

Joseph P. Fischer, January 1980. 

The results of a flight test using the Loran-C receiver and 
data collection system developed at Ohio University are 
described in this paper. An absolute data reference was 



provided by the Lincoln Laboratories DABS (Discrete Address 
Beacon System) facility. 

ANALYSIS OF FIRST AND SECOND ORDER BINARY QUANTIZED DIGITAL 
PHASE-LOCKED LOOPS FOR IDEAL AND WHITE GAUSSIAN NOISE INPUTS 
Paul R. Blasche, March 1980. (Dissertation) 

Specific configurations of first and second order all digital 
phase-locked loops are analyzed for both ideal and additive 
white gaussian noise inputs. In addition, a design for a 
hardware digital phase-locked loop capable of either first 
or second order operation is presented along with appropriate 
experimental data obtained from testing of the hardware loop. 
All parameters chosen for the analysis and the design of the 
digital phase-locked loop are consistent with an application 
to an Omega navigation receiver although neither the analysis 
nor the design is limited to this application. 

DC-TO-DC POWER SUPPLY FOR LIGHT AIRCRAFT FLIGHT TESTING 
Stephen R. Yost, December 1980. 

A DC-to-DC power supply has been designed and fabricated to 
operate the prototype Loran-C receiver and data collection 
system currently in use at Ohio University. The supply is 
designed to operate from an aircraft electrical system. 

BI-DIRECTIONAL COMMUNICATION INTERFACE FOR MICROPROCESSOR-TO- 
SYSTEM/370 

Joseph P. Fischer, January 1981. 

Described is a hardware and software interface to allow two- 
way communication between a microprocessor system and the IBM 
System/370. 

LORAN-C PLOTTING PROGRAM FOR PLOTTING LINES OF POSITION ON 
STANDARD CHARTS 

James P. Roman, February 1981. 

The Loran-C plotting program was designed to plot Loran-C 
lines of position on any standard chart and is used in the 
data-collection system currently In use at Ohio University 
Avionics Engineering Center. 

AUTOMATIC GAIN CONTROL 
James P. Roman, March 1981. 

An automatic gain control has been designed and fabricated to 
operate with the Loran-C prototype receiver and data- 
collection system currently in use at Ohio University. 

A LORAN-C PROTOTYPE NAVIGATION RECEIVER FOR GENERAL AVIATION 
Robert W. Lilley and Daryl L. McCall, August 1981. 

The design, fabrication and evaluation of a prototype Loran-C 
receiver is described. Hardware is complete and microcom- 
puter programming continues for addition of area-navigation 



capability. The receiver is an envelope-processor, offering 
simplicity of RF processor circuitry. 

81 COMMUTATED AUTOMATIC GAIN CONTROL SYSTEM 
Stephen R. Yost, November 1981 

A commutated AGC system for the Ohio University prototype 
Loran-C receiver is described. The circuit design, fabrica- 
tion, and test results are presented in this paper. 

82 A PROTOTYPE INTERFACE UNIT FOR MICROPROCESSOR-BASED NAVIGATION 
SYSTEM 

Stanley M. Novacki III, November 1981 

A command entry and display device designed to allow con- 
venient operation of the Loran-C receiver-processor is 
described. 
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BI-DIRECTIONAL COMMUNICATION INTERFACE 
FOR MICROPROCESSOR-TO-SYSTEA4/37Q 


Described is a hardware and software interface 
to allow two-way communication between a 
microprocessor system and the IBM System/370 


Joseph P. Fischer 

Avionics Engineering Center 
Department of Electrical Engineering 
Ohio University 
Athens, Ohio 45701 




I. 


INTRODUCTION 


This paper documents the design and operation of a bi-directional 
communication interface between a microcomputer and the IBM System/370. The 
hardware unit inter-connects a modem to interface to the S/370, the microcomputer 
with an EIA I/O port, and a terminal for sending and receiving data from either the 
microcomputer or the S/370. Also described is the software necessary for the two-way 
interface. This interface has been designed so that no modifications need to be made 
to the terminal, modem, or microcomputer. This unit is designed to upgrade a 
uni-directional interface already in use [ 1 ] 

II. INTERFACE DESCRIPTION: HARDWARE 

Figure 1 shows the paths of signals between the microcomputer, the modem, 
and the terminal. The hardware interface consists of a four-pole, three-position 
switch and cables and plugs to connect the switch box to the other devices. All 
signals are assumed to be RS-232C (EIA standard). 

In switch position 1, the microcomputer is connected directly to the terminal; 
all communications are between these two only. The modem is isolated in this position 
and it is not necessary to have it connected if no communication to the S/370 is 
desired. In position 2, the serial out from the keyboard is routed to the modem for 
communicating to the S/370. The serial out from the modem goes to the terminal and 
the serial in of the microcomputer. In this position, it is possible to send commands 
and receive responses from the S/370, while the microcomputer reads the data sent by 
the S/370. Thus it is possible to load a program into the microcomputer by displaying 
the object file on terminal. It is necessary to switch to position 1 and issue the 
microcomputer load command prior to typing the file. Position 3 on the switch box 
connects the serial out from the modem to the terminal and to the serial in on the 
microcomputer. In addition, the serial out from the microcomputer is sent to the modem. 
Here, the microcomputer communicates directly with the S/370, the terminal always 
displays the response sent by the S/370. With proper positioning of the half-duplex/ 
full-duplex switches on the terminal and modem, the responses from the microcomputer 
may also be displayed. Note that the serial-out from the terminal is isolated, thus it 
may be necessary to start a program on the microcomputer by pressing the NMI (non- 
maskable interrupt) switch on the switch box. 

Table 1 lists the connection used on the terminal and modem. Connections for 
RS-232C are made through 25-pin D-connectors. Data terminal equipment (DTE) devices 
are supplied with a male (DB-25P) connector while data communication equipment (DCE) 
devices are supplied with a female (DB-25S) connector. Figure 2 shows the detailed 
routing of connections from the connectors on the terminal and modem through the 
switch box. 
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III. 


INTERFACE DESCRIPTION; SOFTWARE 


Full utilization of the bi-directional interface requires a set of programs to 
be run simultaneously on the microcomputer and the S/370. Figure 3 shows a block 
diagram of how the programs would operate for a typical application. Some points to 
be considered in writing the interface software are: 

a. Most microcomputers store character data internally as ASCII. 

b. Serial communications between devices are generally in ASCII 
format. 

c. The I/O routines for the S/370 expect to reoeive ASCII which is 
then converted to EBCDIC, which the S/370 uses for internal 
storage of character data . 

d. The Conversational Monitor System (CMS) portion of the VM/370 
operating system is line-oriented, i.e., no system action is taken 
until a carriage return (hex OD) is received . 

e. The S/370 issues a prompt when ready for another line. 


A typical application for which this interface has been used is transmitting 
data collected by the microcomputer on a cassette tape to the S/370, where it is 
stored on a disk file for further processing. The sequence of events is as follows: the 
data to be transmitted is stored in a buffer in the microcomputer's memory. Generally, 
80 characters comprise one line. Note that one byte consists of two four-bit hexa- 
decimal numbers, each of which is converted to ASCII. Thus if 80 characters are to 
be sent, the buffer is 40 bytes long. After 80 characters are sent, a carriage return 
(hex OD) is sent. The S/370 does the ASCIi-to-EBCDIC conversion and places the 
EBCDIC characters in a user buffer in the S/370 memory. When the S/370 is ready 
to receive another line, it sends a series of control characters. The microcomputer 
reads and recognizes these control characters as the prompt signal to send another 
line. The sequence of control characters currently sent by S/370 is shown in Figure 4. 

Appendix A gives a listing of a MOS Technology 6502 microcomputer program 
(intended to be run using the 'Super-Jolt' micro unit) for reading 40 bytes of data from 
a Memodyne digital cassette tape unit and sending these to the S/370. The data to 
be sent are packed BCD numbers; i.e., one BCD digit occupies four bits, two BCD 
numbers are contained in one byte. Each BCD digit is sent as ASCII by the 'output 
byte 1 routine in the microcomputer monitor program (at address 72B1 (hex) in the 
Super-Jolt (TM) monitor). A carriage return is sent at the end of the line with a call 
to the WRT routine at address 72C6 (hex). 
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Program lines 1 to 78 are initialization steps used for the Memodyne interface 
hardware and to position the tape properly. Lines 79 to 91 constitute the main part 
of the program which builds up the 40-byte buffer then sends the buffer to the S/370. 
This part loops continuously whether or not any data is received. The operator 
should monitor the operation to stop the program when all the data has been transmitted. 
Subroutine READ is called to read a byte from the tape unit. Subroutine W370 sends 
the 40-character buffer to the S/370, sends a carriage return, then looks for a period 
(hex 2E) followed by a DC-1 (hex 11). If this sequence is not done, the S/370 
issues a read-error message. When these two characters are received, control is 
passed back to the main program sequence. 

Subroutine RDT is a modification of the RDT routine at address 72E9 in the 
Jolt monitor. Most serial-read routines on microcomputers are full-duplex; as each 
bit is received, it is echoed back out to the sending device. However, the S/370 
can receive half-duplex only. Thus it is necessary to change the interface method 
through the modem or to re-write the read routine so that the received bits are not 
echoed by the microcomputer. This is the purpose of having a separate read routine. 

If this is not done, read-errors result. The program presented here is shown to 
illustrate one application of the bi-directional interface. Other uses on other micro- 
computers would still use the same basic philosophy. 

The companion program that is run on the S/370 is shown in Appendix B. 

This program is written in IBM 360/370 assembler language [ 2 ] using standard 
CMS I/O routines. Again this program illustrates the application of sending data to 
the S/370 for storage on a disk file. 

The data is read 80 bytes at a time, each BCD character in its ASCII format. 
Each character read is stripped of the upper four-digit mask and is repacked. This is 
done by the translate instruction at line 88 and the PACK instruction at line 90. 

Since one record produces only 40 packed BCD digits, two lines are read before one 
80-byte record is written to the file. A blank line or an incomplete line is filled to 
the end with zeros. Each time a record is written, a counter is incremented which is 
printed at the end of program execution. 


IV. INTERFACE OPERATION 

The example of transmitting data from the microcomputer to the S/370 will 
be continued here to show how the interface may be operated. After the interface 
is properly connected, power should be applied to all units. At this point it is usually 
necessary to load the microcomputer with a program stored on a disk file. Thus the 
switch box should be set to position 2 and the appropriate CMS LOGON procedure 
performed. When the microprocessor object code is ready for transmittal (through 
editing, assembling, simulating, etc.) the switch box should be set to position 1, 
the microcomputer reset button pushed, and a carriage return or other appropriate 
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key to reset the microcomputer typed. Then issue the proper command to set the 
microcomputer for loading hexadecimal data over its serial lines. The switch box is 
then set back to position 2 and the appropriate command is issued to the S/370 to 
load the microcomputer with the object file. Next, the unit is switched back to 
position 1 to verify correct loading, initialize any memory locations and set up the 
NMI vector address to the start of the program. Now the switch box is placed in 
position 2 and the program to receive the data is started and then the unit is set to 
position 3 and the NMI button pressed. 

As operation commences, the prompting period and any other responses 
from the S/370 will be displayed on the terminal. Depending on the setting of the 
half-duplex/full-duplex switches on the terminal and modem, data sent by the 
microcomputer will also be displayed on the terminal. 

When the operation is finished, the unit may be set to position 2 to stop 
the S/370 program then position 1 to stop the microcomputer program. 


V. SUMMARY 

A discussion was presented here of an interface unit and software procedures 
to allow two-way communication between a microcomputer and a central computer. 
This can be used for two-way data transmission, control and other applications where 
bi-directional communications are necessary. As an aid to setting up the software 
for other computer systems, ASCII [ 3] and EBCDIC [4] tables are given in Tables 2 
and 3. 
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DTE 


DCE 


P GND 
S OUT 
S IN 
RTS 
CTS 
DSR 
S GND 

RLSD 

DTR 


1 


1 


2 

2 


3 

3 


4 

4 


5 

5 


6 

6 


7 

7 


8 

8 


20 

20 



Protective Ground 
Serial Output 
Serial Input 

Request to Send 
Clear To Send 
Data Set Ready 
Signal Ground 

Received Line Signal Detect 
Data Terminal Ready 


Table 1 . RS-232C Connections 






000 

001 

010 

Oil 

100 

101 

110 

in 

0000 

NULL 

©DC 0 

1) 

0 

<e> 

P 



0001 

SOM 

DC, 

1 

1 

A 

Q 



0010 

EOA 

DC, 

•• 

2 

B 

R 



0011 

EOM 

dc 3 

# 

3 

C 

S 



0100 

EOT 

dc< 

Utop) 

$ 

4 

D 

T 



0101 

WRU 

ERR 

% 

5 

E 

U 



0110 

RU 

SYNC 

& 

6 J 

F 

V 



0111 

BELL 

LEM 

• 

7 

G 

w 



1000 

o 

UJ 

LL 

So 

( 

8 ■ 

H 

X 


1001 


s, 

) 

. 9 

1 

Y 



1010 

LF 

s 2 

• 

: 

J 

z 



1011 

Vtab 

s 3 

+ 

,* 

K 

! 



1100 

FF 


Icommil 

< 

L 

\ 


ACK 

1101 

CR 

s 5 

■' ; - / 

m 

M 

) 


(?) 

1110 

SO 

s 6 


> 

N 

r 


ESC 

1111 

St 

$7 

' ' : 1 / ■ 

? 

0 

*— 


DEL 


Example; 


0001 


-b, 


100 


b 7 - 


- A 


The abbreviations used in the figure mean: 


NULL 

Null Idle 

CR 

Carriage return 

SOM 

Start of message 

SO 

Shift out 

EOA 

End of address 

SI 

Shift in 

EOM 

End of mossage 

DC 0 

Dcvico control Q) 
Reserved for data 
Link escape 

EOT 

End of transmission 

DC, - DC 3 

Device control 

WRU 

"Who arc you?" 

ERR 

Error 

RU 

"Aro you ....?" 

SYNC 

Synchronous idle 

BELL 

Audible signal 

LEM 

Logical end of media 

FE 

Format effector 

SO 0 - SO; 

Separator (information) 

HT 

Horizontal tabulation 


Word separator (blank, 
normally non-printing) 

SK 

Skip (punched card) 

ACK 

Acknowledge 

LF 

Line feed 

© 

Unassigned control 

V/TA8 

Vertical tabulation 

ESC 

Escape 

FF 

Form feed 

DEL 

Delete Idle 


Table 2 


ASCII Table 






- 1 

* I 

t ^ i 

.5 i S' 

i J * 

it % AO 

i “ 

j #i _ . 

i ■• 

L . " 

m 

*Mill«H 0, 1 

Aliens 2,3 
Heaodectmel Digit 

LfL 

id 

ID 

ID 

00 

01 

10 

II 

00 

01 

ID 

ID 

ia 

ID 

ID 

D 

i *" 

B 

id 

2 

3 

4 

5 

4 

7 

• 

f 

m 

1 

ID 

ID 

ID 

f 

M 


Ifcih 




m 


2 one Punches 

Digit Punches 


E 

i 

B 

Ei 

m 

m 

m 

m 

fl 

IE 

B 

m 

B 

l 

fl 

IE 

Kjj 

IB 



g 

1 

yv:(f 

.'Vi 

TJO'A 

S 
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SOS 




fl 

B 

B 

B 


B 

fl 

B 

| 

1 

Mi 

.B^D 

B 

E 

s 

atu~. 

B 
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fS 

SYN 





m 

B 

• 


B 

B 

m 

IB 

m 



D 

H 

B 

TM 







c 

B 

■a 


a 

B 

■ 

3 

3$ 

ttitX 



D 

u 

9 

ms 

•YP 

PN 




■ 

g 

1 

B 

1 

a 

M 

u 

4 

fl 


BP 

B 

i 

■ 

Nt 

Lf 

IS 




g 

g 

n 

B 

■ 

a 

N 

V 

3 




a 

M 

c-iiia 

1C 

IS 

in 

uc 




■ 

E 

o 

1 

E 

5 

O 

w 

4 

a 


EB 

0 

u 

kk--w 

Ml 

II 

CSC 

COT 





9 

P 

B 

u 

B 

B 

B 

7 

m 


G 

D 


GI 

CAN 



~ 


B 

B 

B 

n 




B 

B 

• 


H 

D 

iS 

a 

Cm 






3 

1 

r 

1 

1 

1 

a 


Bl 

SSs 

B 

B 

m 

usm 

9 

CC 

SM 


D 

fl 

m 


U 


1 

a 




Bl 

fl 

gjj 

i 

PI'7 

® 

a 

CUl 

CU2 

CU3 


s 

m 

3 

B 


B 

B 

■ 


i • • 




0 

shti 


If* 


OQ4 

< 

• 

B 

Q 

B 




a 


■ .< 


i 

g^i 

B 

1 


IGS 

INO 

NAK 

■ 

fl 

■: 

g 
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D 

i 


a 

ACK 


a 

B 

Bl 
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B 

B 
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fli 


EDI 

D 

i 

SI 
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*u 
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? 
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Cord HaU Pottery 


o 

12-0-9-4-1 

© 

No Punch** 

® 12-0 ® 0-1 




0 

12-H -9-1-1 


12 

®n-o ® ii-o-»-i 



0 
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0 

II 
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© 
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Control Character Representation* 





SpeClol Graphic Character* 



ACK 

Acknowledge 

cor 

End ol Transmission 

PF 

Punch Olf 

« 

Cent Sign 

> 


KL 

Soil 

CSC 

Ctcope 

PN 

Punch On 


Period, Decimal Point 

? 


IS 

ftockipoce 

ETS 

End of Transmission Hock 

RES 

Redone 

< 

Lou-lhon Sign 


Grove Accent 

BYP 

Sypou 

ETX 

End of Teat 

RLF 

RevetM line Feed 

( 

lei* Porenthe*!* 



CAN 

Cancel 

FF 

Form feed 

RS 

Reoder Slop 

« 

f lie Sign 

§ 

Number Sign 

CC 

Curvx Control 

FS 

Field Seporalor 

SI 

Shift In 

1 

logical OR 

r< \ 

At Sign 

CR 

Carriage Return 

GE 

Graphic f nope 

SM 

Set Mode 

& 

Ampenond 

• 

Prime, Apostrophe 


Cuitomer Um 1 

HT 

Horljootol Tob 

SMM 

Sferl of Morwol MeivoQe 

1 

Evclamolion Point 

■ 

Equal Sign 

cm 

Cudomer Um 2 

IFS 

Interchonge File Sepo rotor 

so 

Shift Out 

s 

Dollar Sign 

■ 

Quotation Mork 

CU3 

Customer Um 3 

IGS 

Interchange Group Separator 

SOH 

Stort of Heeding 

• 

AtterUk 


Tilde 

DCI 

Device Control 1 

IL 

Idle 

SOS 

Start of Significance 

) 

Right Potentheii* 

( 

Opening Iroco 

DC2 

Device Control 2 

IRS 

interchonge Record Seporator 

SP 

Space 


Semicolon 

s 

Hook 

DC4 

Device Control 4 

IUS 

Interchonge Unit Se pore tor 

Six 

Stort of Teat 

~l 

Logical NOT 

Y 

Fork 

DEL 

Delete 

CC 

lower Com 

SUI 

Subtlihito 

- 

Mlnw* Sign, Hyphen 

i 

Closing fergee 

oie 

De»* Link E scope 

LF 

line feed 

SYN 

Synchronout Idle 

/ 

Sloeh 

\ 

Reverie Sion* 

os 

Digit Select 

NAK 

Negellve Acknowledge 

TM 

Tope Morb 



rl 

Choir 

CM 

End ef Medium 

NL 

New line 

UC 

Upper Com 



I 

long Vert! col Mork 

CNQ 

Enquiry 

NUl 

Nu» 

VT 

Vertlcel Tob 

% 




EO 

Eight Onet 






Underscore 




Table 3. EBCDIC Table. 
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Figure 3. Control Program Flow Charts. 






























APPENDIX A 


Program Listing for Microcomputer Control Program 


FILE: ONLOAD S6502 A OHIO UNIVERSITY AVIONICS ENGINEERING CENTEB 


*********************** *********************************************** 

* 

THIS PROGRAM IS DESIGNED FOR RUNNING ON THE JOLT/MEHODYNE 
SYSTEM FOR RECOVERING DATA STORED ON THE DIGITAL TAPE. 

THE DATA IS READ IN 40 BYTES AT A TIME AND STORED IN A 
BUFFER. THEN THE BUFFER IS SINT TO THE. S/370 OVER THE JOLT*S 
SERIAL LINES. WITH ASCII CONVERSION, 80 BYTES ARE ACTUALLY 
SENT OVER THE SERIAL LINES. 


J.P. FISCHER 


08/1980 


IJNL00010 

UNL00020 

UNL00030 

UNL00040 

UNL00050 

UNL00060 

UNL00070 

UNL00080 

UNL00090 

UNL00100 

UNL00110 


********************************************************************** UNL00120 


UNL00130 

UNL00140 


PIAA 

EQU 

$4000 

ADDRESS OF PIA SIDE A 

UNL00150 

PI AB 

EQU 

$4002 

ADDRESS OF PIA SIDE B 

UNL00160 

BRT 

EQU 

S72C6 

JOLT WRITE DATA TO SERIAL OUT LINE 

UNL00170 

WROB 

EQU 

$72B1 


UNL00180 

MPB 

EQU 

56E02 

PIA B FOR SERIAL I/O WORK 

UNL00190 

HCLK1T 

EQU 

S6E04 

PIA TIMER 

UNL00200 

BCLKRD 

EQU 

$6E04 

SAME AS ABOVE 

UNL00210 

MCLK IF 

EQU 

$6E05 

SOME MORE TIMER STUFF 

UNL0Q220 

MAJCRT 

EQU 

SEA 

UPPER 8 PITS OF EAUD RATE 

UNL00230 

MINCRT 

EQU 

SEB 

LOWER 8 BITS OF BAUD RATE 

UNL00240 

TAPESY 

EQU 

*000000 1 0 PATTERN FOB TAPE SYNC CHECK 

UNL00250 

BOT 

EQU 

X00000100 PATTERN FOR BOT/EOT CHECK 

UNL00260 

LF 

EQU 

*00010000 PATTERN FOR LOAD FOWARD FUNCTION 

UNL00270 

REN2 

EQU 

*00100000 PATTERN FOR BEHIND OPERATION 

UNL00280 

START 

EQU 

*10000000 PATTERN FOB INITIATING START 

UNL00290 

* 




UNL00300 


OBG 

0 


UNL00310 

XTEMP 

BSS 

1 

TEMPORARY FOR X 

UNL00320 

YTEMP 

BSS 

1 

TEMPORARY FOR Y 

UNL00330 

PUFFER 

BSS 

40 


UNL00340 

* 




UNL00350 

* 




UNL00360 


ORG 

S200 


UNL00370 


J SR 

I NIT 

SET UP PIA FOR MEMODYNE 

UNL00380 


LDA 

PI AB 

PREPARE TO CHECK BOT 

UNL00390 


AND 

= BOT 

SEE IF ON LEADER 

UNL00400 


BNE 

NOTBOT 

IF NOT, THEN OK 

UNL00410 


LDA 

PIAB 

GET SIDE B 

UNL00420 


EOR 

= LF 

CLEAR LOAD FOWARD BIT 

UNL00430 


ORA 

* ST ART 

SET START BIT HIGH 

UNL00440 


STA 

PIAB 

AND STORE TO LOAD FOWARD 

UNL00450 

BTLOOP 

LDA 

PIAB 

GET STATUS 

UNL00460 


AND 

3 BOT 

SEE IF STILL ON LEADER 

UNL00470 


BEQ 

BTLOOP 

CONTINUE TESTING UNTIL OFF 

UNL00480 


LDA 

PIAB 


UNL00490 


ORA 

=LF 

SET LOAD FOWARD HIGH 

UNL00500 


STA 

PIAB 

AND REPLACE 

UNL00510 

IDLOOP 

LDA 

PIAB 


UNL00520 


AND 

*BOT 

NOW LOOP UNTIL AT READY POINT 

UNL00530 


BNE 

LDLOOP 

KEEP GOING UNTIL ON HOLE 

UNL00540 


LDA 

PIAB 


UNL00550 
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FILE: ONLOAD S6502 


A 


OHIO UNIVERSITY AVIONICS ENGINEERING CENTEE 



EOR =LF 

SET LOAD FOHARD I0W TO HOVE 

UNL00560 


STA PI AB 

OFF OF HOLE 

UNL00570 

TLOOP 

LDA PIAB 


UNL00580 


AND =BOT 


UNL00590 


BEQ TLOOP 


UNL00600 


LDX =$80 

TIBER ROOTINE 

UNL00610 

TIHXT 

DEX 


UNL00620 


BNE TIH XT 

KEEP LOOPING UNTIL OUT 

UNL00630 


ORA *LF 

NOW RETURN LOAD FOWARD 

UNL00640 


STA PIAB 

HIGH, SHOULD BE OFF OF HOLE 

UNL00650 

NOTBOT 

LDA =0 

CLEAR ACCUM. AND SET 

UNL00660 


STA PIAA*1 


UNL00670 


STA PIAB*1 


UNL00680 


LDA =$F8 


UNL00690 


STA PIAB 


UNL00700 


LDA =0 


UNL007 1 0 


STA PIA A 


UNL00720 


LDA =$FF 


UNL00730 


STA PIAB*- 1 


UNL00740 


STA PI A A* 1 


UNL00750 


LDA PIAB 


UNL00760 


AND =$11110111 

UNL00770 


STA PIAB 


UNL00780 

* 



UNL00790 

* 

NOW INITIALIZE THE 370 AND START 

UNL00800 

* 

SENDING 

DATA. 

UNL008 1 0 

* 



UNL00820 

LFLDS 

LDX =0 

READ 40 CHARACTERS FROM TAPE 

UNL00830 

L80 

JSR READ 

GET A BYTE FROH RECORDER 

UNL00840 


STA BUFFER, 

X SAVE IN OUTPUT BUFFER 

UNL00850 


INX 

DO ANOTHER ONE 

UNL00860 


CPX =40 

DONE 40 BYTES YET? 

UNL00870 


B HE LOO 

IF NOT, DO AGAIN 

U NL00880 


LDY =40 

SEND THESE 40 

UNL00890 


JSR W37 0 

SEND TO SYSTEH 

UNL00900 


JHP LFLDS 


UNL00910 

* 



UNL00920 

* 

INITIALIZATION FOR PIA 

UNL00930 

* 



UNL00940 

INIT 

LDX =0 


U NL00950 


STX PI A A* 1 


UNL00960 


STX PIAB* 1 


UNL00970 


STX PI A A 


UNL00980 


LDA =$B8 


U NL00990 


STA PIAB 


UNL01000 


LDA =$FF 


UNL010 10 


STA PI A A* 1 


UNL01020 


STA PIAB* 1 


UNL01030 


LDA =0 


UNL01040 


ORA =LF 


UNLO 1050 


ORA =REW2 


UNL01 060 


STA PIAB 


UNL01070 


RTS 


UNL01080 

* 



UNLO 1090 

* 



UNL01 100 
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FILE: ONLOAD S6502 A OHIO UNIVERSITY AVIONICS ENGINEERING CENTEI 


*********************************************************************** UNLO 1110 

* * UNL01120 

* THIS IS THE READING PORTION CP THE PROGRAM TO RECOVER * DNL01130 

* DATA FROM THE RECORDER AND PLACE IN THE MICROCOMPUTER * S * UNL01140 

* MEMORY. * UNLO 1150 

* * UNL01160 
*********************************************************************** OHLO 1 170 

* UNLO 1180 

* UNL01190 


READ 

LDA 

PI AB 


U NL0 1 200 


ORA 

=START 

SET START HIGH 

UNLO 1210 


STA 

PIAB 


UNLO 1220 

RDLP 

LDA 

PI AB 


UNL01230 


AND 

=TAPESY 

WAIT UNTIL SYNC IS HIGH 

UNL01240 


BEQ 

RDLP 


UNLO 1250 


LDA 

PIAB 


UNL01280 


EOR 

=START 

SET START LOW AGAIN 

UNL01270 


STA 

PIAB 


U NLO 1280 

INLP1 

LDA 

PIAB 


UNLO 1 290 


AND 

-TAPESY 

WAIT UNTIL SYNC IS LOW 

UNLO 1300 


BNE 

INLP1 


U NLO 13 1 0 


LDA 

PIA A 

GET THE DATA FROM RECORDER 

UNLO 1 320 


RTS 



UNLO 1 330 


* ****** **************** *********** ************************************ 

* 
* 
* 
* 
* 
* 
* 
* 
* 
* 


THIS SUBROUTINE 0 
THE ADDRESS OF TH 
INDEXED BY THE X- 
TO BE SENT IS CON 
BUFFER IS SENT, A 
FOR THE CONTROL C 
TO BE SENT BACK, 
CHARACTER AFTER T 
IS IN THE READ ST 


UTPOTS A LI 
E BUFFER IS 
REGESTER. 
TAINED IN T 
•CR* IS SE 
HAHACTERS B 
THEN WAITS 
HE PERIOD I 
ATE. 


ME OF CHARACTERS TO THE S/370. 

IN PAGE ZERO AND IS 
THE LENGTH OF THE BUFFER 
HE Y-REGESTER. AFTER THE 
NT THEN THE PROGRAM WAITS 
ETWEEN THE *CR* AND PERIOD 
FOR THE CONTROL 
NDICATING THE S/370 


*********************** *********** ************************************ 


UNL01340 
U NLO 1350 
UNLO 1 360 
UNLO 1370 
UNL01380 
U NLO 1390 
UNLO 1 400 
UNLO 1410 
U NLO 1420 
UNL01430 
UNLO 1440 
TJNL01450 
U NLO 1460 
UNL01470 
UNL01480 
U NLO 1490 


W370 

LDX 

= 0 

POINT TO FIRST CHARACTER 

UNL01500 


STX 

XTEMP 

ZERO X-TEHP SPACE 

UNLO 1510 


STY 

YTEMP 

SAVE LENGTH 

UNL01520 

WLOOP 

LDX 

XTEMP 

GET POINTER 

U NLO 1 530 


LDA 

BUFFER, 

X GET A CHARACTER 

UNLO 1540 


JSR 

WROB 

SEND IT 

UNL01550 


INC 

XTEMP 

X*1 

UNL01560 


DEC 

YTEMP 

LESS ONE CHARACTER 

UNL01570 


BNE 

WLOOP 

GO AG AIN IF NOT CONE 

UNLO 1 580 


LDA 

=.*D 

CARRIAGE RETURN 

UNLO 1590 


JSR 

WRT 

TELL 370 THIS IS END- OF- LI NE 

U NLO 1600 

SCANP 

JSR 

RDT 

READ JUNK FROM SISTEH 

UNLO 1610 


CMP 

= $2E 

PERIOD 

UNLO 1620 


BNE 

SCANP 


UNL01630 


JSR 

RDT 

LOOK FOR 

UNL01640 


CMP 

«J1 1 

DC1 

UNLO 1650 
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FILE: UNLOAD S6502 


OHIO UNIVERSITY AVIONICS ENGINEERING CENTEI 


BNE SCANP 

RTS RETURN TO CALL ER 


*********************** *********** ************************************ 


UNL01660 
UNL01670 
UNL01680 
UNLO 1690 
UNL01700 


* UNL01710 

HIGH SPEED REWIND. * UNL01720 

* UNL01730 

***************************************************************** UNLO 1740 

J SR 

I NIT 

UNL01 750 
UNL01760 

LDA 

=$B8 

UNLO 1770 

STA 

PIAB 

UNLO 1780 

LDA 

=LF 

UNL01790 

STA 

PIAB 

UNLO 1800 

ORA 

= REW2 

UNL01810 

STA 

PIAB 

UNL01820 

BRK 


UNL01830 
UNLO 1840 
UNLO 1 850 


********************************************************************** UNLO 1860 

* UNL01870 

MODIFIED JOLT READ ROUTINE. * UNL01880 

THIS ROUTINE IS IDENTICAL TO THE ORIGINAL 'ROT* ROUTINE * UNL01890 

AT ADDRESS S72E9, BUT THIS ROUTINE OPERATES IN HALF- * UNL01900 

DUPLEX RATHER THAN FULL- DO PL EX MODE. * UNL01910 

* UNLO 1920 

********************************** ************************************ UNLO 1930 

UNL01940 


PDT 

LDX 

= 8 


UNL01950 

* 




UNL01960 

RDT1 

LDA 

MPB 

WAIT FOR START BIT 

UNLO 1970 


LSR 

A 


UNL01980 


BCC 

RDT1 


UNL01990 

* 




UNL02000 


JSR 

DLY 1 


UNL020 10 

* 




UNL02020 

FDT2 

JSR 

DLY2 


UNL02030 


LDA 

MPB 

CY = NEXT BIT 

UNL02040 


LSR 

A 


UNL02050 

* 




UNL02060 


PHP 


SAVE BIT 

UNL02070 


TYA 


Y CONTAINS CHAR EEING FORMED 

UNL02080 


LSR 

A 


UNL02090 


PLP 


RECALL BIT 

UNL02100 


BCC 

BDT4 


UNL021 10 


ORA 

-J80 

ADD IN NEXT BIT 

UNL02120 

RDT4 

TAY 



UNL02130 


DEX 



UNL02140 


BNE 

RDT2 

LOOP POE a BITS 

UNL02150 


EOR 

= $FF 

COMPLEMENT DATA 

UNL02160 


AND 

= $7F 

CLEAR PARITY 

UNL02170 


JSR 

DLY 2 


UNL02180 


CLC 



UNL02190 


UNL02200 
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FILE: UNLOAD S6502 


A 


OHIO UNIVERSITY AVIONICS ENGINEERING CENTEi 


DLY2 

JSR 

DLY1 

UNL02210 

* 



UNL02220 

DLY1 

PHA 

SAVE FLAGS AND A 

UNL02230 


PHP 


U HL0 2240 


TXA 


UNL02250 


PHA 

SAVE X 

UNL02260 


LDX 

flAJCRT 

UNL02270 


LDA 

HINCRT 

U NL02280 

* 



UNL02290 


STA 

HCLK1T 

UNL02300 

DL3 

LDA 

HCLK1F 

UNL02310 


BPL 

DL3 

0NL02320 


DEX 


UNL02330 


PHP 


UNL02340 


LDA 

HCLKRD RESET TIMER INT FLAG 

UNL02350 


PLP 


UNL02360 


BPL 

DL3 

DNL02370 

* 



UNL02380 


PLA 


UNL02390 


TAX 


UNL02400 


PLP 


UNL02410 


PLA 


UNL02420 


RTS 


UNL02430 

* 



UNL02440 

* 



UNL02450 


ORG 

SPFFA 

UNL02460 


HEX 

00,02 

UNL02470 

* 



U NL0 2480 


END 


UNL02490 
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APPENDIX 


B 


Program Listing for S/370 Control Program 


FILE: UNLOADSS ASSEMBLE A OHIO UNIVERSITY AVIONICS ENGINEERING CENTS 



TITLE 

•UNLOADSS: 

READS RECORDS FROM MEMODYNE/MICROCOMPUTER 

IN*U NLOO 0 1i) 



TERFACE AND 

STORES ON DISK.' 


UNLOO 020 


PRINT 

NOGEN 



UNLOO 030 


SPACE 




U NLOO 040 

******************** ******* ****************** ************** **** 41*41 ***** UNL00 050 

* 





* UNL00060 

* 

THIS 

PROGRAM IS DESIGNED TO BE RUN ON THE S/370 

IN CON- 

* • U NLOO 070 

* 

JUCTION WITH THE MICRO 'UNLOAD* PROGRAM AND THE 

MEMODYNE/ 

* U NLOO 080 

♦ 

MICROCOMPUTER HARDWARE INTERFACE. RECORDS READ 

FROM 

* U NLOO 090 

* 

TAPE 

BE THE MICRO 

ARE SENT TO THE 370 IN ASCII, 

CP THEN 

* H NLOO 100 

* 

TRANSLATES THESE TO EBCDIC WHICH MUST BE TRANSLATED 

* U NLOO 110 

* 

BACK 

TO HE* BY THIS PROGRAM. 80 BYTES ARE SENT 

AT A 

* UNL00 120 

♦ 

TIME 

(40 EQUIVALENT HEX CHARACTERS) AND 80 HEX 1 

CHARACTERS 

* UNLOO 130 

* 

ARE STORED ON THE 

DISK FILE. 


* nNL 00 140 

♦ 





* UNLOO 150 

* 


J. P. FISCHER 0e/1980 


* U NLOO 1 60 

* 





* UNLOO 170 

******************** ************* *** ************ ****** ***************** [j NX. 00 1 80 


SPACE 

2 



U NLOO 190 

UNLOADSS 

START 

X* E000 • 



U NL00200 


USING 

UNLOADSS, 12 



U NLOO 2 10 


MV I 

FLAGS, 0 

CLEAR ALL FLAG BITS 


UNLOO 2 20 


LA 

1 » 8 (f D 

POINT TO FILE NAME FIELD 


UNLOO 230 


LR 

2,1 

SAVE PLIST ADDRESS 


UNLOO 240 


CLI 

0 (1) , X ' FF* 

BLANK ? 


UNLOO 250 


BE 

NOID 

IF SC, ERROR 


UNLOO 260 


LA 

1 . 8 (, 1 ) 



UNLOO 270 


CL T 

0 (1) , X'FF' 

NO FILETYPE? 


UNLOO 280 


BE 

NOID 

IF NOT, ERROR 


UNL00290 


MV C 

FILEID + 8 (16) 

,0(2) MOVE PARTIAL ID 


UNLOO 300 


LA 




U NLOO 3 10 


CLI 

0(1) , X * FF ' 

NO FILEMODE 


U NLOO 320 


BE 

NOMODE 

IF NOT SUBSTITUTE 'A' 


U NLOO 330 


MVC 

FILSID+24 (2) 

,16(2) MOVE IN NEW MODE 


UNLOO 340 


B 

CHECK 

CONTINUE 


UNLOO 350 

NOMODE 

MVI 

FIL El D * 24 , C * 

A * MOVE IN 'A* 


UNL09360 


MV I 

FILEID*25,C' 

• 


UNLOO 370 


SPACE 




UNLOO 380 

CHECK 

LA 

1,0 (, 1 ) 

MOVE POINTER UP SOME MORE 


UNL 00 390 


CLI 

0 (1) , X'FF* 

SEE IF ANYTHING THERE 


UNLOO 4 00 


BE 

CHECK 1 

IF NOT, CONTINUE 


U NLOO 4 10 


CLI 

0(1) ,C' {• 

SEE IF OPTION 


UNLOO 420 


BN E 

PARMERR 

IF NOT, BAD PARM 


II NLOO 430 


LA 

1 , 8 (, 1 ) 

NEXT FIELD 


UNLOO 440 


CLI 

O 

X 

SEE IF BLANK 


UNLOO 450 


BE 

CHECK 1 



UNLOO 460 


CLC 

0 ( 8 , 1) ,OPTREP SEE IF REPLACE OPTION 


UNLOO 470 


BNE 

BADOPT 

IF NOT, CONTINUE 


UNLOO 480 


01 FLAGS, 1 

SET REPLACE BIT 


II NLOO 490 


SPACE 




UNLOO 500 

CHFCK1 

TM 

FLAGS , 1 

SEE IF REPLACE IN EPFECT 


U NLOO 5 10 


BZ 

OPENF 

IF NOT, GO ON 


UNLOO 520 


FSERASE FSCD=FILEID 


U NLOO 530 

OPENF 

FS OPEN 

1 FSC B = FILEID 

OPEN FOR WRITING 


UNLOO 5 40 


CL 

1 5, F3 6 

SEE IF INVALID DISK 


UNLOO 550 
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FILE: UNLOADSS ASSEMBLE A 


OHIO UNIVERSITY AVIONICS ENGINEERING CENTS 


BE INVDISK 
SPACE 2 

**** *** **** ************************* ********** ****** ****** ******* ****** 

* 


THIS PART OF THE PROGRAM CAUSES A TERMINAL 

READ TO GET THE ASCII CHARACTERS, THEN TRANSLATES 

THEM TO HEX AND STORES ON DISH. 


* 
* 
* 
* 

********************************************************************* * 

SPACE 


RDLOOP 


STRIPE 


Z1 


DONE 1 


SLR 
SLR 
LA 
LA 
L 

LA 

RDTERM 

LTR 

BZ 

WAITT 

TR 

MVC 

PACK 

MVC 

LA 

BXLE 

SPACE 

RDTERM 

LTR 

BZ 

WAITT 

LA 

TR 

MVC 

PACK 

MVC 

LA 

BXLE 

LA 


4,4 

7.7 

2 , W BU F 

4.8 

5 , IBU F80 
3 , IBU F 
I BU F 
0,0 
DONE 

IBUF ( 80) ,TRTBL 
TEMP (8) ,0 (3) 

TEMPI (5) ,TEMP(9) 

0 (4,2) .TEMPI POT 
2,4(2) 

3,4, STR IPZ 


CLEAR RECORD COUNTER 

WRITE BOFFER ADDRESS 

LOOP INCREMENT 

END OF LOOP 

READ BUFFER ADDRESS 

GET A RECORD 

SEE IF NULL LINE 

GO IF IT IS 

WAIT FOR I/O 

CHANGE TO HEX 

GET 8 ZONED BYTES 

REMOVE THE ZONES 
PACKED CHARS. IN OUT BUFFER 
NEXT POSITION IN OUTPUT BUFFER 
CONTINUE UNTIL WHOLE RECORD DONE 


IBUF 
0,0 
DONE 1 


3, IBUF 

I BUF ( 80 ) ,TRTBL 
TEMP (8) ,0 (3) 
TEMPI (5) .TEMP (9) 
0(4,2) .TEMPI 
2,4 (2) 

3,4,Z1 
2.WBUF 


GET ANOTHER 
SEE IF NULL 
GO IF IT IS 
WAIT FOR I/O 
RE-INITIALIZE 


80 CHARS. 
LINE 


POINTER 


FS WRITE FSCB=FILEID 


LTR 

BNZ 

LA 

B 

SPACE 
M VI 
MVC 


15,15 
WRTERR 
7,1 (7) 
RDLOOP 


GET 8 BYTES 

REMOVE ZONES 
PUT IN OUT BUFFER 
NEXT LOCATION 
DO 80 BYTES 
REINTI ALIZ E WRITE 
SEND TO DISK 
SEE IF ERROR 
GO IF THERE IS 
ADD ONE TO RECORD 
PROCESS SOME MORE 


POINTER 


COUNT 


WBUF* 40 , 0 PREPARE TO CLEAR 

WBUF+41 (39) , W8U F*40 REMAINING FIELD 
FSWRITE PSCB=FILEID 
LTR 15,15 
BNZ WRTERR 

LA 7,1(7) ADD ONE TO RECORD COUNT 

SPACE 

*********************************************************************** 
* * 
* NOW CLOSE THE FILE. * 


UNL00560 
UNLOO 570 
UNL00 580 
UNLOO 590 
UNL00600 
UNL00610 
UNLOO 620 
UNL00630 
UNL0064U 
UNL00650 
UNL00660 
UNLOO 670 
UNLOO 680 
UNL00690 
UNL00700 
UNL00710 
U NL00720 
UNL00730 
UNL00 740 
U NL00750 
UNLOO 760 
[JNL00770 
IJNL00780 
UNLOO 79J 
UNL00800 
UNL00810 
UNLOO 820 
UNLOO 830 
UNL00840 
UNLOO 850 
UNLOO 860 
UNL00870 
UNL0088J 
U NL 00 890 
UNL00900 
UNL00910 
UNLOO 920 
UNL00930 
UNLOO 940 
UNL0095J 
U NL00960 
UNLOO 970 
UNLOO 980 
UNL00990 
UNL0100J 
UNL01 010 
U NL0 1 020 
UNL01 030 
UNL01 040 
UNL01 050 
UNL0106J 
UNL01 070 
UNL01080 
UNL01090 
UNL01 100 
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FILE: UNLOADSS ASSEMBLE A 


OHIO UNIVERSITY AVIONICS ENGINEERING CENTS 


* * UML011U 

******************** ************** ************************************* UNL01 120 


INVDISK 


WRTERR 


PARMERR 


SPACE UNL0113I) 

FSCLOSE FSCB=FILEID CLOSE THE FILE UNL01 140 

LINEDIT TEXT= * RECORDS WRITTEN TO FILE.*, *UNL01 150 

SUB= (DEC , (7) ) ,DOT*NO,RBNT=NO UNL01 160 

SLR 15.15 CLEAR RETURN CODE UNL01170 

BR 14 GO TO CMS TJNL01 180 

EJECT UNL01190 

LINEDIT TEXT= ' DMSULD054E INCOMPLETE FILEID SPECIFIED.*. *UNL01200 

DISP=ERRMSG, DOT=NO, RENT=NO UNL01 210 

LA 15,24 RETURN CODE UNL01220 

BR 14 BACK TO CMS UML01230 

SPACE UNL01240 

LR 2,15 SAVE RETURN CODE UNL01250 

LINEDIT TEXT=* DMSULD069E DISK »*..'• NOT ACCESSED.*, *UNL01260 

SUB= (CHARA, FILEID+24) , DISP=ERRMSG , DOT=NO, RENT=NO UNLO 1 270 

LR 15,2 GET RETURN CODE UNL0128J 

BR 14 BACK TO CMS UNL01290 

SPACE UNL01300 

LR 2,15 SAVE RETURN CODE UNL01310 

LINEDIT TEXT=*DMSULD105S ERROR •».." WRITING FILE »• *UNL01320 

ON DISK. • ,SUB= (DEC, (2) , CH AR A , FIL EID*8 ,C H AR *U NL0 1 333 


A, FILEID* 16, CHAR A, FILEID* 24) ,DISP=ERRNSG, 

DOT=NO,RENT=NO 

15,100 RETURN CODE 

14 RACK TO CMS 


LA 15,100 RETURN CODE 

BR 14 RACK TO CMS 

SPACE 

LR 2,1 SAVE PARM ADDRESS 

LINEDIT TEXT=* DMSULD070E INVALID PARAMETER •».... 

SUB= (CHARA, (2) ) ,DISP=ER RMSG, DOT=NO, RENT=NO 
LA 15,24 RETURN CODE 

BR 14 BACK TO CMS 

SPACE 

LR 2,1 SAVE OPTION ADDRESS 

LINEDIT TEXT=*DMSTJLD003E INVALID OPTION *• 

SUB= ( CHARA , (2) ) ,DISF=ER RMSG, DOT=NO,RENT=NO 


♦UNL01 340 
UNL01 350 
UHL01 360 
UNL01 370 
UNL01 380 
UNL01 390 
♦UNL01 400 
UNL01 410 
ONL01 420 
UNL01 430 
M NL01 440 
UNLO 1 450 
♦UNL01 460 
TJNL01 470 



LA 

15,24 RETURN CODE 

UNLO 1480 


BR 

14 BACK TO CMS 

UNL01 490 


EJECT 


UNL01 500 


DS 

OD 

UNL01 510 

F36 

DC 

F * 36 * 

UNL01 520 

IBUFRO 

DC 

AL4 (IBUF+79) 

UNL01 530 

OPTREP 

DC 

CL8 * R EP • 

IT NL01 540 

FILEID 

FSCB 

«* * * • , BUFFER=WRUF,ESIZE=80 

UNL01 550 

TEMP 

DS 

XL8 

UNL01 560 


DC 

C’ 1 * 

UNL01 570 

TEMP 1 

DS 

XL5 

UNL01 580 

FLAGS 

DS 

XL1 

UNL01 590 

WBOF 

DS 

XL80 

UNL01 600 

IBUF 

DS 

XL132 

UNLO 1 610 

TBL 

DC 

X’OO* 

UNL01 620 


ORG 

TBL+8* 16 

UNL01 630 


DC 

XL1 6 ' 00 FAFBFCFDFEFPF0F0FQF0F0F0F0F0F0 ' 

UNLO 16 40 


ORG 

TBL*1 1*16 

UNLO 16 50 


DC 

XL 16* F0F1F2F3F4F5F6F7F8F9F0F0F0F0F0F0 * 

UNLO 1660 

TRTBL 

EQU 

TBL-X *40* 

UNL01 670 


END 

UNLO A D $S 

UNL01 680 
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TECHNICAL MEMORANDUM OU NASA 78 


LORAN-C PLOTTING PROGRAM FOR PLOTTING LINES 
OF POSITION ON STANDARD CHARTS 


The Loran-C plotting program was designed to plot Loran-C lines of 
position on any standard chart and is used in the data -col lection 
system currently in use at Ohio University Avionics Engineering 
Center. 


James P. Roman 

Avionics Engineering Center 
Department of Electrical Engineering 
Ohio University 
Athens, Ohio 45701 




I. 


INTRODUCTION 


The NASA Tri-University program at Ohio University is currently involved in 
the development of a low-cost Loran-C navigation receiver for use in general aviation 
aircraft. This paper describes a set of programs designed to be run on the IBM System/370 
computer at Ohio University. These programs are used to plot Loran-C lines of position 
(LOP) on any common map or standard aviation sectional chart. The Loran-C plotting 
program JRPLOT FORTRAN uses a standard Cal comp-compatible plotting subroutine 
package for the Hewlett-Packard 7203A graphic plotter. 

This paper gives a general description of the features of the Loran-C plotting 
program. This program involves a simple add/subtract method to calculate the LOP. 

Refer to Figure 1 . Included is a description on how to use the program and some methods 
of operation . 


II. FEATURES OF THE LORAN-C PLOTTING PROGRAM 

The program will accommodate any scale of map desired. (Note: the larger the 
scale of the map the more distortion will occur.) The program was designed for standard 
aviation sectional charts; any larger scale than 1 :500,000 is not recommended. 

Plotting may be done on any size chart within the limitations of the Hewlett- 
Packard 7203A graphic plotter (10" high by 15" wide). 

Four station pairs are calculated in the prograrrfs execution where, for the 9960 
chain: (Master) control for W, X, Y, and Z is Seneca, NY 42° 42 1 50.6" N and 

76° 49' 33. 9" W 

Block address 1 is the W-pair 

Caribou, ME 46° 48' 27.2" N and 67° 55' 37. 7" W 

Block address 2 is the X-pair 

Nantucket, AAA 41 ° 15* 1 1 .9" N and 69° 58* 39.1"W 
Block address 3 is the Y-pair 

Carolina Beach, NC 34° 03« 46.0"N and 77° 54' 46.8"W 

and block address 4 is the Z-pair 

Dana, IN 39° 51 ' 07.5" N and 87° 29' 12.1"W 

For best results only plot two LOP sets on a single chart. 

The time difference for each line of position is placed to the top or side of the 
chart, depending on the angle of the LOP, along with the station pair identifier. 
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III. ACCURACY 


There are two sources of error in the system. Although the latitude/longitude 
conversion subprogram takes into account the curvature of the earth, the plotting 
routine is purely linear, therefore, the larger the scale of the chart the greater the 
percentage of error. The error for large scale sectionals is approximately 1 10 psec. 
The error for small scale geographic survey charts is 1 5 p sec. The other source of 
error is in measurement of the parameters listed below. These errors can be reduced 
by methods listed in Section IV depending on the accuracy the programmer wishes to 
achieve. 


IV. EXECUTION 

1 . The center point of the chart must be measured as accurately as possible. 
Then, the latitude and longitude must be taken from that point. The proper form for 
entering into the computer is: 

integer degrees integer minutes floating point seconds 

/ / 7 - / 

Wspace^space^.* 

2. Then measure the number of inches per degree of latitude and measure 
the number of inches per degree of longitude. See Figure 2. 

3. Then enter the increment that the lines of position should be spaced apart. 
Recommended increments are: 

50.0 for standard sectional charts 
2.0 for geographic survey charts. 

4. Enter the actual dimensions of the chart to be plotted (see Figure 2). The 
chart should be no larger than 11x15 inches and no smaller than 5x5 inches. These are 
the practical limitations of the Hewlett-Packard 7203A plotter. When the program is 
finished executing, the chains will be located as follows: 

Block Address Station Pair 

1 W 

2 X 

3 Y 

4 Z 

It is recommended to plot all four chains on a blank sheet of paper the same size as the 
chart. 
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V. 


IMPROVING ACCURACY 


This procedure should be necessary only when there is a need for extreme accuracy. 
To make the Loran-C plotting program more accurate, three time difference positions are 
needed before execution of the program: the time difference of the center point, the 

time differences for a point left of center, and the time differences for a point right of 
center. Then it is a simple matter to align the grid with the known time-difference positions. 
By adding or subtracting from the latitude position of the center point, the grid will shift 
north or south respectively. By adding or subtracting from the longitude center point 
position, the grid will shift west or east respectively. Another method of adjusting the grid 
would be to add or subtract from the inches per degree parameters. Then the grid may be 
expanded or contracted respectively. 


VI. SUMMARY 

The Loran-C plotting program is a system of plotting routines and conversion 
subprograms. The program is designed to accommodate a wide range of mapping needs. 
The program may be easily modified to meet the specific needs of the current experiment. 
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Figure 1 . Flow Chart. 
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Figure 1 . (Continued). 
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APPENDIX A. Listing of JRPLOT Program 


P( |C; J1PLOT FORTRAN C OHT 0 UNIVF® ST TY AVIONICS ENGTNFFR TNG CFNTFR 


rr. cr.crrrrccr.cc plotting ioran-c curvfs on sfcttonal fharts cc.rcrcr 

n|MFNS ION T0!4),TH! 4),0RP0S! 7) , POS! 7) 

DIMENSION TOCO! 7),X?(?,?00) ,Y2!2,?00> 
n T MENS I ON RIIF ! 1 I 
PH I R , GAMP 
rnWMHM PH I R, GAMP 

COMMON/CHAIND/OFL (?) , A «? , Aft , AO (81 , OM { R ) , C S f 8 > 

INPUT SECTION 


V/RTTF (4, 1 1 ) 

11 format (•+*,* fmtfo lATiTUDc OF CFNTPP POINT r OPM H*« UU ifV.H •) 
PFAD(5,4) !OFG,M!N,SFC 
WP T Tr I A, 12 J 

1 ? format! • • , • fnTFR LONGITIJOF OF r ENT FP point*) 

RFAO! * 5,4 ) I 0 FG 1 .MINI , SFC 1 
4 FORMAT! 18 , IX, 17 , IX, FO.l) 

WR T TF ( 4 , 10 ) 


JRPOOOl 0 

jrpdoo?o 

JRPOOOOO 

JRP00040 

jppoooso 
JRPOHOAO 
JRP00070 
JRPO00R0 
JRPOOOOO 
JRPOOIOO 
JRPOOl 10 
JRP00120 
JRPOOIOO 
JRPOOl 40 
JRPOOISO 
IRP00140 


10 cpRMAT!* » f »FMTFR NUMOER OF INCHES PER DFGRFF L AT I TUOF * ) 

RFAO( S, S ) A 
WPTTF(4, 14) 

FORMAT! • • , • FNTFR NUMRFR OF INCHES P p P OFGRFF LOMGITUOF*) 

RFAD(5,S)0 
5 FORMATIF7.4) 

WR I TF ( ft , l S I 

IS FORMAT! •+*,* INPUT I NCRFMFNT RFTWFFM L.O.P..*) 

RFAO! 5, S) x INC 
X I NO 7= X INC 

IF! XINC.GT. 10.0)XINC7=10.0 


R= R* (- 1 . 0 ) 

YTP = inFCMMIN+SFC/ftO.O) /60.0 
X 1°! = I0FG1 ♦(MINI ♦SFCl/ftO.O) /ftO.O 


CA! L PLOTS ( OUF , 1 , 11 ) 

MAP DIMENSIONS 


W° I TF (6,1) 

l FORMAT! «+*,' FNTFR HORIZONTAL DIMENSION*) 

PFAD!S,?)XLAT 
? FORMAT! Fft. 7) 

WR I T c ( ft » 8 ) 

0 format! •«■*, 'FNTFR VERTICAL DIMENSION* ) 

R F AD ( 5 , ’ ) XI 0NG 
XLAT=XLAT-0.07S 


JRPOOl 70 
JRPOOl RO 
JRPOOIOO 
JRP00700 
JRP00210 
JRP00770 
JRP00780 
JRP0OP40 
JPP007SO 
JRP00740 
JRP00770 
JRP007R0 
JRP00790 
JPR00300 
JRPOOOIO 
JRP00020 
JPPOORR0 
J RP00040 
JROOPOSO 
J RP 008ftD 
J R P 00 8 70 
JRPOOORO 
JRPOOIOO 
JRPOD4DO 


X10NG=XL0NG-1 .0 


JPP00410 


PH I R= ! IOFG*-! MIN+SFC /60.) /60. ) *8. I 4 1 ft o? ft 5/ 1 80. 

GA MR= l IDFG1MM!Nl*SFCl/60. ) /60. 1*0. I 4 l 50265/ IRO. 
YCP = XL 4T/2 . 0 


XF p 1=XL0NG/ ?• 0 


rr ft 


cnNVFPSlON TO time miff. ccccrrccccccccrr rrrrrr 


CAM TDS(TO) 

opodS < l ) = I OE0+! MIN«-SFC/ftO. J /60. 

ORP0S ! ?)=IOFG! +IMTN 1+SFCl /60. ) /ft O. 


JRPD0470 
JP.P00480 
J RPD044D 
JRP00450 
JR P 004ft 0 
JRP00470 
JPP0D48D 
J RP00490 
JRP00500 


WR TTElfi, 10) IDFG , MI N , SEC, I DEG l , M IN! , S Ff l 
10 format!* ' , • L A T I TUOF= *,I?,' ',17,* » , FF.2 , »LONGITUDE= » , T2 
>,* *,T7,* ' , FF . 2 ) 

IF! A.GT.70.0IGOTO 45 
0O 40 K= 1,4,1 


JRP00S1 0 
JRP00570 
JRPO058O 
JRO0054O 
JRP00550 
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-jrpldt 


FORTRAN r 


OHIO UN ? VFR S I TY AVIONIGS FNGIMFFR TNG GFNTFP 


c TLF: 



v 1 = TP ( K ) / I. 0 

JRP00660 


K 1 =K1 * 1 0 

JRP00970 


r=TO(K)-Kl 

JRPO09RO 


lF(C-9.0)?0,30,30 

JRP00S90 

?n 

TP(K1=TD1K1-C 

JRP00600 


goto 40 

JRP00610 

in 

Tn(KI-TO(KU(lO.O-CI 

JRP00600 

An 

CPNTTNiJC 

JRP00630 

A 9 

WP. TTF(6,501 (TO(KI , K = l ,4 ) 

JRP006A0 

sn 

format ( • » , »TOS=» , APIS.?) 

JRPO06S0 


nn 0999 17=1,9,1 

JPP00660 

rrc r rrrrrrrrrr.rrr.rrrrrr.frccrcrrrr.rcc<'r.rrrrrrrrrr.rr,r.rrrrr.rr.cr.rrrrrr.rrr.r. 

JRPO0670 

rrrrcrrrrrrrrrcrrr. data ali nc at i on rcrrrccrrrrrrrr.rrrrrcrccrrcrr r.rrrrrrr jopooaro 


IF( T7.0F.9 JCOTn 3000 

IPR006O0 


J 7=0 

J RP 00700 


I F f 17-0)90, 100, 110 

JRP0071 0 

on 

1 = 1 

JRP00700 


J = P 

JRP00730 


I 0=1 

| R P 00 7 AO 


) 0 =o 

JRPOO790 


PALL OATAWX 

,J R P 00760 


OPTO 130 

JRP00770 

l nn 

I = o 

JPP 00780 


J = 1 

J RP 00790 


T ?= 0 

JPP00800 


00=1 

J PP 008 1 0 


P M. L oat AW X 

JRPOORPO 


GOTO 130 

J R P 008 3 0 

110 

TFM7.F0.A)G0TP 100 

JPP008A0 


T = 3 

JRP00890 


J=o 

JRP00860 


r ■>=? 

JRP00870 


J 0= 1 

JRP00R30 


'■All. PATAXY 

JPP00890 


GOTO 130 

JRP00900 

l?o 

T = A 

JP D 009 1 0 


J = 3 

JRP009P0 


I o= ? 

JRO00930 


J0=1 

JPP009A0 


P.AI L 0ATAY7 

JPP00990 

1 in 

T F< ORPOSl?) .GT.76.0)PALl NFWOAT ( 17 , I , T ?, J, J ? ) 

JPP00960 


tp=o 

JRP00970 


I o = 0 

J RP 00980 


T H ( T?)=TO( I ) 

.IRP00990 


TH( J?)=TDl J) 

J RP 01 000 


T PFA|= 3 

JRP01010 

rrrr.r.rrrcrr.ccrcrr.r.rrrrr.crccc.cr.ccrcr.ccccr.c.crccrr.rrrccrrcccccr:cccrr.cc''r.r.r.rjP'>oio?o 

rrrrrcrrrrrcrrcrrrrcrrrrrrr plotting sfgttpn i rr.rrrrrrrrrrcrrrrrrrrrr.r. 

rj poo 10 30 

1 AO 

r ALL rnNVFRl TH, POS, DROPS, I FR RPR , VC P , XC » 1 , Y 1 P , X T P 1 , X , Y, 13 , A ) 

jRPOl 060 


rr=ir+i 

JRPOIOSO 


YOl l, I r ) = Y 

JRP01060 


X?( i, JC)=X 

JPP01070 


I 3 = 0 

J RP 0 108 0 


PALL TFSTl XL at, XLDNG.X, V, T3, IFRRPR) 

JRPOl 090 


f F( T3.0E.I )GOTn 160 

JRPOl 100 


101 


r ? i_ c : 


JPPLDT FORTRAN C 


OHIO UNIVERSITY AVIONICS FNSINFFRINO CENTER 


TALI PLOT ( X, Y, I PEN ) 

TH( J 2 I =TH( J 2 ) +XINC 2 
T PFN =2 

in if( 11-11140,900,000 

crrcrcrrccrr cccrcccc.r.ccc. plotting, section 2 


JRPOlllO 
JRPOI 120 
JRPOlllO 
JRPOI 1 4-0 

ccccccrcccccccccrrccrrcrrrccjRpni no 

JRPOI l 60 


ooo n=o 

npROS Cl.) = 1 OFG+( MINE'S rr/6 0.0 1/60.0 
npofJM 2)=IDFG1 +(MINI+SFC1 /60.0) /6 0. 0 
TH( J2)=T0( J) 

T P f N-s 3 

OOTO 91.*' 

910 TH(J?I=TH( J21 — XTNC2 

oil gall convfr < th, pns, opphs, t fprgip , ytp, yr p i , yt p, xt pi. , x,y,i , a > 

T n= tp + i 
Y 2 ( 2 , T 1 )=Y 

x 2 (?,tp)=x 
I 1= 0 

GALL TFSTI XI AT, XI ONC, ,X,Y,IT, IFPP0P) 

I F ( T 3 .GF .1 1 GOTO qio 
GAL!. PLOT! X, Y, IPFN) 

r pfn= ? 

ooo IF( TO — I 1910, 1000, 1000 


J PRO 1 1 70 

jrpoi mo 
JRPOI 190 
J RP 01 200 
JRPOI ?1 0 
JPPOl 220 
JRPOI 210 
JRPOI 2*0 
JRPOI 210 
IRP 01260 
JRP 01270 
JRR 012 P 0 
JRPOI 290 
JRPOI 100 
JRPOI 110 
JRP 01120 


crGrfrGcccc.GrrrcrccGccrcr numbering section rrrrccrcrrrcccccrcccccrcccrr jrpoi 1*0 


1000 UFTGHT=0.1 

I F ( rr .LT.20I goto too*; 

OGITGt 1 006 

10 01 0F|.TAX=X2( 2, ( TR/?))-y?( 2 , ( 10/1)1 

nFLTAY = V?< 2, ( TR/2) I-Y2I ?, ( IP/1) ) 

I P= 2 

SOTO 1007 

i oo*> oeltax=x 2 ( i,iir/ 2 ))-x?ii,( ic/m 
nr i tay=y2( i, « ir/2) >-Y2( i, ( iG/i) > 

TP=1 

in07 AMr,LE=ATAM2( OF! TAY,OFl. TAX) 

A NO LF = A 1 S ( f ANGLE*lflO. ) /1.141 19 2 6D 

I F { ( ( ANGLF.LF.41. I.OR. f AM0LF.GF.111. ) ) . OR . ( ( ANGI E.GF. 115. I 
>. 4N0. (ANGLE. I F ., P 2 1 . 1 I I TP011=l 
I F ( ( ( ANGLE.GT.41. 1 . AMO. ( ANGIF.LT. 111. 1} .G1P . (< ANGI E . GT . ? 2 1. ) 
>. AMD. ( ANGLE.LT. 111. ) ) ) IPOS 1=2 
1000 MPI. ACF = 0 
ANGLF=1 .0 
mghap=o 

FNI)M=TH( T 2 1 

TRCOI 1 1=101+17 

IF! IPOS 1-2) l 02 o,I 01 1,1011 

mu vnnwN=i.o 

CALL NU9RFR c X2 ( 1 , IC ) , ( XL AT-XOOWN) , HFT GHT , FNUN, ANGLE , NPI. AC F ) 
CALL SYMBOL (1X2(1 ,10+0. 21, ( XLAT-l X0OWN-0.2 ) ) , HE I GHT, mCO{ l 1 
> , N C H A R ) 

GOTGl 1010 

1020 TAIL NUHRFR(0.?,Y?( 2,1 PI ,HF1GHT,FNUM, ANGI F.NPLACF) 

TAIL SYMpnLlO.fl, Y?f 2, IP) , HF I GHT, I PC DC 1 1 , ANGL E , NC H A R ) 

1070 TG=0 
TR=0 


JPPOl 110 
JRPOI 160 
JRP 01170 
JPP 0 L 1 B 0 
JRPOI 190 
J RP 01 400 
JRP 01410 
JRP 01420 
JRPOI 410 
JRP 01440 
JRPOI 410 
JRP 0 I 460 
JRPOI 470 
JRP 014 R 0 
jppOl 490 
JR o 0 L 100 
JRPOlllO 
JRPOI 120 
JRPOI 510 
JRPOI 140 
JRPOI 110 
JRPOI 160 
JRPOI 570 
JRP 011 P 0 
A NGL F JRP 01190 
JRPOI 600 
JRPOI 610 
JRP 01620 
JRPOl&IO 
JRP 01 640 
JRP 01610 
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FILF: JRPLOT 


FORTRAN C 


OHIO UNTVFRSITY AVIONICS FNG I NE C R I NO CFNTFR 


TF{ J7.FQ. noOTP 7010 JRP01660 

CFrrrcrccccrccccrrcr.rcccccccccccrccccccrrccrcr.rrr.rcrrcrccrc.crrrrrrrrrrrrjPP0ift70 


increment up 

1010 TH( I ? J = TH( T7)+X INC 
THI J->)=TD( J) 

TPEN=3 

CM L CONVER ( TH , POS, ORPOS , I FR ROR , YC P . XT PI , Y I P , X I P 1 »X,Y,R,A) 
13=0 

TAIL TFST ( XLAT, XLONG,X, Y, H, IFRROR ) 

IF( TO-l )140, 7000, 7000 

7000 th( T?)=TD( I > 

J7=l 

M=o 

?010 TH( I?) =TH( T7J-XTNC 
T H ( J?)=TO{ J) 

TPFN=3 

CALI C ONVER ( TH, POS, nPPOS, TER ROR, VC P,XC Pi ,YIP,XIP1,X,Y,R,A) 

r 0 = 0 

TAIL TFST ( *LAT,XLONG,X,Y, 13, IFF ROR) 

TF ( TO-lI 140, 7999, POOR 


r ** r 

7000 

ooon 


JPP016R0 
JRP01 690 
JPP01.700 
JRP0I710 
JRP01770 
JRP01730 
JRPOI 740 
JRP017S0 
JRPOI 740 
J RP 01 7 70 
J RPO] 780 
JRPOI 790 
JRPOI ROO 
JRPOI fl 1 0 
JRP01870 
JRPOI 800 
JRPOI 840 
JRPOI 880 
JRPOI 860 


CALI PLOTIO. 0,0. 0,-3) JRP018R0 

FAIL PL OT ( 0. 0, O. 0,999) JRO01R90 

STOP JPP01900 

FND JRPOI 910 

S IJ8 POUT INF TOS(TO) JRP019?0 

Fr.ccccrrrcrcrccccccrcr. oata fop timf oiff. conversion ccr.r.ccccccccrcccccjppoiooo 

DIMENSION TO ( 4 ) J R P 0 1. 9 4 0 


PFAI «8 "HI M, GAMM, PH I W.GAMW, PHI X, GAMX, PHI Y,GAMV, PH I Z ,GAM7 
OATA PH I m,GAMM/0. 745 5007 761 ,1.040870774/ 

DATA PHI W, GAMW/O. 81 69491 S90 ,1.185559303/ 

DATA PHI X, GAMX /O. 77 00060 971 , l.7?1 345098/ 

DATA PHI Y,GAMV/0. *>94*505733 8,1 .38984 0019/ 

OATA PH I 7, G A MZ/ 0.69 8 8 47640 9, 1.82 69 7 8009/ 

OATA C0RW/1.379724F4/ 

OATA CDOX/ 7. 696991 F4/ 

OA^A C0BY/4.777161F4/ 

DATA COR7/8. 716208F4/ 

TOM=ARC (PHIM.GAMM) 

TO ( 1 )=r.ORW-*-ARC (PHIW.GAMWI— TOM 
TO | ? )=CORX+ARC ( PHTX, GAMX)— TDM 
TO(O)=r08Y4-APC(PHIY, GAMY) -TOM 
T0(4) = CDRZ+ARC( PH 1 7 , G AM Z ) -TOM 
RETURN 


JRP01980 
JRPOI 960 
JRPOI 970 
JPP01980 
JRPOI 990 
JRP07000 
JRD07010 
JRP07070 
J RP 07030 
J RP 07 04 0 
JRP07080 
J RP 07060 
JRP07070 
JP.P070RO 
J RP 07090 
JRP071 00 


r.rrrc.crcrcc.rcccrccccc testing RnuTiNF rrccrrcrcrcccccccrccrrrccrrc.r.rrrrrjRpo?i 70 

FNO JRP97130 

S URP n U T I NE TFSTIXLAT, XIONG, X,Y, TO, TERROR) JRP07140 


TFIY.GF.IXLAT-l .5) ) 13=1 

IF1Y.LE. 1.078) T 3= l 

IF ( X.GE. ( XL PNG— 1 .0) 113=1 

IFTX.LF. 1 .8)10=1 

TFI IERROP.F0.-l ) T3 = 2 

RFTIJRM 


JRP07180 
JRP071 60 
J RPO? 1 70 
JRP071 80 
JRP07190 
J RPO 7 700 
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TTir; JRP1.0T F DR TO AM C OHIO UNIVFRSITY AVIONICS r MO I MFFP TNG CP NT PR 

r.rrrcrrr,ccccccccccccccccr.cr.cccc.r.rrccrcrrrc.cccrcccccccccccrccccccccrccrcc.iRPO??io 
crcccrrrrc.rrrcrrcccrc.cr.ccrc oata fop cat. , lonc. conversion rrrr rrfrrrrr jrpo? 2 ?o 


F*>n 

JRP077T0 

SUBROUTINE OATAWX 

JRP07740 

f PMMON/f HA I ND/OFL ( 2 ) , A 5 , Afi , A 0( fl 1 , OM ( 51 , r S C 8 1 

JRP07750 

OFl (I 1=1 l. GET 

,IP°0?740 

net f ?) = 75. OFT 

JRP02770 

A 5= 7 . 1.77S404F4 

JPP07750 

A4=7. 170T751E4 

JRP07790 

AD ( 1 } = 47 .0 

J«»0?TOO 

A 0 ( 7)=74.0 

JRP07T10 

Tl=44.0 

JRP07T70 

AO( 41=47.0 

JPP07TT0 

APf 51 =47.0 

JRP07T40 

*0(41=74.0 

JRP07T50 

AOf 7)=41 .0 

JRP07T4O 

A0( 51 =40.0 

JRP07T70 

OM( l) = 47.0 

JRP07T80 

OM( 71=49.0 

JPP07T90 

OM ( T 1 =45 • 0 

JPP07400 

PM( 4)=5*.0 

1RP0741 0 

HMf S)=4?.0 

J PP 07470 

0M( 41=49.0 

JRP074TO 

ON ( 71=19.0 

JRP07440 

OM( 01=50.0 

JRP07450 

f S( 11=50.47 

JRP07440 

f SC 71=04.44 

JRP07470 

fS(T1=77.B4 

JRPO7480 

fS(4)=T9. 16 

JRP07490 

CS( 51=50.47 

JRP07500 

CS( 41=04.44 

JRP07510 

CS(71=I 1 .95 

JPP07570 

CS( 51=40.51 

JRP075T0 

p FT UR N 

IRP07540 

FMO 

JRP07550 

SIJBROUTTNF DATAXY 

JRP07540 

m^MON/CHA I NO/ OFL (7>,A5 f A6,AO(R),0‘*C8),rS(F1 

JRP07570 

OFL ( | 1=75. OFT 

JRP075R0 

PEL f 71 = T9.0FT 

JRP07590 

A R = 7 . I 75740F4 

JRP07400 

A4=7. 1 71045F4 

JRP07610 

AO ( 1.1=47,0 

JRP07470 

A0( ?1 =74.0 

JRP076T0 

A 0( T 1 =41 . 0 

.( RP 07640 

AO( 41=49.0 

JRPO7650 

*0(51=47.0 

J pp 07440 

A0( 41=76.0 

JRP07470 

*0(71=T4.0 

JPP07450 

AO(81=77.0 

JRP07490 

OMI 11=47.0 

JPP07T00 

nM( ?)=49.0 

JRP07710 

OM ( 71=1 5.0 

JPPO7770 

0M( 4 1 = 55.0 

JRP077TQ 

nM ( 51=47.0 

JRP07740 

0M( 41=49.0 

JPP07750 
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f ft p : jrplot 


r r)P TR AN 


r. 


nHTT MN1VFRSTTY HVinNKS C N 9 TMFFP 1NG CFNTFR 


0Y(7)=3.0 

JRP02760 

hmi «) = •;/♦. o 

JRP02770 

r<!( l 1=59.47 

JRP02780 

FM 21=34.44 

JRO02790 

rsf 3>*ll-*>** 

JRP02800 

r. M*1 =40. 51 

JRP02810 

05(51=50.47 

JRR02820 

I’M 61=33.44 

.IR'»02830 

05(71=45.61 

JRP02840 

F5( 81=47. ?0 

JRP02850 

ppriiPN 

J R P 07 860 

FNO 

JPP02870 

sunnniiT TNF oatay7 

JRP07880 

COMMON/CHA f M 0/0FL (2 » , A5, A6 , AO( 8 1 ,OM (8 1 ,C5 ( 31 

JPP07890 

nF( ( 1 1=39. OF3 

JRP02900 

nFI. ( 2 ) =54. 053 

JRP02910 

A5=2. 128260F4 

JRP07920 

A 6 = 2. 121 04 5F4 

JRP02930 

AD( 1 1=47.0 

JRP 02940 

A0(2>=76.0 

JRPO2950 

AD( 31 =34.0 

J PP02960 

A0( 4) =77.0 

JPPOP970 

AO( 5}=47.o 

JPPO7980 

A 0 ( 61=76.0 

JPP02990 

A0( 7} =39.0 

JRP 03000 

A 0 ( a ) = B 7 . o 

J R. d 0 3 0 1 0 

DM( 1 ) =4? . 0 

JRP 03070 

OM( ?)=49.0 

JRP03030 

OU( 31=03.0 

JRP 03 040 

OM< 41 = 54.0 

JRP03050 

OM( 51=42.0 

J RP03060 

6) =40.0 

JRP03O70 

OMf 71=51 .0 

J RP 03080 

OVf 81=79.0 

JRP03090 

rs» 1 1=50.47 

J RP 03 1 00 

rs ( 7! = 34. 4 /, 

JPP031 10 

rs( 31=45.95 

JRP071 20 

C5(41 =46.7 A 

JRP031 30 

r 5 ( 51=50.47 

JP.P03140 

r<;(6>=34.44 

JRP031 50 

CM 71=07.46 

JRP03160 

05(81=12.14 

JRP031 70 

PFTURN 

JRP03 1 80 

FNn 

JRP03190 

5U8P n lJT (NF NEWOATt T7,T,T7,J,J71 

JRP03200 

rOMMON/CHA! N070FL (21 , AS , A6 , ADI P 1 ,0M (8 1 , C5 ( 8 » 

JRP03210 

I F ( 17-2110,70,30 

JPP03220 

1 = 1 

JRP03730 

J=? 

JRP03240 

T 7=1 

JRP03250 

.12 = 2 

JRP 03260 

CALL OATAWX 

JRP03770 

00T040 

JOP03780 

r =2 

JRP 03 290 

J = 3 

J RP03300 
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e II C; JPPLHT pORTP AM C 


Ohio untvfrsity avidnics fnotnffp ing center 


12=1 

J?=? 

'■*1.1 OAT 4W 
r.riTn 40 

oo >p t rr. fo.4)goto os 
T - 0 
J = 4 
17=1 
J 2=2 

nil OATAY7 
OPTO 40 
0 K T =4 

1 = 0 
T 7=7 
J 7=1 

TAIL OAT AY Z 
4 O PPTtJRN 

rrcrcrrcrcrccrccrccc conversion routine rr.rrr.rr.rccccccrcccr.rccrcrr.rcr.ee 

FMO 

SijprPUT IMF rnMVFP(TH, n OS,PRPnS, TFPRPP T vrp f YCP1 » YTR, XIR1 , X, Y,R f A) 
or MENS TOM TH ( 4) , POS ( ? ) , DRPOR ( 7) 

C A I 1 TPOOSl TH, PDS ,ORP0 S» terror ) 
nn in K=l,7,i 
10 0RPOS{K)=P0S(K ) 

v =( (dost ?)-* tpi )*R)> xrpi 
y=( (POSI i)-ytp)*a)*ycp 

70 PFT'.JPN 


JRP00010 

JRP00070 

JRPOOOOO 

JRP00040 

JRPOOOSO 

JRP00060 

JRP00070 

JRPOOOSO 

JRPOOOOO 

JRO00400 

JPP00410 

JRP00470 

JRP00400 

JRP00440 

JRP004SO 

JRP00460 

JRP00470 

rrrrccc':rr r rrrrcrrrjRooo4so 
crcccr.r.ccrcrr.rccccr jr»oo4ro 

JRPOOSOO 

YTP.XIPl, X,V,R f A) JPPOOSTO 

JRPOOS ?0 
JRPOOSOO 
JRPOOS40 
JRPOOSSO 
JRPOOS60 
JRP00S70 
JRPOOSSO 

rrrrc.crcrcrcrrrcrr jrpoosqo 


FMO 


JRP00600 
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APPENDIX B. Listing of TDPOS Program 


F 

n p: 

mans FOPTRAN r OHIO UNTVFRS ITY AVIONICS 

FNCiT NEFP TNG C FNTFR 



SOP ROUT INF TOPOS ( TH , POS , op POS , I FPROP ) 

Tnpoooi 0 



OIMTNSION lOSVIP) tOSVl? J t ANGf 8) , AOf RI , CS( fl) , PDSI 7) , 

7. IWO( 7) , FTOPOO07O 



: W'M 71 ,0^1 ( ? I » T OP ( ? ) , OROI ?) ,nP.H( ?» ,TH( 7) , A I II ) ,B( 111 ,C( 111 

,0( m ,FITOPOOOOO 



! M ) iff (1 1 ) *™(? ) , RIF>*(? » , REOCL ( ? ) ,RAnp {?) , BFTA I?) .OMGl 7 J, 

T WO ( 7 ) , T DTOo 0004 0 



: p 0 f 7) ,OPP0S( ?> 

TOPOOOSO 

r 

- 


TOP 00060 



OAT A A I /74.030S/,A7/-0.4O7S«/, AS/3. 4677AF-1/,B1 /0.S104P3/ 

,n?/-0.OiT0P00070 



: 1 40?/,ns/0.001 760/, RO/l. 74S3?9F-7/,PM/?.9nRRfl7F-4/, RS/4.R481 77F-6/TOP000R0 



: , PI /S. 1 41 *9?/, A4/ 2. 9969 1 2E7/ 

TOP00O90 

r 

- 


TOPOOIOO 



rOMMOM/CHATNO/nFL ,AS,46,AO,OM,CS 

TOP 00 1 1 0 

r 

- 


top on i 70 

'■* 

- OF0IN t]mf DfFppppMfp TO POSITION f ONV FP S T nn. 

TOP00I 00 

r 

- 


TOPO0140 



on 1 1=1,7 

TITPOO 1 SO 



TOPOm =0RPOS ( T ) 

T0P001 60 



op o ( t ) = innot T ) 

TOPOO 1 70 


1 

POM( I ) = ( or POS ( I )-ORO( I ) ) *60.0 

TOPOOl PO 



TFDROR^l 

TOP 00 19 0 



00=1 4S*AS-A6*A6)/( AS*AS ) 

TOP007O0 



A 1 4=1 . 0-A6 / A c 

T OP 00 7 1 0 



A *0=1 1 .OfAl 4+Al4*A1 4) 

70P00770 



4M = I ASO-1 .0) 

TOPOOPOO 



AS?={ A l 6* A l 4 ) / 7 . 0 

TOP 00740 



A ^T = - AS 1/7.0 

T0PO07S0 



A s 6 = ( A 1 4* A 1 6 ) / 1 6 . O 

T QP 00 7 60 



A SS= ( A14*Al4)/«.0 

T OP 00 770 



ASft = A I 4* A I 4 

T0POO7P0 



ACTrfl'IA*! . ?S 

THP00790 



4SR=AS6/4.0 

T0P0OR00 



00 l 7 R K=l,fl 

TOPOOOIO 



TF( An(K I ) 1 ?/,, 1 76, 1 7 6 

TOPO007O 


174 

ANGI k ) =PD*AO< K )-RM*nM( K )-P S*GSf K ) 

TOPOOOOO 



GO TO 1 ?P 

T0P00040 


I ? 6 

AMO <K ) =0 0* AO I K ) +RM*OM(K » +»S*rS ( « ) 

T OP 000 S 0 


170 

PONT T Ml J r 

T OP 00 06 0 



A I ?= ( A N G ( l )-ANG(6 ) 4- AA'G ( 7 ) - A^'G ( 6 ) I 

T OP 00 0 70 



A 1 7 = ARS I A1 2) 

topooopo 



T F ( A 1 7-0.0000! ) 7, 7, A 

Tnpoonno 


7 

A I 1 =— I . 0 

T 0° 00400 



GO TP g 

T OP 004 l 0 


0 

A l 1 = 1 .0 

T0PO047O 

c 

- 


T0°00400 

r 

- APPROXIMATE POSITIONS A NO STATION C OORO INAT FS . 

TQP0044O 

r 

- 


TOP 004S0 


9 

f ( I ) = ANO< 1 ) 

TOP 004 6 0 



Ff?)=ANGI?) 

T 0° 00470 



rr ( n = ANGm 

TOR 004R0 



rri 2)= ANGI 4) 

TOP 00490 



r (0)=SIN(E(ll) 

TOP 00 SOO 



F ( 4 )=CPS ( F ( l 1 ) 

T OP 00 S l 0 



F(S)=F(R)/F(4) 

T0P00S70 



r f R > - ( F ( F l )*( I..0-A1.4) 

TOP00S00 



af7 = atani f < n ) ) 

TOP00S40 



F(6)=SIN( A67) 

TOP00SS0 
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r I | r 


TOROS 


F OR TR A M C 


nHTO UNIVERSITY AVIONT'S F^G l NFFP T MR TFMTFP 


r J7)=r n S( A6?) 

rr n) = siM(rr n n 

rnft)=cnsicri n ) 
rr(<?)=rr, ( 3 )/rrfA) 
rr< ft)=icrift) )+( i.o-aia> 
AA?=ATAM<rC( R) ) 
r r ( ft j IN ( Aft?) 
rr (7)='- f iS( A62) 

T = l 

r,n rn ftoo 

1 5 c(0) = AV5 

F ( J 0) = Aft* 

FI l l ) = Aft7 

no I? j = i,u 

a r ji = f( i) 

17 R(J)=CP(J) 

FI 1 )= ANOI ft) 

F ( •> ) = ANG( ft ) 
rr. ( 1 ) = A M G I 7 ) 
rr i ? ) = ANO t ft ) 

Fm=M'iiFiin 
F(A) = cns(F(i) ) 
r (R)=C( 3)/EI4) 
c(ft)x(F(S) )*U.0-A14) 
a ft 7= A T AN ( FI ftl ) 
r ( ft) = S INI Aft? 1 
r<7!=CnS(A6?) 
rr < 3) = s INI re 1 1 ) ) 
rr ( ft i =ms irr ( n ) 
rr. i«>)=rr. m/rr. (A) 
rr(R) = <rr<«>n*(l.O-AlA) 
Aft?=ATAN(Cr (ft) 1 

cr i ft) = s ini Aft?) 

rr ( 7)=cnsi Aft?) 

1=7 

rn TO ROO 
19 E ( 9 ) = A"* ft 
F ( 10) = A4ft 
F 1 1 | ) = A'*7 
OP ?1 .1 = 1,11 

n n = Fi.n 

71 ni.n = rr(j) 

TMt 1 ) = A( 101+AI 1 l) 

T’a t ? ) = r i i n ) +r ( l l ) 
nn ftft m=i,? 
r>FTAf M)=TMIN) 

1FOFI I M 1 = ft FT AIN) «-O c L I M ) 
ftft ftLFNIM)=ftFTAI M )<-RPOFI (M) 

I OS VI l ) =09999 
TOSVt ?) =99999 
1 TF»=0 

n? SOR=OROl 1 ) 4-ORMI 1 ) +Oftnt 7 ) FOR All 7) 
)E(SnP)ft3,ft4,ft3 
ftft on RO K=l,7 

I f( ORO f K )) 37, 3ft, 34 


TORflOftftO 

ropoos7o 

Toonoftfto 

roROOftoo 

TOP0O600 

toroomo 

rnpno 6 ?o 

T0P00ft30 
TOROOAftO 
TOPOOftftO 
T0° OOftftO 
TOP09A70 
TOPOOftftO 
TOPOOftPO 
r pP 00700 
T OP 0071 0 
T OP 00 770 
TOP 00700 
T OR 00 740 
TOPO07S0 
Tnp 007*0 
TOP00770 
TOR007RO 
TOR 00790 
TOROOPOO 
THPOOftlO 
TOROOftRO 
TOPO0R3O 
TOPOOftftO 
TOPOOftftO 
TOPOOftftO 
TOPO0R70 
TOPOOftftO 
TOR 00ft9O 
T OP 00900 
TOPOO<M O 
T OR 009 7 0 
TOP00930 
TOPOOQftO 
T OR 009ftO 
TOOOOQftO 
TOR 00 9 70 
T0P009 ftO 
TO™ 00090 
T0P01 000 
TOROTOl.0 
T ORO 10 7 0 
T0PO1 ORO 
T0P01 O'VO 
T0R01 OftO 
TORO) OftO 
mooi 070 
TORO) OftO 
TOP01.O90 
TOROl 1 00 
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I 


r- 1 1 r : TO^ns FORTRAN 


n Hin UNTVFPSTTY AVIONICS FNGTN^FP IMG rFNTFP 


G? P A0» (K )=40*0PO( K )-PM*OPM( K ) 

Gn TO ?Q 

G4 P AGP ( K ) = P O*0RO { K ) +P M*DRM ( K 1 
GO rnMTIMUF 

F ( l ) = P AO*> ( 1 > 

F ( ? i = PAOP (7) 

A?S = - I .0 

q 4 Fni*si‘iifiin 

f ( a )= cns i e i i i ) 

C{6)=F(G1/F{4) 

P(«) = (F(«5n*n .n-Ai 4) 

A'»? = 4TAN( F< S } ) 

F{6)=SfM{A6»l 
r(7i«r.GSC aa-m 
on Sf J = |,S 
sa rr(j»=G(j) 

1 = G 

Gn to son 

°0 r 1 = AGP 
r?- A44 
r T= AA *5 
r mi = A47 

on 4? .i=i f s 
9 7 rc(.J> = r(j> 

1=4 

r.n to poo 

05 G4=AGS 
r s = A A 4 

0 6= A A S 

f l 0 4* A 4 7 
on so J=1,A 
so r r i j i = p ( j i 

T = S 

on fn s oo 
SS f 7 — AGS 
rn= A44 
r<?=A4S 
r 107= A47 
rr< a» i ) s?,oo,sg 
ss rin=r7 

ru=rq 

r i ?=C9 
r i io=r i 07 

GO 63 J = 1,A 

43 re ( J ) = A » .1 ) 

1 = 6 

go rn son 

6 S r 7=AGS 

rp=A44 
T0= A4S 
r 107=A47 

r i g=th(?»-c ( loi-r c 1 1 )-r i oi «-r 104 -op i (?) 

r i7=r 1G*A6 

r i G = THf 1 )-M 10 )-A( l l)-r 11 0«-C. 10 7-GFI. ( 1 1 

r ?2=r 1S*A4 


moot i io 

T 0 p 01 170 
TGP0U30 
TOPOI 140 
TGP01 1 SO 
TGP0U60 
TGP01 1 70 
TOPOI ISO 
TG»01 190 
TOPOI 700 
TG°0 1710 
TGP01730 
TrtPOl 7G0 
TOOOl 740 
TG ra 01 7 SO 
TOP 01 260 
TOPOI ?70 
TOPO) 2 SO 
TOPOI 790 
TOPOI GOO 
TOPOI G 1 0 
TOPOI G 70 
TOPOIGGO 
TGP01G4O 
TOPOI GSO 
TOPOI GAO 
TOPOI G70 
TOPOI GSO 
TOPOI GOO 
TOPOI 400 
TOD01 41 0 
TOPOI 470 
TOPOI AGO 
TOPOI 440 

TOPOI 4S0 
TOPOI 460 
TOPOI 470 
TOPOI 4S0 

TOPOI 490 

ropot son 

TOPOI S10 
TOPOI S?0 
TOPOI SG 0 
TOPOI S40 
TOPOI SSO 
TOPOI 560 
TOPOI S70 
TOPGISSG 
TOP01S9O 
TOPOI AGO 
TOPOI 61 0 
TOP 01 6^0 
TOPOI AGO 
TOPOI 640 
TOO016S0 
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r ii. r 


Tr>on$ 


FORTRAN r 


OH 1 0 UNTVFRSITY AVIONirS FNGINFFRTNG CFNTFR 


no 4R40 Y=l,7 
W=nMG(V)/»r> 

Iwn(M)=w 
Fwn=TWO(M) 
nwn=w-Fwo 
FWM=0WO*8O.0 
FWM ( m ) = 4RS ( FW M I 

TF{FWM(M)-S9.99aS)4840,4Bl0,48l0 

4810 r WM ( M ) =0 . 0 

IF ( IWD(M) ) 4R20,4R30,4880 
4 8 0 0 TW r )( M ) = IW r >(M)-l 
on t n 4840 

4880 TWO(M)=TWO(M)+l 
4840 onNTINUF 
r,n T n OOO 

r* 

r - CM r H1 A T I ON OF INVFRSF 1/4° T API F S . 

r _ 

800 A so=-rr(?) 

A60=-F( 7) 
r88=AS9-460 
C 86=ARG((". 8S) 

10(086- 81)801, 808,802 
SO? M 4=8.0*81-086 
00 Tn 808 
SO! Al6=r86 

80S I F( A l 6 » 804 , 807, 806 
SOT A 1 4 = 0. 00000008 
S04 A17=8IM(A14) 

A 1 8 = rno ( A1 6) 
a ip=c ( 6 ) *rr ( 6 ) 

A ?o=f { 7 ) *TC ( T) 

4? 1 =419+470*418 

427=( f 4i7*rrf 7n**7 + ir.r( 6) *f< 7)-«m 6 >*rr c 7 )*a 18) **?» **o. s 
A28={ A70*A1 7)/A?8 
4?4=1 .0-478*473 
A.8S = AR8TN( A 2?) 

4?6=A?S*A7S 

427=1.0/48? 

A 7 8 = A ? 1 / A ? 8 
4 ?0 = A 24* A?4 

430=14 8 0*A?6)+A1O*1 481*477- 4 8 7* 476*427) 

481=424*1 4 83*478+ 44 3* 477*471+ 482*476*478) 

488=410*41 9* (~48?*4?1*A?8) 

A 8 8= A 29* ( 484* A?8+A84*A 71* A??-A8?*A?6* 478-488*422*1 A 71 **8)) 
434= 4)9*474* <48?* 476*477+467*477*421*471) 

4?S = ( 4 30 + A 3 1 +4 32 + 433 + 4 34 ) *4 6* 44 

486=1 A81*A?6)+A1 0*( - A87*A7?-41 4*A) 4*A?6*A77) 

4 87 = 4 ?4*(-A8 7* A 78+488*427* 4? 1*4 14*414*4? 4* 428) 
A3R=(A36+A37)*A?3+A«6 

489=STM( A38) 

A40 = C.nS 1 A88 ) 

A4t=irr(61*F{7)-A40*F{6)*rrf 7) )/{439*rr(7> ) 

T f ( 44 1 ) 810,809,810 
809 441=0.00000008 


T0P07710 

TOD02770 

T0P07730 

T0P02740 

TOP07780 

T0P07760 

TOP 02770 

TOP07780 

TOP02790 

TOP07300 

T0P0?3! 1.0 

TOP02870 

TOP0778O 

T0P0?340 

T0P07380 

TOi»07840 

T OP 07870 

TOP073RO 

TOO07390 

TOP074O0 

T OP O? 4 l 0 

T0P07470 

T0P07430 

T0°0?440 

T0P07480 

T OP 0 8/.60 

T OP 0? 470 

TOP07480 

T0P07490 

T0P07800 

T0O0781 O 

T0P0?8?0 

T 00 07 8 80 

Tr)t»o? S40 

T OP 0? 8 SO 

T0P07860 

TOPOOS70 

T OP 07 88 0 

TOPO?890 

TOP07400 

T OP 0 76 1 0 

T0?0?6?0 

TOP07680 

T0°0?640 

TOP0768O 

TOP08660 

T0P07670 

T0° O7680 

TOP07690 

T OP 0? 700 

T OP 0? 7 l 0 
TOr>0?7?0 
T0PO278O 
TOP0?740 
TOPO87S0 
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r Tl r : Tfiont; 


F no t P A N r 


nMjn IJMJVFOSITV AVIOMir < ; PMf, I NFFR I ^G FFMTFR 



C.?1=Ct-r 17 

TOPOl 660 


r?4=r4 

TOP01 67(3 


r?s=r7+r?j 

T0P016B0 


c?6=rio 

T 0° 01 690 


r?7=l ir?-r6)*(r?6-r?6>-Mr?3-r?4)Mri1-rB) ) 

TOPOl 700 


p?9=< (f?-r , i)*(ri?-ro) + (ri-r.6)*(C ( »-r n n 

T0P01710 


C70=r 27/c?o 

TOP01770 


r ?fl=( r?7-r •>A+ r 10*m-rA) 1 / (T 5-C21 

T OP 01770 


on m l 3o 

T0P01740 

6? 

rnsTHi?)-rnoi-fn n-rion-rio4-oci 

TOPOl 760 


r 1 7=r 17*A4 

TOOOl 760 


r 1 R = TM( 1 ) - A( 10 )-A< 1 U-r 1 07*-r. 1 0A-0F(. ( l 1 

TOPOl 770 


r?7 = f l R*A4 

TOPOl 7B0 


f ? 3=r i — r 1 7 

TOPOl 790 


F?4=f4 

TOPOl BOO 


r?6=r 7-C7? 

TOPOl BIO 


r 77 = r r?*{ r?6-r ?4 ) *r ?7* (F6-fb )+r.R*r?A-r 6*r ?6 ) 

TOPOl B?0 


r ?9= r r?* ifa-cp) «-F3* < r 8-rs i+f 6*c9-c b*F6 ) 

TOPOl fl 30 


nn = r ?T/r?q 

topoi B40 


r 7« = ( r?7-F2'**no* (cvcmi / (rs-r? i 

TOPOl B 60 

1 10 

T 7 l = ( AS*A4*f l . 0— A 10)1/11. 0— A 1 0* F ( 7 ) .*=> 

T OP 0 1 B60 


r37=| AA*AA)/(1 .0-MO*Fn|*Fnn**A.‘i 

TOPOl R7ri 


rv)*mn/ni j 

TOPOl RRO 


F74=(-F?B/( r 7?*F{4) ) 1 

TOPOl B90 


F( ll = Fl 1 )+F,3? 

TOP 01 900 


F(?1 = F(?)«-F.34 

TO" 01. 01.0 


l F( A7R ) I 77,99, 1 74 

TOPOl 070 

n? 

a?r = i . o 

TOP 01 97 0 


nn to r a 

T0O01940 

r _ 


T0°01 960 

r - rnfjVFOMON nnNF, RFTIJPN TO 01 ST ANFF-B CAR TNG RnUTINF. 

T0P01960 

r _ 


TOPOl 970 

OOO 

TF( TPSVll ) .MF. Twn(l) )GO TP 7717 

T0P019B0 


I r f I06V( 7) .NF. !WO( 7 ) )F,n TO 7717 

TOP01O90 


irtABSIOSVn )-FWM(l ) I.GT.O.DGO TP 7717 

T OP O’OOO 


TF1ARS(0SV(?)-FWM(?)).GT.0.1 ) GO TP 7717 

T0P070 I 0 


t pp. m » = t no ( n * i o 

TOP07070 


on R7Q 1=1,7 

T0P07070 


7 iwm I) =IW0( n 

T OP 0 ? 040 

° 3 9 

PPSC T 1 = 7 IW01 I ) + FWM1 n /60.0 

T0°0?060 


P FTIJON 

TDP07060 

r - 


T0P07070 

r - rnMTINUE ITFRATIPNS. 

TOPO70B0 

r __ 


TOP 07090 

7 7 1 ? 

np 771? *1=1,7 

TOP0? 100 


ppn(M)=9.r) 

T0P07110 


rv>M(M) = 9.o 

TOP 07 1 70 


PSV ( M ) = C WM 1 M ) 

T0P0?1 70 

771 ? 

106VCR) = TW0(M) 

T0P071 40 


TTFR=ITFR+l 

TOPO? 160 


1F1 1TFR.LT. V0O)Gn TO 8? 

TOPO? 1.60 


I frp.°p=-i 

T0P07170 


P FT UP M 

TOPO? 1 BO 

1 3* 

PHG m=F(ll 

T0P07190 


PMG(?)=F(?1 

T0P07200 
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run 


TOPOS 


FORTRAN 


r 


OHin DMIVFRSTTY AVIONICS FNOTNFFR IMP, CFNTFR 


510 A4R=1.0/A41 
A4T=ATAN( A4?) 

TF(r*5)51S, 514*514 

M4 T f f CRS- p?) 511 , 5l?,5l2 

511 TF ( A4I ) 5RO ,5*1,5?! 

SRO A 40 = P I +■ A41 

r.n to Sri 

FIR T F f A41) SI 7, SIR, SIP 
SIS iF(rRS + nr ) Sll , Sll , SI f- 
FI 6 TFIA4 1)517, SIR, SIR 
SI? A4R=ft-A40 
on TO F?1 

SIR 44R=R . 0*P I-A4R 
SRI »4R = A4R*r>! 

A4R = A4R-R . 0*P I 
T r ( 440) 5RR,5RO,5RR 
SR? A4R = A4R«-R.0*PI 
SRR A44-=S r M ( 441 ) 

A4S-rns r 441 ) 
446=415/1*00.144 
IF ( 446-1 00.0) SRF , S?6,5?6 
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AUTOMATIC GAIN CONTROL 


An automatic gain control has been designed and fabricated 
to operate with the Loran-C prototype receiver and data 
collection system currently in use at Ohio University. 


James P. Roman 

Avionics Engineering Center 
Department of Electrical Engineering 
Ohio University 
Athens, Ohio 45701 




INTRODUCTION 


The NASA Tri -University program at Ohio University is currently involved 
in the development of a low-cost Loran-C receiver for use in general aviation aircraft. 
An automatic gain control (AGC) has been designed and built to operate with the 
prototype Loran-C receiver. 

Since there are such extreme distances between Loran stations, the signal 
strengths coming into the user are at different magnitudes. It is advantageous to have 
the signal magnitudes equal; therefore, the automatic gain control was designed for 
the front end of the prototype Loran-C receiver. 


II. CIRCUIT DESCRIPTION 

The automatic gain control is a three -transistor circuit (see Figure 1) which 
requires a constant D.C. voltage of 8 volts. Tests conclude that this value may be in 
the range of 4 to 12 V without change in circuit performance. Transistors Q] and Q 2 
are cascaded to pass and amplify the input signal . The gain of Qi and Qj is controlled 
by Qg which itself is controlled by an external AGC voltage between Oand 8 volts D.C. 

The integrated circuit used is an RCA CA3028A, an 8-pin chip, which is a differential 
cascode amplifier designed for use in communications operating at frequencies from 
D.C. to 120 MHz. The integrated circuit has been balanced for AGC capabilities, and 
has a wide operating current range. The maximum input current at pins 1 and 5 is 0.1 mA. 

The absolute maximum dissipation at T a < 85° C is 450 mW. At T Q > 85° C the inte- 
grated circuit is derated linearly 5 mW/°C. The ambient temperature for operation is -55° C 
to +1 25° C and -65° C to +1 50° C for storage. 


III. TEST RESULTS 

1 . Gain Vs. Frequency and Phase Angle Vs. Frequency (Figure 2) 

The frequency response and phase angle is measured with the input signal 
voltage held constant at 50 mV, and is an average of all AGC voltages from 1 V to 8 V. 
The band width of the automatic gain control is 20 KHz to 2 MHz, with a standard 
deviation of no greater than 1 0.4 dB. The phase angle increases linearly at frequencies 
between 70 KHz to 120 KHz, from 0° to +21 .6° respectively. 

2 . Gain Vs. AGC Voltage (Figure 3) 

This test was performed at a constant frequency of 105.4 KHz and a constant 
input signal voltage of 50 mV. The gain is approximately -2 6 dB from 0 to +2.5 volts, 
between 1 and 2.5 volts the gain increases rapidly from -22 dB to +2.5 dB, between 3 and 8 
volts the gain increases from 2.5 dB to 25 dB. Distortion and gain compression occurs at 
8.4 volts AGC and loss of gain occurs at AGC voltages greater than 12 volts. 
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3. Gain Vs. Signal Voltage (Figure 4) 


This test was performed at a constant frequency of 100 KHz and is an average 
of AGC voltages from 3 to 8 volts. The operating recommended input voltages are between 
5mV and 75 mV. In this range there was a constant gain with a standard deviation of only 
± 0.109 dB. 


4. Power Dissipation and Circuit Resistance 

This test was performed at a frequency of 100 KHz with an input signal voltage 
of 50 mV. The circuit draws 50 mA constant current. The D.C. power dissipation 
was 59 mW (Max. Pq q dissipation 450 mW). A.C. current at pins 1 and 5 is 
approximately 61 .2lp A at 100 KHz. The A.C.power dissipation is approximately 
3.035 x 1(T^ watts with a power factor of .9 leading calculated error for A.C. power 
is 4.1%. 


Circuit Resistance 

Input Resistance Rj = 140 ft 

Output Resistance R^ =740 ^ 


IV. SUMMARY 

The automatic gain control was designed specifically to operate with the prototype 
Loran-C receiver and data collection system. The use of the automatic gain control is 
intended to eliminate error which occurs when signals are received at different magnitudes. 
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A LORAN-C PROTOTYPE NAVIGATION RECEIVER 
FOR GENERAL AVIATION 


The design, fabrication and evaluation of a proto- 
type Loran— C receiver is described. Hardware is 
complete, and microcomputer programming continues, 
for addition of area-navigation capability. The 
receiver is an envelope-processor, offering simpli- 
city of RF processor circuitry. 


Robert W. Lilley 
Daryl L. McCall 


Avionics Engineering Center 
Department of Electrical Engineering 
Ohio University 
Athens, Ohio 45701 


I. 


INTRODUCTION 


The Avionics Engineering Center, Ohio University Department of 
Electrical Engineering, has pursued the techniques required for 
Loran-C navigation with application to the general-aviation pilot. 

The goal has been to produce prototype equipment for flight evaluation 
which will provide enroute navigation in both latitude-longitude and 
rho-theta coordinates and to evaluate the non-precision approach capa- 
bilities of such equipment. 

For this prototype project, single-chain, master-dependent 
operation was chosen as a demonstration mode, with three stations 
tracked. The prototype hardware design will, however, permit cross- 
chain, master-independent navigation. The number of stations tracked 
simultaneously can be expanded. These extended operating modes are 
implemented through programming utilizing existing shared tracking- 
loop hardware as described below. 

The prototype Loran-C receiver has been flight-tested using a 
variety of flight paths, with and without simultaneous ground radar 
position data collection. Results are presented later in this paper; 
further flight evaluation is planned, and will be reported separately. 

The following sections describe major receiver elements, 
drawing upon the work and publications of project team members who 
contributed to the design. 

II. RECEIVER OVERVIEW 

As configured for laboratory and flight evaluation (Figure 1), 
the prototype Loran-C receiver utilizes an aircraft ADF sense antenna 
or similar unit, connected directly to the wide-band preamplifier/ 
coupler. Signal levels for the linear RF processor are stabilized by 
a commutated, sampled AGC element, under control of the receiver com- 
puter. 


The RF unit performs analog signal processing and conversion to 
TTL-compatible output pulses corresponding to Loran-C envelope events. 
For the prototype receiver, a commercially-available microcomputer is 
utilized for both sensor and navigation processing, plus AGC control. 
Computer control and data recording for experimentation are provided 
by a hand-held ASCII terminal and either analog or digital cassette 
units. 


Pilot control of receiver functions is effected through the 
panel-mounted keyboard and video display unit, supported by a video 
processor board with independent memory. The receiver computer is, in 
this prototype implementation, supported by a mathematical function 
processor chip, aiding in the coordinate conversion from Loran-C to 
geodetic coordinates, and for rho-theta conversions. 
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ANTENNA PREAMPLIFIER/ COUPLER 


The Loran-C prototype utilizes the wide-band preamplifier of 
the type reported by Burhans [1, 2, 3] for connection to either a 
1-meter (or larger) whip antenna for laboratory tests or to the 
existing ADF sense antenna aboard the test aircraft. This pream- 
plifier, illustrated In Figure 2, provides -4.9 dB voltage gain at 
100 kHz, matching a High-Z antenna to a 1000 ohm receiver input 
impedance, and a dynamic range from 0.2 to 10,000 microvolts rms at 
the input terminal. Preamplifier 3 db points occur at 10 Hz and 
8 MHz. 


The preamplifier schematic is shown as Figure 3. Input surge 
protection is provided by the NE-2 bulb at the input terminal. The 
low-noise JFET drives an open-collector output stage, the load 
resistor for which is contained in the receiver RF module. Both DC 
power and signal use the signal coaxial cable connecting the preampli- 
fier to the receiver RF section, eliminating power-supply ground-loop 
problems . 

IV. LORAN-C RF PROCESSOR 

The Loran-C receiver RF processor is based on an auto- 
correlation envelope detector. The unit is driven by the low- 
impedance output of the broadband, unity-gain antenna preamplifier. 

See Figures 4, 5 and 6. 

The input circuit is a broad-tuned transformer with a 40kHz 
bandwidth, centered at 100kHz. The output of this transformer is 
placed across a IK ohm potentiometer, the wiper of which controls the 
amount of signal provided to the trap circuitry. Note that this 
voltage-divider circuit has created a passive RF gain control, to 
which the operator has access in the prototype design. 

To improve signal— to— noise ratio (with respect to Interfering 
frequencies) of the Loran-C signal, it is desirable to eliminate 
strong signals close to the 100 kHz region before the RF is actively 
amplified. This is accomplished by passing the RF through a pair of 
notch filters. These narrowband filters are centered at 88 kHz and 
119 kHz, to eliminate the 88 kHz, 116 kHz, and 122 kHz interfering 
frequencies affecting the 99600 U.S. Northeast Chain. These inter- 
fering frequencies are listed in the May 1980 edition of the Loran-C 
User Handbook by the U.S. Coast Guard [4]. 

Current work is underway to investigate distortion products 
from multiple RF signals and harmonic mixing or multiplying to produce 
new frequencies that interfere with Loran-C. Research results and 
documentation of this work will be available from the Avionics 
Engineering Center through the NASA Joint University Program. 

Once the interfering frequencies are minimized, the RF signal 
is amplified by a tuned RF amplifier. This TL072 operational 


amplifier circuit offers +30 dB gain with a 30 kHz bandwidth, centered 
at 100 kHz. The RF signal flow now branches into two separate paths, 
to be manipulated separately to create the RF signals required for the 
LM2 111 FM Detector and Limiter. One path drives a T-notch filter 
tuned to 100 kHz which delays the RF signal and adds the delayed 
reproduction to the actual incoming RF (delay-and-add) . The delayed- 
and-added signal is then amplified by a broadband, +30 dB gain TL072. 
This signal path is terminated at one of the RF inputs to the balanced 
demodulator of the LM2111. The second path drives the RF limiter of 
the LM2111 (the limiter has +60 dB gain), which in turn provides the 
signal for the second RF input of the balanced demodulator and the 
carrier zero-crossing detector. 

The LM2111 IC provides a double-balanced active multiplier used 
as an envelope detector. The demodulated signal is fed to a lowpass 
RC filter (20 kHz bandwidth), creating the Loran— C envelope, from 
which the digital output is derived. An envelope level detector pro- 
duces a digital pulse at the zero crossing of the rising edge of the 
Loran-C envelope. This pulse, correlated with the zero crossings of 
the RF carrier (from the output of the LM2111 limiter) produces the 
digital pulses to be used for signal processing. These pulses are 
stretched to 70 psec. to permit the tracking loop search routines to 
operate efficiently. 

V. TRACKING LOOP HARDWARE 

The Loran-C prototype receiver achieves time-difference 
measurement by use of a software-controlled, shared, digital loop. 

The block diagram for this loop circuit is shown in Figure 7, and 
its schematic in Figure 8. A 1 MHz clock drives the 6-digit BCD 
free-running counter portion of the M0STEK 50395 integrated circuit, 
providing the receiver time base. The 6-digit comparator produces an 
EQ pulse each time the counter and the 6-digit register are identical. 

In operation, the microcomputer loads the register with the 
count corresponding to the desired sample time for loop operation, 
while the counter continues to run. Equality of counter and register 
produces an EQ pulse, which is latched as an interrupt request (IRQ) 
for the microcomputer. At the time the EQ pulse is received, the 
Loran-C digital envelope signal is also latched, and its value made 
available to the microcomputer. The microcomputer may then reload the 
register for the next sample point, a process which has been measured 
to require 450 psec. The loop, therefore, is able to detect each 
envelope pulse, at the 1 msec. Loran— C interval, with sufficient guard 
time to insure correct operation. 

The schematic shows interconnection between the tracking loop 
and the Super-Jolt (TM) microcomputer, which is used in the prototype 
receiver for evaluation. Pinouts are detailed in Figure 9. This 
8-bit microcomputer is based upon the M0S Technology 6502 CPU chip, 
running at a 1 MHz clock rate. In the prototype receiver, the com- 
puter and loop clocks are obtained from the same TCXO for convenience, 
but computer software is entirely asynchronous, interrupt-driven code. 
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The digital loop circuit board is shown in Figure 10. The 
MOS Technology 50395 is visible as the 40-pin chip, supported by a 
7474 latch for IRQ and Loran-C input latching, and one monostable 
multivibrator, 74123, to stretch the Loran-C digital RF processor out- 
puts to 70 nsec. for loop use. The remaining circuitry provides level 
shift services, to interface the MOS loop chip to TTL input/output 
lines for the microcomputer. 

The entire loop circuit is operated by one 6520 PIA interface, 
which is an integral part of the Super-Jolt microcomputer. 

VI. TRACKING LOOP PROGRAMMING 

Initially, tracking loop software lias been developed to 
demonstrate correct hardware operation and to provide a basis from 
which coordinate-conversion and area-navigation routines could be 
developed. This sensor processor software consists of signal acquisi- 
tion, tracking and time-difference generation segments, operating on a 
single Loran-C chain. In the current implementation, three stations 
are tracked, one of which must be the master station. 

Loop routines are initiated at receiver start-up, by user 
choice of Loran-C chain. This selection, made by entering the group 
repetition interval (GRI) , causes the loop routine to add the GRI, in 
microseconds, to the loop register upon receipt of each loop interrupt 
request (IRQ). The result is a series of interrupts, at the GRI rate, 
with samples of Loran-C Input data occurring with each interrupt. 

Loop arithmetic is continuous; that is, the counter is allowed to 
overflow at 10^ counts, with no resulting effect on sample rate. 

Each IRQ causes loop software to read the state of the Loran-C 
latch and to clear the Loran-C and IRQ latches. While in this 
acquisition mode, the goal of the loop routines is to find correlated 
signals at the Loran-C latch, compared with the GRI samples. If no 
occurrence of five contiguous Loran-C ones in 32 GRI frames is found, 
the acquisition segment modifies the sample time by adding 6,500 nsec, 
to the register. This addition effectively delays the sample comb by 
6,500 nsec., and the test for Loran-C signals is repeated. After an 
unsuccessful test at the 6,500 nsec, increment, the test is repeated 
for 36 mini-increments of 33 nsec. each. 

With each Loran— C pulse stretched to 70 nsec., acquisition is 
generally accomplished in less than ten seconds. Once one Loran-C 
pulse is found, acquisition code passes control to station-track code, 
which immediately subtracts 2,000 nsec, from the original register 
contents and repeats the acquisition test. If no pulse is found, the 
register is increased by 1,000 nsec., and the test repeated. If a 
pulse is found, 2,000 nsec, is subtracted again. In this manner, the 
first pulse of each Loran-C station is acquired. 

Fine tracking begins at this point, with the station-track 
routine subtracting one microsecond when the Loran-C data is high, and 



adding one microsecond when the Loran-C data is low. A loop filter is 
inserted at this point, in the form of an up-down software counter, to 
provide for optimum loop characteristics. Filter parameters are taken 
from previous theoretical [5] and experimental [6] work, and were 
chosen as a 5-bit register. When this register overflows, one micro- 
second is added to the loop register. Upon filter register underflow, 
one microsecond is subtracted. The result is a low-pass loop filter, 
with time constant of 4 GRI time periods. Using the first-pulse posi- 
tion in the loop register, determined by the search/track process just 
described, the station-track segment then generates eight sample 
pulses, by 1,000 psec. increments of the loop register, spaced one 
millisecond apart, and permits correction of the loop register value 
at each of the eight pulses. 

After successful lock to the station, the loop register is 
incremented by 9035 iisec., and the master 9th pulse is sought by the 
acquisition test. If found, the current station is labeled as the 
master, for use by background routines in computing time differences. 
The register content, representing the time of occurrence of the 
station first-pulse, is stored for retrieval at the next occurrence of 
the station's signal. 

The loop software then generates a search pulse immediately 
after the last pulse for the current station, and the search/track 
process is repeated. The third station is located in the same manner. 

In addition to the basic pulse-tracking function, the loop 
routines also produce a signal-quality number, useful in generating 
user warnings and assessing receiver operation. A software counter is 
incremented by one, for each station, when the loop routines must 
subtract one microsecond to retain lock. This counter is then cleared 
by background routines after display. In the current implementation, 
the counter is active for ten GRIs, resulting in a total pulse count 
opportunity of 80. In ideal conditions, with no noise present on the 
Loran-C input, the counter should reach one-half the total pulse 
count, since alternate additions and subtractions of one microsecond 
would be required to retain lock on each of the station pulses. 
Significant deviations from this value indicate the potential for loss 
of lock, and may be used to trigger re-acquisition of the station. 

The background routines, which operate with IRQ enabled, exe- 
cute when the signal-processing software just described is completed 
for each Loran-C pulse. A BCD buffer for each Loran-C station con- 
tains the loop register value generated in the previous GRI; these 
values are averaged over ten GRI periods, and the slave values 
subtracted from the master after averaging. The results are displayed 
as the time differences. The signal quality numbers are displayed 
also. Time difference data is placed in buffers for use by the 
latitude-longitude conversion routines. 

The initial implementation of tracking loop programs has met 
the desired objective. The receiver logic has been demonstrated, and 
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loop parameters measured. It has been determined that the single 
microprocessor receiver will be capable of five-station tracking plus 
coordinate conversion, and that ample guard time exists between 
Loran-C envelope pulses for signal processing to take place, with the 
microcomputer operating at 1 MHz. 

A counter/comparator IG offering faster digit-strobe operation 
would be a definite aid, as the scan oscillator on this IC is the 
limiting factor in register loading by the microcomputer. 

The software-controlled tracking loop implementation has 
demonstrated some interesting by-products, in that the loop has appli- 
cations in time synchronization and navigation audio processing not 
contemplated at its inception. Use of the loop in IRIG-B time 
synchronization, for example, has been accomplished with complete 
success. [7] 

VII. VIDEO OUTPUT 

A prototype video interface [8] , designed specifically for 
large-character output for cockpit use, has been used throughout the 
receiver development and evaluation program. Figure 11 shows this 
video circuit board. With the forthcoming addition of rho-theta area- 
navigation software, this video interface will be exchanged for a 
smaller circuit board, able to display both character and graphics 
data. 


Two-page output permits display of Loran-C time-difference 
data, signal quality and housekeeping data on one video page, and the 
latitude-longitude and range-bearing waypoint data on the other. With 
the graphics interface available, CDI information will be displayed on 
both pages, driven by the bearing-distance coordinate conversion 
routines . 


Figure 12 shows the receiver package, with video monitor 
and keyboard installed. The Loran-C RF processor will be enclosed 
under the chassis for isolation from the digital circuitry and the 
video monitor oscillator, and the digital circuit boards will surround 
the monitor chassis. The package is standard general-aviation width, 
for mounting in the vertical stack. The power transformer visible at 
the rear of the unit is installed for bench testing only, and is not 
part of the final prototype, which will operate on 14VDC. 

VIII. LABORATORY AND FLIGHT EVALUATION 

Receiver tests run with a Loran-C simulator have consistently 
provided receiver time differences within +1 psec. of simulator out- 
puts. These tests have also defined the need for receiver AGC applied 
station-by-station, due to TD offsets observed for variations in rela- 
tive signal strength among stations. This AGC circuitry is currently 
under test. 


Flight evaluations, made without AGC circuitry at Ohio 
University Airport, Athens, Ohio on March 9, 1981 are illustrated in 
Figures 13 and 14. These graphs are plots of receiver latitude/ longi- 
tude outputs, converted to range and bearing from a waypoint. 

For Figure 13, the waypoint is the calculated latitude and 
longitude of the runway center point. The flight path begins at point 
(1), on takeoff roll over the waypoint. The path proceeds on climb- 
out, heading 240° and through pattern turns at 800 feet AGL to 
crosswind and downwind legs. Downwind is extended to seven miles 
northeast of the airport, where the aircraft is turned inbound, across 
the UNI non-directional beacon, 5.3 nautical miles from threshold. A 
normal low approach is then executed, heading 243, flying visually 
along the runway centerline. 

Following the low approach, a tight turn is made to a close 
downwind leg, followed by an outbound segment (2) approximately 30° 
south of centerline. A perpendicular cut across the UNI beacon (3) is 
then executed, followed by a turn (4) back to the beacon, and a left 
270° turn to the Inbound approach course. A second low approach along 
runway centerline is then executed, followed by landing and taxi 
operations . 

The data presented is averaged by the receiver over a ten-GRI 
time period; approximately one second on the 99600 chain. Positioning 
repeatability over the UNI NDB and along the approach track are 
observed to be excellent. Long-term variations along the approach 
course tend to be well-correlated, with a peak variation of less than 
600 feet. Note, however, the offset to the north, resulting in a 
track parallel to the desired runway centerline track, of approximate- 
ly one-half mile. This offset has been determined using the labora- 
tory simulator to be due to signal-strength characteristics among 
Loran-C stations in the local area. 

Figure 14 shows Loran-C data, with one-second TD averaging and 
subsequent latitude/longitude conversion, presented as distance and 
bearing from the runway. For this flight, a receiver waypoint was 
determined by placing the aircraft stationary at the runway center 
point, and entering the resulting receiver latitude and longitude 
value as the waypoint. The flight then consists of a takeoff and 
climbout to 1000 feet AGL, on runway heading. Variations noted on 
this plot are a maximum of 0.4 nmi south of course, and 0.3 nmi to the 
north. 


It should be noted that these flight evaluations were local, 
short flights to assess basic receiver operation and raw data stabil- 
ity. Current plans call for documented flights, including ground 
radar tracking for position reference, as soon as AGC circuitry is 
fully tested. 
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IX 


CURRENT WORK IN PROGRESS 


Tracking loop software changes are contemplated for five- 
station tracking with master independence. Decoded warning messages 
will be provided when the master 9th— pulse blink code is available. 

All valid time differences will be made available to the coordinate- 
conversion and area-navigation software, so that composite position 
fixes will be possible. 

Initial flight evaluations have shown encouraging results, 
especially with regard to repeatability. Absolute accuracy in the 
first field tests suffered due to variations in signal strength among 
stations. A computer-controlled AGC, acting on each Loran-C envelope 
pulse and commutated among active stations, has been designed and is 
under test as of this writing. Initial results show marked improve- 
ment, with receiver bias reduced over 75% from operation without AGC, 
on the 9960 chain in southeastern Ohio. 

The latitude/longitude and range/bearing coordinate conversion 
circuitry and software have been used routinely in receiver evalua- 
tion; the documentation for this portion of the receiver will appear 
[9] as a M.S. thesis in the near future. 

Additional flight evaluations are planned, using ground-based 
radar for position reference data. Techniques and procedures for dif- 
ferential Loran-C are also under study, to determine whether receiver 
cost may be minimized by permitting a differential setting periodi- 
cally during a flight, or prior to an approach. 
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Figure 3. Preamplifier schematic. 
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Figure 4. Loran— C autocorrelation envelope detector. 
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Figure 5. Mini-L80 Loran-C detector circuit. 



Figure 6. Prototype Loran-C RF processor. 
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Figure 7. Tracking loop block diagram. 
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Figure 13. Flight test, Loran-C receiver prototype, March 9, 1981. 
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Figure 14. Flight test, Loran-C receiver prototype, March 9, 1981. 

Ohio University Airport, PA-28-180, N56241, 99600 GRI. 
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COMMUTATED AUTOMATIC GAIN CONTROL SYSTEM 


A commutated AGC system for the Ohio University 
prototype Loran-C receiver is described. The 
circuit design, fabrication, and test results 
are presented in this paper. 


Stephen R. Yost 


Avionics Engineering Center 
Department of Electrical Engineering 
Ohio University 
Athens, Ohio 45701 



I. 


INTRODUCTION 


This technical memorandum deals with the commutated automatic gain 
control (AGC) system that has been designed and built for the Ohio 
University prototype Loran-C receiver (refer to Figure I). The current 
version of the prototype receiver, the Mini L-80, was tested initially in 
the summer of 1980. The receiver uses a Microcomputer Associates”' Super 
Jolt microcomputer to control a memory-aided phase-locked loop (MAPLL). 
The microcomputer also controls the input/output, latitude/longitude con- 
version, and the recently added AGC system. (For a more detailed 
description of the receiver operation see ref. 1.) 


The Ohio University receiver uses an envelope generator and zero 
crossing detector to produce a "Loran pulse" which is used by the MAPLL 
to track the Loran station. It was observed that significant errors in 
the time differences occurred when very strong or weak stations were in- 
cluded in the Loran chain being tracked (ref. 2). Experiments with a Loran 
simulator revealed that this error was caused by a phase shift and that 
this phase shift was due solely to the signal-to-signal ratio of the sta- 
tions being tracked. For example, a 10 dB signal-to-signal difference 
produced approximately a 10 microsecond error while a 20 dB difference 

produced approximately a 25 microsecond error. (These results were 

obtained with an Epsco Loran simulator.) 

To reduce or eliminate this error, a commutated AGC was proposed 
at the December 1980 NASA Joint University Program meeting (ref. 3). (See 
Figure 2.) The AGC samples the peak of the envelope for each station 
and stores the resulting voltage on a capacitor. This stored voltage is 
then used to control the gain of the input RF signal. The microcomputer 
switches the AGC sample from station to station. This paper describes 

the circuit designed for the AGC and will also present bench and flight 

test results. The AGC circuit described actually samples starting at a 
point 40 microseconds after a zero crossing determined by the software 
lock pulse ultimately generated by a 30 microsecond delay and add network 
in the receiver front-end envelope detector. Thus this sample point will 
be at about the peak of the ground wave signal and not necessarily at the 
peak of the envelope delayed by strong skywave contamination. Throughout 
this report the reference to "peak of the envelope" has this restricted 
meaning. The whole idea of AGC control is to adjust the level of each 
station signal such that the early portion of each envelope rise is about 
at the same amplitude in the receiver envelope detector. 

The final design is an expansion of the original proposed design, 
(refer to block diagram Figure 3, as well as Figures 4, 5, and 6.) It 
consists of three major parts: A) The sample circuit, B) the DC gain cir- 
cuit and, C) the AGC amplifier circuit. 

II. SAMPLE CIRCUIT 

The sample circuit is a two-stage saraple-and-hold system with three 
separate channels; the switching of these channels is controlled by the 
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microcomputer. All integrated circuits are CMOS, and the supply voltage 
is +12 volts. This supply voltage is necessary because the sample 
voltage levels are in the range of 5 to 7 volts. It is important to note 
that this design has not yet been optimized for minimum chip count. 

The amplitude of a station's envelope varies with that station's 

signal strength (refer to Figures 7 and 8) . The two cascaded 4047 
monostables create a delayed pulse to sample the peak of the envelope. 

The first 4047 triggers on the equals pulse (zero crossing; refer to 
Figure 9), and delays 40 microseconds. The second 4047 triggers on the 
negative edge of the 40 microsecond delay output and produces a 20 
microsecond sample pulse (refer to Figure 10). The delay and duration of 
the sample pulse is optimized to sample the peak of the envelope for the 
weak as well as the strong station. 

The envelope from the RF front-end is passed through a voltage 

follower, 1/2LM353, and into the 4051 analog demultiplexer. The 4051 has 

on-chip address decoding; therefore, the two control lines are decoded to 

0 0 search mode 

0 1 channel one on 

1 0 channel two on 

1 1 channel three on 

Note that the control lines are ANDed with the sample pulse so that each 
channel samples the peak of only one station’s envelope. Each Loran sta- 
tion consists of eight pulses, and, therefore, the station is actually 
sampled eight times (refer to Figure 11) . These eight 20 microsecond 
sample pulses are of sufficient duration to charge the 25 microfarad 
capacitor to the desired final value (the voltage of the peak of the 
envelope). Each station is sampled on a separate AGC channel every group 
repetition interval (GRI). The 4016 analog switch is configured for the 
negative of the logic of the 4051; therefore, the two- to three-line 
decoder is needed. The 4016 controls the "hold" for each channel. When 
a receiver channel is on, the corresponding switch for that channel in 
the 4016 is off. This allows the first stage capacitor to be charged 
while a constant voltage from the previous GRI charge is outputted from 
the second stage capacitor to the AGC amplifier. The .68 microfarad 
capacitor in the second stage of the sample-and-hold is large enough to 
hold a constant voltage for .1 second, the maximum GRI value. The second 
4051 is enabled in the same manner as the first 4051 except that it is 
not pulsed. In addition, a +5 volt signal is applied to the 00 channel 
on the second 4051 to serve as a receiver gain setting used during 
Loran-C station search. The multiplexed output of the AGC sample voltage 
is represented in Figure 12. 

The software changes to implement the procedure outlined above are 
minimal since the receiver operating software tracks each station indi- 
vidually. When the microcomputer starts its search routine, the 00 chan- 
nel is activated to output the constant +5 volt search voltage to the AGC 
amplifier. When all three stations are being tracked, the commutated AGC 



is activated and the envelopes for each station are sampled separately. 

It is not important to Identify which Loran-C station is on a certain AGC 
channel since the microcomputer provides the necessary sychronization as 
a function of basic receiver operation. 

III. DC GAIN CIRCUIT 

The three AGC sample voltages are equal to the peak voltages of the 
envelopes of their respective stations. Therefore, a strong station 
stores a higher sample voltage than a weak station. The voltage- 
controlled amplifier designed for this AGC system requires a higher 
control voltage to amplify the weak stations; therefore, an inverter cir- 
cuit was designed for this purpose. The circuit actually has a two-fold 
purpose, inverting and increasing the gain of the sample voltages. An LM 
353 dual op-amp was chosen for this circuit. Referring to Figure 5, one 
can see that one-half of the chip is used for a voltage follower while 
the other half serves as the gain/inverter. The trimpot on the input to 
the second op-amp controls the gain of the AGC output voltage while the 
other trimpot controls the DC level of the output. The adjustment of 
these trimpots for proper AGC operation will be explained in Section V of 
this report. 

IV. AGC AMPLIFIER 

The AGC amplifier utilizes a CA3028A differential cascade amplifier 
(refer to Figure 6) . The accompanying circuit has been optimized for 
this particular AGC application. (Refer to ref. 4 for a complete circuit 
description as well as test results.) 

V. OPERATION 

As stated previously, the gain and DC level trimpots of the DC gain 
circuit must be adjusted properly for optimum AGC performance (refer to 
Figure 13). As this graph indicates, the best AGC amplifier performance 
lies in approximately the three to eight volt range. The adjustment pro- 
cedure is as follows: 

1. Set the gain to unity (adjust the 10K ohm trimpot to its full 
value) , 

2. Set the DC level to +8 volts and allow the receiver to track 
all three stations and, 

3. Increase the gain from unity until the lowest AGC voltage is 
equal to approximately 4 volts. 

The properly functioning AGC will be similar to Figure 14. This adjust- 
ment procedure outlined above allows for the use of different antenna- 
preamp combinations which may possess different DC components in the 
input RF. It is important to note that once the DC gain circuit is prop- 
erly adjusted to match a certain antenna-preamp combination, it need not 
be adjusted further. 



VI 


TEST RESULTS 


The first test of the AGC system performance was the use of the 
Epsco Loran simulator to provide different signal-to-signal ratios and 
record the results. An example of such an experiment appears in Figures 
15, 16, and 17. A 10-millivolt input signal was used, which is charac- 
teristic of the signal strengths encountered with "live” Loran signals. 
The GRI was set at 99,600 microseconds and the time differences (TD's) 
were set to resemble those received off the air at Clippinger Labs, 
Athens, Ohio (TDY =*42, 594. 3 microseconds, TDZ»56,775.9 microseconds). 
These results are approximated, observed time differences: 


Figure // 

Attenuation 

AGC 

TDY 

TDZ 

15 

none 

on 

42,594.5 

56,776.0 

17 

20dB station Y 

on 

42,605 

56,776 

16 

20dB station Y 

off 

42,622 

56,776 


CONCLUSION: The AGC has little or no effect on "perfect" Loran signals, 

meaning that there is no degradation of performance with the AGC in opera- 
tion. Also, twenty dB of signal-to-signal difference is an extreme case 
which might be encountered only at the limit of a Loran coverage area. 


Following a number of simulator tests, the receiver was tested with 
a live signal. The first step was to obtain an accurate value for the 
correct time differences as recorded at Clippinger Labs. The Loran chain 
used was the U.S. Northeast, GRI=99,600 microseconds. Four receivers were 
tested and the results are as follows: 


Receiver TDY 

Texas Instruments 9900 42,594.4 

Trimble 10A 42,594.3 

TDL 302 42,595 

OU Mini L-80 (no AGC) 42,600 


TDZ 

56,776.0 

56,775.9 

56,776 

56,775 


To obtain enough data points for a good statistical sample, sixty minutes of 
data was collected with the Ohio U. Loran receiver on October 20, 1981 
from 4:00 to 5:00 p.ra. , thirty minutes without AGC, and thirty minutes 
with AGC. Each of these thirty-minute segments was broken into ten- 
minute blocks for a total of six blocks, 550 to 600 data points each. 

The atmospheric conditions were: light cloud cover with moderate 

spherics activity observed on an oscilloscope. A statistical analysis 
package available on Ohio University's IBM 370/158 was used to obtain the 
results shown in Figures 18 and 19. Most of the accuracy displayed is not 
significant but the trends are evident. 


CONCLUSION: The addition of the AGC improved the value of TDY by approxima- 

tely four microseconds. The overall accuracy of the receiver is approaching 
+1 microsecond. One point of special interest is the greater variance of 
TDY with the AGC on. This is due to the occasional sampling of cross— rate 
interference. 
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The most important test of a prototype navigation receiver is a 
flight test. Ohio University's DC-3 flying laboratory made two Loran data 
collection flights on August 29th and 31st, 1981. The flights were center- 
ed around south-central Ohio and in areas of light thunderstorm activity. 
The operational Loran receivers were: TI 9900, Trimble 10A, and 0. U. 

Mini L-80. The results of the August 29th Columbus to Albany via Zanes- 
ville leg are presented graphically in Figure 20. (August 31st data is 
omitted from this report because the results are essentially the same.) 

The 0. U. receiver data was hand-collected and it appears with the flight 
path plotted by the TI 9900. 

CONCLUSION: The plot shows a very close alignment of the two paths. Note 
the slight (less than .5 nautical miles) north bias of the 0. U. receiver 
path. These plots were obtained with latitude/longitude data and not time 
difference data. The latitude/longitude conversion employed in the 0. U. 
receiver does not use any overland propagation delay corrections in the 
calculations. This bias is due mainly to lack of propagation delay 
correction rather than large time difference errors. 

VII. CONCLUSIONS 

The addition of AGC to the 0. U. Loran-C receiver has improved the 
accuracy of the time difference calculations to within approximately + 

1.5 microseconds of the observed time differences for a given position. 

This translates to an improvement of absolute accuracy of approximately 
0.5 nautical mile. Tests of Ohio University's receiver with and without 
the AGC have indicated these results. The majority of error now present 
in the positional data supplied by the Ohio University receiver is due to 
the lack of propagation delay corrections. 

Two additional refinements could improve the performance of the AGC 
system further: 1) a filter to reduce the effect of cross-rate inter- 

ference on the sampling of the envelopes, and 2) an AGC amplifier with 
more dynamic range for an even greater signal-to-signal gain. Complete 
software control would eliminate the adjustments outlined in Section V, 
thus the receiver would require no manual gain adjustments. Other soft- 
ware development could allow for the tracking of all the Loran stations 
in a particular chain. The three most suitable signals would then be 
used to obtain positional data. 

VIII. SUMMARY 

A commutated automatic gain control system has been designed and 
constructed specifically for the Ohio University prototype Loran-C 
receiver. The AGC is designed to improve the signal-to-signal ratio of 
the received Loran signals. The AGC design does not require any analog 
to digital conversion and it utilizes commonly available components. The 
AGC system consists of three major parts: 1) the sample circuit, which 

samples the peak of the envelope of the Loran signal to obtain an AGC 
voltage for each of three Loran stations, 2) a DC gain circuit to control 
the overall gain of the AGC system, and 3) an AGC amplifier to amplify 
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the input RF signal. The performance of the AGC system has been observed 
in bench and flight tests and it has improved the overall accuracy of the 
Ohio University receiver considerably. 
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Figure 1. Ohio University Loran-C receiver. 



LORAN-C RECEIVER COMMUTATED AGC 


• To be Added ha Loran-C Low- Co* t Prototype 

• Commutated, Sampled AGC 

• Avoid* Fronh-End Phase Shift Problem 

• Five Chip* Required for Breodboard 

• Permit* Present AGC for Search Mode 

• Minimum Load on Computer 


Figure 2. Proposed commutated AGC 
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Figure 3. Ohio University Loran-C receiver block diagram with AGC. 



Figure 4. Schematic of the commutated AGC circuit designed for the Ohio 
University Loran-C receiver. All of the IC's are CMOS. This 
circuit is designed to sample the peak of the envelope of a 
Loran pulse to obtain an AGC voltage. 
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Figure 7 . 


Strong station envelope. 


Live signal, M station (Seneca, NY). 



Figure 8. 


Weak station envelope. Live signal, Y station (Carolina Beach). 


149 




Figure 9. 


One ysec equals pulse at zero crossing. 



Figure 10. Twenty-ysec sample pulse at peak of envelope. 
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Figure 11. 


Loran chain as received at Clippinger Labs. 

GRI = 99600 ysec. Note cross-rate interference. 



2 volts/Div. 


GND. REF. 


Figure 12. AGO sample voltages. 
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2 volts/Div. 




Figure 16. Ten-xnV signal, 20 dB attenuation on station Y, AGC off. 
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Figure 17. Ten-mV signal, 20 dB of attenuation on station Y, AGC on. 
Note improvement from Figure 16. 
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Figure 18. Bench test results without AGC. 
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Figure 19. Bench test results with AGC . 


fiPE VGR 

« 



HRSBH CO. JRtKSON C3. 
© 0 


HNN VeR 

< 

Figure 20. DC-3 flight test, 29 August 1981. TI-99C0 flight path is the solid 

OU Loran— C is the line with X*s. Note the slight north bias 
of the OU Loran-C data. 
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A PROTOTYPE INTERFACE UNIT FOR 
MICROPROCESSOR-BASED LORAN-C RECEIVER 


A command entry and display device 
designed to allow convenient operation 
of the Loran-C receiver-processor is 
described. 


Stanley M. Novacki III 


Avionics Engineering Center 
Department of Electrical Engineering 
Ohio University 
Athens, Ohio 45701 




I . INTRODUCTION 

This paper documents an inexpensive data/command entry and 
display system being developed by the Ohio University Tri-University 
group. This system is designed to operate in place of a separate 
ASCII terminal. Also described is the software to interface this unit 
to the 6502-based navigation receiver currently under development at 
Ohio University. 

II. HARDWARE IMPLEMENTATION 

See Figure 1 for an overview of the command entry and display 
logic. In order to retain the use of some of the DEMON(TM) monitor 
facilities provided by the SuperJOLT(TM) microcomputer, an ASCII 
encoded keypad consisting of the decimal digits, decimal point, and 
nine letters has been designed as shown in Figure 2. A printed cir- 
cuit board was prepared to produce the appropriate X-row Y— column code 
appropriate for each character (Figure 3). This X-Y code is input to 
a General Instruments AY-5-2376 keyboard encoder [1]. Figure 4 shows 
the encoder along with an NEC p82I2 octal latch, which stores the 
value of the pressed key until the microprocessor can poll the 
keyboard and read the data. 

The latch holding the 7 bits of character data and a key- 
pressed strobe is made available to a 6530 Versatile Interface Adapter 
(VIA) on the J2 connector of the SuperJOLT microcomputer. Port A of 
the 6530 is configured as an input: lines 0 through 6 carry the ASCII 
value of the key, line 7 the key-pressed strobe which acts as a flag 
during the polling process to indicate that new data is present. Line 
2 of the "B" port is set as an output line; after data has been read 
from the ”A" port, line B2 is toggled to clear the latch so that new 
data can be read. 

A Sony AVF-3250A 4-inch black-and-white monitor designed for 13 
VDC operation was chosen for the display due to the ease with which it 
could be integrated into a standard avionics-size enclosure as shown 
in Figure 5. The monitor accepts standard NTSC composite video 
signals, requires 14 watts at 13 VDC and weighs approximately 1.9kg. 

The VDM-1 Video Display Module [2] allows for the use of alpha- 
numerics and graphic primitives in a 16x16 format along with a 256 x 
256 coarse graphics mode and a 512 x 512 high resolution mode. This 
versatility allows for a variety of alphanumeric, graphic, or 
combined-mode displays. Such capabilities allow for receiver output 
to be displayed in CDI, HSI, or other analog data formats easily 
recognized and interpreted by the pilot. 

III. SOFTWARE INTERFACE 

A monitor routine is being written in 6502 assembly language to 
perform data input and output between the command/display unit and the 
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SuperJOLT microcomputer. The DEMON(TM) monitor routine can provide 
all services necessary for system initialization and data input; 
however, it also requires a full communications terminal and does not 
offer the range of display formats offered by the VDM graphics unit. 
The attendant reduction in size and weight coupled with the ease of 
operation made possible by tailoring the monitor to a specific appli- 
cation make the development of a custom software interface highly 
desirable. 

There are four tasks currently envisioned for the interface 
monitor. They consist of: 

1. Selecting the mode of operation; for example, direct-route 
or multi— waypoint navigation. 

2. Provide user prompts for the data input needed for the 
specified mode of operation. 

3. Provide data conversion from a user-oriented format to a 
microprocessor-oriented one. 

4. Select the display mode of the processor output. 

Tasks 1 and 2 are fairly obvious and no further elaboration will be 
given here. Data format conversion is required because the current 
versions of the time-difference to latitude-longitude and area naviga- 
tion routines require information such as waypoint location to be 
given in a particular 32-bit floating-point format as used by the 
Advanced Micro Devices Am9511A arithmetic processor [3]. This format 
is illustrated in Figure 6. The format conversion typically consists 
of stripping the ASCII zone bits and performing a BCD-to-binary con- 
version and then "floating" the 8-bit integer number into the 32— bit 
format. This data format change is greatly simplified by the arith- 
metic facilities of the Am9511A: software multiply-and-add routines 
are replaced by presenting data to the math chip and giving it the 
appropriate operation codes. This decrease in the size of the inter- 
facing software and attendant improvements in program legibility make 
software maintenance much easier, especially in terras of code optimi- 
zation for faster execution as well as increased memory resources to 
allow for more sophisicated I/O routines. 

Perhaps the most significant difference between the interface 
monitor and the terminal monitor is that the former cannot be 
interrupt driven. The processor's principal function is monitoring 
the LORAN-C pulse train and deriving time-differences from them to 
drive the navigation routines. An interrupt to service something as 
irrelevant to the pulse tracking and TD measurements as a change of 
display formats would force the receiver to retrack the pulse trains, 
reducing the total time the processor is devoting to actual navigation 
duties. As an extension of the data display currently used in the 
prototype receiver, the keyboard will be polled as part of the general 
housekeeping software. The information from the keypad is stored 
until the processor has the opportunity to implement the command. 



IV. 


SUMMARY 


An ASCII keypad with a CRT display capable of alphanumeric and 
graphics-mode operation is being developed to provide specialized data 
entry and display for the Ohio University LORAN-C receiver/processor. 
This unit is being developed to replace conventional communications 
terminals so as to simplify receiver operations to a level typical of 
current avionics systems. 
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Figure 2. Proposed ASCII Keypad Layout. 














Figure 3. Code Assignment Chart - AY -5- 2376 Keyboard Encoder 
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Figure 6. AM951 1 A 32-Bit Floating-Point Format. 


164 



R. F. PROCESSING 


Ralph W. Burhans 

Avionics Engineering Center 
Department of Electrical Engineering 
Ohio University 
Athens, Ohio 45701 



LORAN-C RF PROCESSING 
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AUGMENTED DIFFERENTIATOR FOR PULSE ENVELOPES 


Circuit adapts conventional communications receiver with wideband AM capa- 
bility to detection of pulse signals such as Loran-C. Also enhances reception 
for surveillance and observation of HF over-the-horizon radar signals or others 
where time difference estimates between pulse returns are of Interest. 









SWEEP 

2.0ms/division 

TOP 

Direct Receiver Envelope 

BOTTOM 

Differentiated Envelope 



SWEEP 0.5ms/division 

TOP Direct Receiver Envelope 

BOTTOM Differentiated Envelope 


20:00 EST 11/4/81 SENECA, N.Y., MASTER ON 99600 psec. GRI 


DIFFERENTIAL LORAN-C ENVELOPE SIGNALS 

As observed with YAESU FRG 7700 Communications Receiver, 12 kHz bandwidth, 
slow AGC, total input signal plus noise level = -66 DBM (20 DB over S9). Signal 
lag and subtraction circuit provides differential edge or zero crossings from a 
standard, unmodified communications receiver. Needed for driving external post 
detection microcomputer or digital tracking system for frequency standard or 
simple navigation experiments. 
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10:17AM 26.430MHz 7 minutes 10:19AM 26.320MHz 9 minutes 



10:21AM 26.320MHz 11 minutes 10:24AM 26.320MHz 14 minutes 


H F QTH RADAR RETURNS AT LAUNCH OF COLUMBIA, SHUTTLE II, November 12, 1981 

Top trace is amplitude of received signal. Bottom is differentiated rise time 
from scattered edges of ionosphere returns. Note, slight increase in amplitude 
and intensity of scattered returns progressing to 14 minutes after launch. 
Probable cause is ionized trail from main engine and booster exhaust. The 
pulse GRI was 100 milliseconds (lOms/di vision). Variable 120Hz interference 
present. USSR "Woodpecker" changed channels and direction at about 10:25AM. 
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SOUTH 





22:00 EST llA/Bl 

23.300 MHz GRI 100ms , lOms/div. 

Ums main bang with strong secondary 
scatter at lS-2?ms, 27-3?ms, and ^0-60ms 
after start of main bang. 

TOP amplitude 
BOTTOM DE/DT 



09:37 EST 11/2 Al 

27.900 MHz GRI 100ms , 10ms /div. 

Long PN code sequence, no FM 

TOP amplitude 
BOTTOM DE/DT 



l'>:37 EST llAAl 

1B.S7S MHz GRT looms, hms/div. 

Nominal U to 6ms single pulse with much 
secondarv scatter to 36 ms after main 
bang starts. 

TOP amnlitude 
BOTTOM DE/DT 


USSR HF "WOODPECKER" EXAMPLES WITH FRG 7700 RECEIVER 
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HIGH RESOLUTION VEF SPECTRUM 


Example of reception of weak VLF signals using an HP 35B1A’ Wave Analyzer 
which can detect signals with a very narrow bandwidth of only 3 Hz. The 
graph shows the 50 Hz frequency shift keying format transmitted by VLF 
station "GBR" on 16 kHz from Rugby, England as received in Athens, Ohio. 




PROGRAMMING FOR LORAN-C RECEIVER PILOT DISPLAY 


Fujiko Oguri 

Avionics Engineering Center 
Department of Electrical Engineering 
Ohio University 
Athens, Ohio 45701 



PROGRAMMING FOR LORAN-C RECEIVER PILOT DISPLAY 


PROBLEM 


PROGRESS 


Compute: 

Given: 

Using: 

For: 


Range, Bearing, Ground Speed, & Cross Track Error 
Lat and Long of waypoints and receiver 
6502 CPU and 9511 A math chip 

Final output & pilot display for OU Loran-C Receiver 
Earth Model Chosen 

Kayton & Fried "Avionics Navigation Systems", 

John Wiley & Sons, NY 1969, p. 160. 


Basic equations used check with FORTRAN program for precision. 

Range - aS - £[ < S ln S +sln8 1 > 2 (i=§fi§) + (iiltii^i) 2 (l+eos9K fH-alnO) ] 


(0) Bearing 



cos 0^sin( A-A^) 

cosgsinS^ - singcosB^cosC A~A^) 


) 


where: a =* semimajor axis 

f = flattening of reference ellipsoid 
“ geodetic latitude of the receiver 
(j)^ “ geodetic latitude of the waypoint 
tan0 “ (l-f)tan<f> 

tanBi = (1-fHan^ 

0 ” angle from receiver to waypoint 
A =* geodetic longitude of the receiver 

A^ “ geodetic longitude of the waypoint 


PLANNED WORK 


Convert to assembly language program for 6502 and 9511 chips 


Check precision 

anticipate 0.162° mean bearing error 
0.0106 nm mean range error 


Minimize Time Only 200,000 ysec.(0.2 sec.) available for these computations 
since lat-long and Loran-C receiver TDs take up remainder of 
computer time. 
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After computing range and bearing, other calculations are 
relatively simple. 



PROGRAMMING FOR LORAN-C 


COMPU TE RANGE-BEARING-ETC. 

FROM LORAN-C DATA 

FOR MULTIPLE WAYPOINT POSITIONING 


PROGRESS 

EQUATIONS EVALUATED 

PLANNED EFFORT 

CONVERT TO MACHINE LANGUAGE 
CHECK ACCURACY 


MINIMIZE COMPUTER TIME 




COMPREHENSIVE INVESTIGATION OF LORAN-C 


FOR GA APPLICATION 


Stephen R. Yost 

Avionics Engineering Center 
Department of Electrical Engineering 
Ohio University 
Athens, Ohio 45701 



COMPREHENSIVE INVESTIGATION OF LORAN-C FOR GA APPLICATION 


I. ESTABLISH PERFORMANCE CRITERIA 

A) FAA ENROUTE, TERMINAL, AND APPROACH ACCURACY 
REQUIREMENTS 

B) SUMMARIZE PREVIOUSLY COLLECTED O.U. LORAN-C DATA 

II. DATA ANALYSIS TECHNIQUES 

A) INFORMATION SEARCH TO DETERMINE THE MOST CONCLUSIVE 
METHODS FOR STATISTICAL ANALYSIS 

1) DETERMINE OPTIMUM SIZE OF A STATISTICAL SAMPLE 

B) MULTIPLE, SIMULTANEOUS LORAN-C RECEIVER DATA COLLECTION 
1) UTILIZATION OF CASSETTE TAPE DATA RECORDERS 

C) USE OF OHIO U COMPUTER SYSTEM FOR DATA REDUCTION AND 
ANALYSIS 

III. GROUND MONITORING 

A) LONG-TERM ON~SITE MONITORING OF LORAN"C POSITIONAL 
DATA 

1) EXAMINE DAILY AND POSSIBLY SEASONAL LORAN-C COVERAGE 
GRID WARPAGE 

2) DETERMINE LORAN-C WAYPOINT VALUES FOR RUNWAY 
THRESHOLDS 

B) SELECTION OF TEST SITES 

1) AREAS WHICH ARE COVERED BY MORE THAN ONE LORAN-C 
CHAIN PERMITTING CROSS-CHAIN POSITION COMPARISON 

2) EXAMINE LORAN-C PERFORMANCE IN FRINGE COVERAGE AREAS 

IV. FLIGHT TESTING 

A) USE OF AN ABSOLUTE DATA REFERENCE TO SUBSTANTIATE LORAN-C 
DATA 

B) SIMULATE APPROACHES USING PREVIOUSLY DEFINED LORAN-C 
WAYPOINTS 

0 EXAMINE THE EFFECTS OF NON-GEOGRAPHIC IRREGULARITIES 
ON LORAN-C FLIGHT PATH DATA 

V. CONCLUSIONS 

A) FORMULATE SPECIFIC HARDWARE AND SOFTWARE IMPROVEMENTS TO 
TO THE OHIO U LORAN-C RECEIVER 

B) COLLECTIVELY EXAMINE THE EFFECTS OF OVERLAND PROPAGATION 
DELAY TO LORAN-C SIGNALS 

C) ADDRESS THE POSSIBLE INTEGRATION OF LORAN-C NAVIGATION 
TO CURRENT ATC SYSTEM 



Princeton 

University 




INVESTIGATION OF AIR TRANSPORTATION TECHNOLOGY 


AT PRINCETON UNIVERSITY, 1981 


Professor Robert F. Stengel 
Department of Mechanical and Aerospace Engineering 

Princeton University 
Princeton, New Jersey 08544 


SUMMARY OF RESEARCH 

The Air Transportation Technology Program at Princeton Univer- 
sity, a program emphasizing graduate and undergraduate student 
research, proceeded along seven avenues during the past year: 

• Investigation of Fuel-Use Characteristics of General 
Aviation Aircraft 

• Investigation of a Dead-Reckoning Concept Incorporating a 
Fluidic Rate Sensor 

• Experimentation with an Ultrasonic Altimeter 

• Development of Laser-Based Collision Avoidance Systems 

• Flight Path Reconstruction from Sequential DME Data 

• Application of Fiber Optics in Flight Control Systems 

• Voice Recognition Inputs for Navigation/Communication 
Receiver Tuning 

Principal investigators Robert Stbngel and Larry Sweet have 
worked with Prof. H.C. Curtiss and Richard Miles in advising the 
students participating in this program. 

Navigation research conducted during earlier years provided a 
logical stepping stone to consideration of flight management sys- 
tems for general aviation aircraft, including fuel-minimizing 
guidance between origin and destination points. The effective- 
ness of such guidance logic is, of course, dependent on the accu- 
racy of aircraft fuel-use models, and this, in turn, requires 
adequate mathematical description of engine characteristics. In 
his doctoral research, Richard Parkinson has developed a cruise 
performance model that can be developed from operating handbook 
data, and he has developed a detailed mathematical model of the 
fuel-use characteristics of general aviation reciprocating 
engines. Interim results are documented in a technical paper{l}, 
and his thesis is near completion. Mr. Parkinson began his 
research under the guidance of Prof. Dunstan Graham, who retired 
from the university in 1980; he has continued with the supervi- 
sion of Prof. Curtiss and Sweet. Graduate student Eugene Morelli 
is continuing this research. 
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A simple approach for dead reckoning navigation that is 
enhanced by the use of a fluidic angular rate sensor has been 
studied by M.S.E. candidate Robert Ellis. The fluidic sensor is 
used to correct the northerly turning error of conventional flux 
gate magnetic compasses. (The northerly turning error results 
from the effect of "dip angle" on sensed magnetic heading when 
the aircraft is banked to turn.) Limited flight tests were con- 
ducted during the year, and final results of experimental studies 
will be reported in Mr. Ellis's thesis. Mr. Ellis began his 
research with former principal investigator Dunstan Graham; Prof. 
Sweet is his current advisor. 

Amy Snyder, a junior in mechanical and aerospace engineering, 
completed an investigation of the feasibility of using an ultra- 
sonic transducer as an altimeter for the landing approach through 
touchdown. The transducer, initially developed for use in an 
automatic range-finding camera, demonstrated better than 1% accu- 
racy for ranges of 0.27 to 10.67 m (0.9 to 35 ft) in laboratory 
and wind tunnel tests, with transverse air velocities of up to 
45.7 m/sec (150 ft/sec) and with sensor angles of 0 to 25 deg{2). 
Senior Philip Chu is continuing this research, with the objec- 
tives of doubling the sensor's range and conducting flight tests 
of the system. Prof. Stengel has been Ms. Snyder's and Mr. Chu ' s 
advisor . 

The investigation of laser-based collision avoidance systems 
continued, resulting in several papers, reports, and' theses {3-8}. 
This research is being conducted in parallel with a NASA Ames- 
sponsored study of short-range tracking between two cooperative 
aircraft using low-power lasers. Demonstration equipment has 
been assembled and will be tested in flight within the near 
future. This research has involved Prof. Richard Miles and Larry 
Sweet, graduate students Steven Webb, Edward Wong, and Gregory 
Russell, and undergraduates Leonard Blackburn and Maged Tom eh . 

The Flight Research Laboratory is engaged in a continuing 
study of flying qualities criteria for single-pilot instrument 
flight operations with separate sponsorship from NASA Langley 
Research Center{9}, and elements of the flight path reconstruc- 
tion algorithms, which have common application for the study of 
air transportation technology, have been developed under this 
grant. Techniques for applying extended-optimal smoothing to 
sequential DME data and on-board inertial and air data have been 
developed and applied to flight test data by Ph.D. candidate 
Aharon Bar-Gill. In this approach, a single DME receiver is 
switched between alternate DME stations to update the state esti- 
mate, which is processed after the flight has been completed. 
Because all the measured data are applied to estimate the state 
at each point in the flight path, the result is somewhat more 
accurate than that which would be obtained from the extended-op- 
timal filter. Mr. Bar-Gill is advised by Prof. Stengel. 
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Fiber optic data transmission promises to provide decreased 
weight, improved resistance to electro-magnetic interference, and 
increased reliability in future active flight control systems. 
Graduate student Kristin Farry is nearing completion of an appli- 
cation of fiber-optic technology in a multi-microprocessor sys- 
tem, with the combined sponsorship of this grant and the Office 
of Naval Research. Princeton's Variable-Response Research Air- 
craft (VRA) has dual angle-of-attack and sideslip angle vanes 
mounted at each wingtip, previously connected to the central ana- 
log and digital control logic by standard wiring. Dual micropro- 
cessors have been installed at the wingtips; the sensors are con- 
nected directly to these units, which in turn communicate with 
the central microprocessor via fiber-optic links. The local pro- 
cessors provide scaling and instrument-error correction, elimina- 
tion of roll rate effects on sensed angles, and analog/digital 
conversion before sending the data to the flight control computer 
unit. The fiber-optic multiprocessor system will be tested in 
flight shortly, and the results will be presented in Ms. Farry* s 
M.S.E. thesis. Ms. Farry is working with the guidance of Prof. 
Stengel . 

Voice recognition of pilot inputs could play a major role in 
future air transportation, not only for commercial aircraft but 
for general aviation aircraft as well. This capability has par- 
ticular significance for single-pilot instrument- flight opera- 
tions, where one pilot is required to perform all the tasks nor- 
mally carried out by two or three persons in the larger aircraft. 
The tasks that could be carried out using voice command are var- 
ied; as a generality, these are characterized as the jobs which a 
captain might ask the co-pilot to do, e.g., tuning radios, main- 
taining contact with air traffic control, holding altitude, and 
so on. Frances Koo, a senior in electrical engineering and com- 
puter science, is conducting a project in which voice recognition 
will be used to tune the navigation and communication receivers 
in Princeton’s Avionics Research Aircraft (ARA) . The ARA ' s 
receivers will be tuned by digital signals from a microprocessor, 
which will receive inputs from the pilot via the voice recogni- 
tion board. Ms. Koo is being advised by Prof. Stengel. 

In addition, the research on OMEGA-Dead Reckoning hybrid navi- 
gation completed earlier has been documented in the M.S.E. thesis 
of Ralph Nichols { 10 }, completed this year. Mr. Nichols began his 
research under the guidance of Prof. Graham and completed his 
thesis with Prof. Stengel as his advisor. 

The NASA grant supporting student research in air transporta- 
tion technology has inestimable value in helping educate a new 
generation of engineers for the aerospace industry, and it is 
producing research results that are relevant to the continued 
excellence of aeronautical development in this country. 
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R.C.H. Parkinson, "An Operational Model of Specific Range for 
Microprocessor Applications in Piston-Prop General Aviation 
Airplanes", AIAA Paper No. 81-2330, Nov 1981. 

The paper describes a cruise performance model for conven- 
tional general aviation airplanes of fixed design, powered by 
naturally aspirated spark ignition piston engines and con- 
stant-speed propellers. The model computes the airplane 
operating point which maximizes specific range and is suita- 
ble for airborne microprocessor implementation, with the 
assumptions that the engine ignition timing is optimized for 
maximum engine torque and that fuel-air mass ratio can be 
controlled . It appears that fuel savings of 20% to 26% could 
be achieved through use of this procedure. 

Amy V. Snyder, "Ultrasonic Altimetry: Feasibility and Design 
Principles", Princeton M.A.E. Junior Independent Research 
Project Report, Princeton University, May 1981. 

The report examines the feasibility of using a Polaroid 
Ultrasonic Ranging System as an altimeter for aircraft on 
landing approach. The ranging device generates pulses in the 
50 to 60 KHz band using an electrostatic transducer that 
functions as both an emitter and a receiver of the sound 
energy. Initial wind tunnel tests indicated that the flow of 
air over the unshielded transducer generated sounds in the 
operating band that precluded ranging; however, a megaphone- 
like baffle was found to restore ranging capability with air 
velocities up to 45.7 m/sec (150 ft/sec). This technique 
shows great promise for providing landing approach data at 
low cost, with high reliability, and for aircraft of all 
classes . 

Richard B. Miles, "Laser Beacon System for Aircraft Collision 
Hazard Determination", Applied Optics, Vol . 19, No. 13, July 
1, 1980, pp. 2098-2108. 

A laser beacon collision hazard determination system is 
capable of simultaneously determining range, bearing, and 
heading of threat aircraft. Calculations demonstrate that 
threat aircraft may be observed at >10 km under good visibil- 
ity conditions, the conditions under which the overwhelming 
majority of mid-air collisions take place. A wide variety of 
detection systems may be chosen based on cost, detection 
range, and sophistication. Preliminary tests demonstrate 
that accurate range measurements are possible under daylight 
cond itions . 

Leonard A. Blackburn II, "A Display for a Laser Aircraft 
Position Locator System", Princeton M.A.E. Senior Independent 
Research Project Report, Princeton University, May 1981. 



The purpose of this study was to develop a system that 
displays the relative (three-axis) position between a 
research aircraft and a following helicopter. Signals would 
be displayed in such a way that the helicopter pilot could 
hold position relative to the other aircraft to within one 
foot in distance and one degree in azimuth and elevation 
angles. The study included the selection and construction of 
a data acquisition system, the design of a suitable display, 
and the analysis of position control and pilot-helicopter 
dynamics using root-locus techniques. Helicopter dynamics 
were simulated using an eleventh-order, six-degree-of- freedom 
model; the pilot was modelled by a conventional compensatory 
transfer function. Particular attention was paid to the 
pilot gains and lead-lag compensation required for closed- 
loop stability. The analysis indicated that stabilizing the 
system would require substantial pilot effort at the data 
sampling intervals anticipated . 

5. Maged Tomeh, "Non-Linear State Estimation Algorithms for the 
Development of an Air-Collision Avoidance System Using Laser 
Beacons", Princeton M.A.E. Senior Independent Research Pro- 
ject Report, Princeton University, May 1981. 

An extended Kalman filtering algorithm based on a nonli- 
near dynamic system model was applied to the collision avoid- 
ance problem. The purpose was to assess the level of estima- 
tion accuracy that could be obtained in determining the 
flight path of threat aircraft and to evaluate on-board 
implementation using steady-state filter gains. Over 50 sim- 
ulations of near-collision flight paths were computed . The 
filter was found to possess extreme sensitivity to initial 
conditions and an inherent instability in estimating one of 
the unmeasured state variables. Suggestions for improving 
the filter's characteristics centered on the use of a dis- 
crete linear Kalman filter. 

6. Steven G. Webb, "Aircraft Position Measurement Using Laser 
Beacon Optics", Princeton M.S.E. Thesis, M.A.E. Report No. 
1535-T, Princeton University, Sept. 1981. 

The thesis investigates a system to precisely measure the 
relative position between two aircraft utilizing a laser bea- 
con, an optical detector array, and on-board digital computa- 
tion. The laser beacon consists of two orthogonal fan-shaped 
narrow-width beams with 166-deg coverage, each rotating at 
four revolutions per second. Each of the four detectors in 
the receiving array is composed of two compound parabolic 
concentrators which collect the incoming laser beams over a 
150-deg field of view, collimate the light so that it can 
pass through an interference filter, and concentrate the beam 
on a photodiode. The diode converts the detected beams to 
pulses that are transmitted through a bandpass filter and 
processed for tracking purposes. Initial tests have demon- 
strated the system's potential for use as a position 
measuring system. 
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7. L.M. Sweet, R.B. Miles, S.G. Webb, and E.Y. Wong, "Wide Field 
of View Laser Beacon System for Three-Dimensional Aircraft 
Position Measurement", ASME Paper No. 81-WA/DSC-9, Nov 1981. 

This paper describes a new wide-field-of-view laser beacon 
system for measurement, in three dimensions, of aircraft or 
other remote objects. The system provides aircraft collision 
hazard warning independent of ground-based hardware, and it 
has applications in flight research, helicopter-assisted con- 
struction and rescue, and robotic manipulation. Accurate 
information describing the relative range, elevation, and 
azimuth of the aircraft is generated by the sweep of a low- 
power fan-shaped rotating laser beacon past an array of opti- 
cal detectors. System reliability and performance are 
enhanced through bandpass filtering of the pulse signals, 
digital logic designed to mask spurious signals, and adaptive 
modulation of trigger threshold levels. 

8. L.M. Sweet, R.B. Miles, G.F. Russell, M.G. Tomeh, S.G. Webb, 
and E.Y. Wong, "Digital Detection and Processing of Laser 
Beacon Signals for Aircraft Collision Hazard Warning", AIAA 
Paper No. 81-2328, Nov 1981. 

A low-cost collision hazard warning system suitable for 
implementation in both general and commercial aviation is 
presented. The laser beacon/photodetector system provides 
direct measurement of relative aircraft positions, using 
optimal nonlinear estimation theory. The measurements 
resulting from the current beacon sweep are combined with 
previous data to provide the best estimate of aircraft prox- 
imity, heading, minimum passing distance, and time to closest 
approach . 

9. A. Bar-Gill, W.B. Nixon, and G.E. Miller, "Flying Qualities 
Criteria for Single-Pilot Instrument Flight Operations", 
Princeton M.A.E. Report No. 1528, Dec 1981. 

Preparations for flight test related to the development of 
flying qualities criteria for Single-Pilot Instrument Flight 
Rule (SPIFR) operations are described. The principal objec- 
tive is to examine the effects of aircraft dynamic charac- 
teristics (representative of long period stability and ease 
of longitudinal trim variations) on pilot opinion, pilot 
workload, and flight technical error during simulated IFR 
missions in general aviation aircraft. The Avionics Research 
Aircraft (ARA) has been modified and adapted for SPIFR opera- 
tions. Aircraft configurations to be tested using the ARA’s 
variable-stability fly-by-wire system have been chosen and 
matched in flight. The mission matrix has been designed. 
Microprocessor software for the on-board data acquisition 
system has been flight tested . Flight path reconstruction 
procedures are in a final stage of computer program 
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development. Work has begun on algorithms associated with 
the statistical analysis of flight test results. 

10. Ralph A. Nichols, Jr., "Implementation and Evaluation of an 
OMEGA-Dead Reckoning Hybrid Navigation System" , Princeton 
M.S.E. Thesis, M.A.E. Report No. 1523-T, Princeton University, 
June 1981. 

This thesis reports on the results of flight tests in 
which navigation is accomplished by combining data from an 
OMEGA radio receiver with heading and true air speed measure- 
ments. The implementation and evaluation of the system, 
which was installed in Princeton's LASA 60 utility aircraft, 
are discussed, including fabrication, hardware configuration, 
and software design. Flight test procedures are described, 
the theoretical basis for analysis is developed, and flight 
test results are presented . These flight tests showed 10% to 
40% reductions in rms position errors over that obtained from 
OMEGA data alone. Raw OMEGA position errors typically ranged 
from 0.8 to 1.8 nm (rms). OMEGA signals were found to con- 
tain large amplitude, low frequency noise; flight test 
results were consistent with pre-flight predictions for noise 
of this type. The hybrid system improved navigation continu- 
ity with poor OMEGA signal reception. OMEGA lane count could 
be maintained during as much as an hour of signal loss. It 
was concluded that OMEGA navigation does not provide suffi- 
cient accuracy for general aviation operations over the con- 
tinental United States, although improvements could be made 
in the navigation algorithm. 
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FLYING QUALITIES CRITERIA FOR GA 
SINGLE PILOT IFR OPERATIONS 


Aharon Bar-Gill 
Princeton University 



BACKGROUND 


o Over 200,000 GA aircraft operational in U.S. only 
o Over 40,000,000 flight hours per year 

PROBLEM 


GA accident statistics investigated by NASA (TM 78773} 
and attributed mostly to the SPIFR flight regime - due to 
coupling of the following effects: 

o Low frequency aircraft dynamics deficiencies (no 
design criteria in FAR #23,181, wind shear) 
o High workload environment 
o IFR pilot experience 

o Limited GA instrumentation (compared to airliners, e.g.) 
DIRECTIONS FOR SOLUTION 

o Autopilots (wing levelers, etc., e.g., NASA CP 2170) 
o Training (GA simulators, including IFR conditions, become 
available; e.g.. Aviation Week and Space Tech ., 11/30/81) 
o Instrumentation improvements (advanced displays, DME 
sequencing, Loran C, etc.) 

o Aircraft low frequency dynamic response improvement - 
to date little has been done in this area 

MEANS 

o Theoretical work 
o Ground simulators - done a lot 

o Actual flight testing plus extensive statistical analysis 
much less 
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Modern control theory algorithms are employed to identify 
the most important - in the SPIFR context - aerodynamic con- 
figurations and to implement these configurations on the test 
vehicle . 

The experimental vehicle is the Ryan Navion, which has 
been modified and instrumented to a level of a fly-by-wire 
system, capable of in-flight simulating a wide range of aero- 
dynamic configurations. 

The sensor package includes a strapdown IMU, and airspeed, 
altitude, and aerodynamic angles sensors. The navigational 
package provides VOR/DME measurements. Control manipulation 
time histories are recorded as well. 


SPIFR RESEARCH - FRL, PRINCETON 

o IDENTIFICATION OF KEY AERODYNAMIC CONFIGURATIONS 
o IMPLEMENTATION ON THE ARA IN-FLIGHT SIMULATOR 
o MISSION MATRIX DESIGN 
o EXPERIMENTAL SYSTEM 
o DATA REDUCTION 

o OPTIMAL FLIGHT PATH RECONSTRUCTION 
o REPRESENTATIVE RESULTS 
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LOGICAL FLOW CHART OF THE SPIFR RESEARCH 



QUANTITATIVE AND 
STATISTICALLY SIGNIFICANT 

SPIFR FLYING QUALITIES 
CRITERIA. 
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Formulation of the output command algorithm has been pre- 
sented by Professor R. F. Stengel at the TRI-UNIV Conference 
in Boston, September 1980, and is discussed in the context of 
the SPIFR research in a report entitled Flying Qualities Criteria 
for Single-Pilot Instrument Flight Operations (M.A.E. Report No. 
1528, Princeton University, December 1981). Implementation of 
the implicit model following algorithms may also be found in 
this report. 

Ranges of variation in the aerodynamic parameters are in- 
tended to reflect possible trends in GA aircraft design. 


IDENTIFICATION OF HIGHER-PRIORITY AERODYNAMIC CONFIGURATIONS 


0 OUTPUT COMMAND ALGORITHM 
0 COMPLETE DYNAMIC SIMULATION 

CRITERIA FOR REPRESENTATION OF CHOSEN CONFIGURATIONS 


0 INPUT AND OUTPUT STEADY STATES 
0 TRANSIENT RESPONSE CHARACTERISTICS 

IMPLEMENTATION OF CHOSEN CONFIGURATIONS ON ARA IN-FLIGHT SIMULATOR 

0 IMPLICIT MODEL FOLLOWING ALGORITHM 
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The experimental matrix is designed to provide a 
range of low-frequency aircraft flight characteristics. 
The numbers indicate how many replications will be made 
for each configuration and test subject. The pluses and 
minuses indicate whether the parameter variation will be 
positive or negative. 

EXPERIMENTAL MATRIX DESIGN 

Tradeoff between 

- number of configurations, 

- NUMBER OF REPLICATIONS, 

- NUMBER OF EVALUATION PILOTS 
(UNDER THE CONSTRAINT OF ~ 30 FLIGHT HOURS) 


CONFIGURATION 

PRINCETON 
TEST PILOT 

NASA 

PILOT 

GA 

PILOT 

NO, OF 

MISSION RUNS 

NOM. 

2 

2 

2 

6 


2+ 

2+ 

3++ 

7 

z u 

2+ 

2+ 

i|++ 

8 

Zw 

2+ 

2+ 

3_t 

7 

H u 

2+ 

2+ 


8 

^5T 

3_t 

3_+ 

4++ 

10 

Z 6E 

2+ 

2 ± 

A++ 

8 

z iT 

! z - 

2 i 

1 


8 


L= 62 
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FLIGHT PLAN FOR EACH PILOT IS OF FOLLOWING STRUCTURE: 


TAKEOFF no. 

MISSION no. 

CONFIGURATION NO. 

TRACK NO. 


1 

1 (nom) 

1 

1 

2 

2 

2 


3 

6 

3 


A 

1 

2 

2 

5 

3 

3 


6 

i 

6 (replic.) 

i 

4 


PILOT EXPERIENCE PARAMETER (PE): 


PILOT NO. 

PILOT 

IFR LOGGED 

FLIGHT HOURS/TOTAL TIME 

pc- JUER-TJME 
600 

1 

B. Nixon 

600/5000 

1.00 

2 

D. Peoples 

270/2500 

0.45 

3 

NASA 

i 



o Track and configuration randomization 
o Learning curve effect 
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Navigational stations engaged during this flight test 
series are Robbinsville , Colts Neck, Solberg and McGuire. 


The semi-circles represent covariance contours along 
which the flight technical error in the VOR/VOR or DME/DME 
navigational inodes is constant. GDOP effect (Geometrical 
Dilution of Precision) and the covariance contour analysis 
are discussed in M.A.E. Report No. 1528 (Princeton Univ. , 
Dec. 1981) . 

SPIFR MISSION SIMULATION 
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Each variant has a different vertical profile in addition 
to the horizontal projection differences. The evaluation 
pilot is TAC - vectored on the McGuire AFB tower frequency. 

SPIFR FLIGHT PATH VARIANTS 
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SPIFR ON-BOARD DIGITAL DATA ACQUISITION SYSTEM 
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I 


Subjective performance indicators are the Pilot Opinion 
Ratings (PORs) . CHR is the Cooper-Harper Performance scale 
and the workload scale has been developed by the M.I.T. man- 
machine laboratory . 

We are examining both flight test segments characteristic 
of the SPIFR regime and overall mission performance. This 
provides additional insight and augments the statistical basis. 


EVALUATION SHEET 

MISSION VARIANT # 
CONFIGURATION # 

PILOT 

DATE 

SPEED RETRIMMING 
CHR 

WORKLOAD 

COMMENTS 

HOLDING PATTERN 
CHR 

WORKLOAD 

COMMENTS 

GLIDE SLOPE TRACKING 
CHR 

WORKLOAD 

COMMENTS 

OVERALL MISSION 


CHR 

WORKLOAD 

COMMENTS 
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OBJECTIVE DATA PREPROCESSING 



0 CONVERSION OF 16-BIT BINARY DATA INTO DECIMAL INTEGERS 
0 PHYSICAL TIME VECTOR FORMATTING 
0 CONVERSION INTO VOLTAGE 
0 CONVERSION INTO ENGINEERING UNITS 





As the mathematical model is nonlinear, the extended 
Kalman filter is employed. The stabilized Kalman filter 
formulation is implemented to provide numerical robustness. 

Bias-type states account for directional gyro bias and 
drift, for the DME range error and for the low frequency wind 
gusts. The geoid geometry and Earth rotation effects are 
incorporated in the model . 

The system model has been divided into angular and 
translational submodels, which may be addressed sequentially, 
thus augmenting the numerical robustness. The filter tuning 
methodology takes advantage of the model division. It applies 
theoretical considerations and interacts with the SPIFR generic 
simulation . 


OPTIMAL FLIGHT PATH RECONSTRUCTION 


o Extended, stabilized Kalman filter 

o Rauch-Tung-Striebel optimal smoother 

o Bias-type states and navigational aspects 
(mathematical model extension) 

o Analytical derivation of state and observation Jacobian matrices 
o Division into lower-order submodels 
o Sequential complementary filter tuning 
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The optimal smoother algorithm improves both the state 
estimates and the associated covariances. Bias-type states 
may cause significant estimation errors if their wide-range 
variations are not accounted for. 


DESIGN AND ANALYSIS OF THE OPTIMAL SMOOTHER 

o Filter sensitivity analysis 

o Implementation issues (data compression, nonsimultaneous 

DME lock-on, outliers) 

o Interaction with preflight and flight-testing procedure 
(biases extraction, mission starting point) 
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With regard to the pitch angle and yaw rate time 
histories, "derivative" states are noisier than "integral" 
states, as may have been expected. 

OPTIMAL FLIGHT PATH RECONSTRUCTION - INPUTS 

- Data record length 

- Time vectors of codes identifying the ground 

STATIONS 


station 

CODE 

STATION 

N Wx 

w m° y' 

(FT) 

1 

RBV 

12.1 

29.7 

-250. 

2 

COL 

18.7 

09.6 

-120. 

3 

GXU 

00.6 

35.8 

-120. 

4 

SBJ 

35.0 

44.5 

-190. 

0 

FORR. 

20.8 

36.6 

-109. 


Channel dedicated to station switching timing 


NAVl 

nav2"-\^ 

RBV 

COL 

GXU 

SBJ 

RBV 

EE 

ED 

EB 

E 7 

COL 

DE 

DD 

DB 

d7 

GXU 

BE 

BD 

BB 

b7 

SBJ 

7e 

7d 

7b 

77 


- Approximate initial conditions 




orno -<rn x —4 






With regard to the true airspeed time history, note 

demonstration of characteristic filter lag, which is 

corrected by the smoother. This is a reconstructed measure 

T 

ment from the [u,v,w] states. 


TRUE AIRSPEED RECONSTRUCTION (Sept. 17, 1981) 






This post-flight flight path reconstruction approach may 
be useful in evaluating navigational system performance (on 
board the TCV, e.g.) or to investigate statistically causes 
of flight path deviations if recording altitude, airspeed, and 
DME data in the "black box" . 


SPIFR MULTIPLE REGRESSION ANALYSIS (PHASE I) 


SUBJECTIVE INDICATORS (PORS) 

OBJECTIVE INDICATORS (GS-, H~, V-TRACKI NG DEVIATIONS 

CONTROL MANIPULATIONS 



OVER 

FLIGHT SEGMENT 
OR WHOLE MISSION 
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WORK ACCOMPLISHED TO DATE 


o Integrated flight testing and flight path reconstruction 

METHODOLOGY DEVELOPED 

o High accuracy in trajectory estimation achieved with an 

INEXPENSIVE EXPERIMENTAL SETUP 

o Part of flight test series flown 
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AIR DATA MEASUREMENT USING 
DISTRIBUTED PROCESSING AND 
FIBER OPTICS DATA TRANSMISSION 


Kristin A, Farry 
Princeton University 




Ill I I II II 


A DUAL REDUNDANT AIR DATA MEASUREMENT SYSTEM INCORPORATING DISTRIBUTED 
PROCESSING AND FIBER OPTIC DATA BUSES 


OVERALL RESEARCH GOAL: 

GAIN EXPERIENCE WITH --DISTRIBUTED PROCESSING 

"FIBER OPTIC TECHNOLOGY AND 
"REDUNDANCY MANAGEMENT 
IN THE AIRCRAFT ENVIRONMENT, 


THIS PROJECT: 

DEVELOPMENT OF AN ANGLE OF ATTACK AND SIDESLIP DATA COLLECTION SYSTEM WHICH 
FEATURES: 

—TWO INDEPENDENT MICROPROCESSOR-CONTROLLED DATA 
COLLECTION AND CALIBRATION UNITS 

—TRANSMISSION OF DATA TO THE CONTROL SYSTEM ON A 
FIBER OPTIC DATA BUS 

—SOFTWARE- IMPLEMENTED ERROR DETECTION AND RECOVERY 
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V 



AIR DATA MEASUREMENT SYSTEM 


This figure shows a functional breakdown of the entire 
system. 

Two air data sensor vanes — measuring angle of attack 
and sideslip angle — are mounted at each wingtip of the 
aircraft. Their analog output signal is converted to 
digital format by a multiplexed A/D converter. Next, any 
signal noise and other unwanted contributions such as upwash 
and roll are eliminated from the raw data by software- 
implemented filtering and calibration, so that each wingtip 
sensor unit's data is independent of the other's, thus 
giving the system dual redundancy. Finally, the data are 
transferred via a fiber optic data link to the aircraft’s 
digital flight control system, which uses the system's dual 
redundancy to determine the data's validity before using it 
as inputs to the flight control system. 
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EQUIPMENT INSTALLED IN AIRCRAFT 


This figure shows the location of the equipment in the 
aircraft . 


WINGTIP MICROPROCESSOR PACKAGE: 

Card cage and microprocessor bus — 

Prolog STD (56-pin, tightly structured bus). 

CPU and memory — 

Mostek MDX-CPU2 (4 MHz Z80 CPU with byte-wide 
interchangeable RAM, PROM, and ROM) . 

A/D converters — 

Data Translation DT2742 (12-bit converter with 8 differ- 
ential input channel multiplexing capability) . 

Fiber optic interface -- 

Optelecom 2100 series fiber optic transmitter and 
receiver mounted on a Prolog Utility I/O card. 

Power supply -- 

Power Products DC-DC converters (28 vdc to 5 vdc and ± 

12 vdc) . 

FIBER OPTIC LINKS: 

Cables — 

Siecor 155 (200 micron core all-glass fiber in rugged 
sheathing) . 

Connectors -- 

Epoxy Technology SMA (all-metal epoxy-bonded connectors) . 


MICROPROCESSOR DIGITAL 



WINGTIP MICROPROCESSOR 
PACKAGE WITH 

--CPU AND MEMORY 
— A/D CONVERTERS 
—FIBER OPTIC TRANSMITTER 
AND RECEIVER 
—POWER CONVERTERS 


DATA COMMUNICATIONS STRUCTURE 


This air data measurement system uses two microprocessors 
loosely coupled to a third, the micro-DFCS, which is the 
executive processor of the system. 

The traffic management scheme chosen for data exchanges 
between the microprocessors is a modified polling system. 
The micro-DFCS interrupts the wingtip microprocessors 
whenever angle-of-attack or sideslip data is needed and 
whenever roll rate and velocity data is available. The 
wingtip microprocessors cannnot interrupt the micro-DFCS, 
however, since uncontrolled delays in some flight control 
routines might be undesirable. This bus protocol minimizes 
the data communications overhead for all processors, the 
timing sensitivity, and the propagation of wingtip 
microprocessor problems to the micro-DFCS. 

To further simplify timing, each data point transmitted 
is accompanied by a name. Transmission of a parity bit 
enables the receiving processor to identify most 
transmission errors. 

The code used in transmitting data is Manchester Bi- 
Phase, also known as self-clocking return- to-zero (RZ) code. 
The beginning of each bit interval is marked by a 
transition. If the bit is a one, a second transition will 
occur in the middle of the interval; for a zero bit, the 
signal will remain at the same level throughout the 
interval. This code requires twice the bandwidth of more 
commonly used non-return- to-zero code schemes, but 
eliminates many of their timing problems and lengthens 
transmitter life. 


TRAFFIC MANAGEMENT — POLLING 

MDFCS INTERRUPTS THE WINGTIP MICROPROCESSORS WHENEVER DATA IS NEEDED 
OR AVAILABLE, 


DATA COMMUNICATIONS FORMAT — 



CODE TYPE — 

B I -PHASE SELF-CLOCKING 
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FLIGHT TESTING 


Flight testing of this system will be divided into two 
phases: hardware evaluation and total system evaluation. 

The primary goal of the first phase is verifying that all 
the sytem hardware performs satisfactorily in flight. These 
tests will be done with a minimum of software so that 
hardware faults can be easily separated from software bugs. 
Of course, the communications hardware cannot be adequately 
tested without the data communications software, so the data 
communications techniques just described will be tested 
simultaneously. A secondary task scheduled in this part of 
the testing is collecting some data to check the 
analytically estimated calibration coefficients to be used 
in the system later. 

The second phase of the tests, total system evaluation, 
will focus on verifying the data filtering and calibration 
software as well as the system's fault tolerance. 
Performance criteria for this phase include the accuracy of 
the calibration and the degree of isolation of each sensor 
unit . 


OBJECTIVES 

1) HARDWARE TECHNOLOGY EVALUATION 

—TEST OF HARDWARE AND DATA TRANSMISSION TECHNIQUE 

—COLLECTION OF DATA TO CHECK CALIBRATION COEFFICIENT 
ESTIMATES 

2) TOTAL SYSTEM EVALUATION 

—TEST OF DATA FILTERING AND CALIBRATION SOFTWARE 
—TEST OF SOFTWARE- IMPLEMENTED FAULT TOLERANCE 



PROGRESS SUiiMARY. 


PREVIOUS — 

— FIBER OPTIC LINKS DESIGNED, ASSEMBLED, AND GROUND TESTED 
—MICROPROCESSOR, A/D, AND POWER HARDWARE ACQUIRED 

— SOFTWARE PARTIALLY COMPLETE 

THIS PERIOD — 

— FIBER OPTIC/ MICROPROCESSOR INTERFACES DESIGNED, ASSEMBLED, AND 
GROUND TESTED 

—MICROPROCESSOR, A/D, AND POWER HARDWARE GROUND-TESTED OUTSIDE 
THE AIRCRAFT 

—ALL HARDWARE INSTALLED IN THE AIRCRAFT AND GROUND-TESTED IN POSITION 
—TRANSMISSION SOFTWARE GROUND TESTING IN PROGRESS 
—CALIBRATION AND FILTERING SOFTWARE STILL IN DEVELOPMENT 

FUTURE W0RK-- 

— EVALUATE ALL HARDWARE AND DATA TRANSMISSION IN FLIGHT 
—FINISH DATA FILTERING AND CALIBRATION SOFTWARE 

—GROUND TEST SYSTEM WITH FILTERING AND CALIBRATION SOFTWARE ADDED 

— FLIGHT TEST COMPLETE SYSTEM 
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INPUT/OUTPUT MODELS FOR GENERAL 
AVIATION PISTON-PROP AIRCRAFT FUEL ECONOMY 

L, M. Sweet 
Princeton University 
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A FUEL-EFFICIENT CRUISE PERFORMANCE MODEL 
FOR GENERAL AVIATION PISTON ENGINE AIRPLANES 


SPECIFIC RANGE 


R* 

GROUND MILES/LB FUEL 
INSTANTANEOUS FUEL ECONOMY 


FIXED WING GA AIRPLANES 

• FIXED DESIGN 

• NATURALLY ASPIRATED SI PISTON ENGINES 

• CONSTANT SPEED PROPELLERS 

• MAXIMIZE R* USING AVIONICS AND CONTROLS 



Figure 1 


Block diagram: Cruise performance model 

(standard and nonstandard atmospheres) . 








DEGREES 



1.0 



HELICAL MACH NUMBER M AT X = 0.98 

a ■ ANGLE OF ATTACK OF SECTION ZERO LIFT LINE 
AT X -- 0.7 

X * FRACTIONAL RADIUS ON PROPELLER BLADE 

fcomj* COMPRESSIBILITY CORRECTION TO PROPULSIVE 
v EFFICIENCY 


Figure 4 . 


Propeller 

f 

comp 


performance : 
plot. 


Compressibility correction 


TABLE 1. 

EQUATIONS FOR THE STANDARD ATMOSPHERE 


These equations approxioate the characteristics of 
the U.S. Standard Atnosphere 1962 in the troposphere. 


- H 

P 

■ 1 - 6.87239 x 10 
. 0 5. 25581 


-6 


atnos 


. £_. 4 

p o 6 
- 288.15 6 

« 29.92 6 


( 1 - 1 ) 

( 1 . 2 ) 

(1.3) 

(1-4) 

(1.5) 

( 1 . 6 ) 
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TABLE 2 


EQUATIONS FOR AIRFRAME-PROPELLER-ATMOSPHERE CRUISE PERFORMANCE 


V e - 1.6889 V £ 

C L - 2W/p o V* S 

P AUX " P AUxV^ 

C 

v « v Jm o/n 

w w8 o 

c 

x - v e> /k 7W 

Y « N^WqO/W 

s * . 5252.1 JL |jo 

3 E YZ K 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 


ip 

AB 

B 

A 


AB* 

- 8 



(2.8) 

C cos ip sin 

M'- 

^[(l^Cjdh)**-!] 

(2.9) 

B* ♦ 

AB 



(2.10) 

V 1 

WZ 

♦ tc L 

"3h) 

f co»p 

n> I1V ‘ 

(2.11) 

1 + 

5252.1 

(1 

* * C L dh > f co»pL YZ 

J(2. 12) 



(2.13) 

(2.14) 

(2.15) 


Note: B is required for computing ( COB p- 


NATURALLY ASPIRATED SI ENGINE 
MBT IGNITION TIMING 




CORRECTED AIR MASS FLOW RATE 


Figure 5. Corrected part-throttle engine performance. 
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NATURALLY ASPIRATED SI ENGINE WITH MET IGNITION TIMING 


Figure 6. 


Figure 7 






Wide open throttle performance: Engine shaft speed 

versus engine brake torque Qp 

M ^WOT 


NATURALLY ASPIRATED SI ENGINE 
WITH MBT IGNITION TIMING 



. Wide open throttle performance: Engine shaft speed 

versus inlet manifold absolute pressure. 


Ewqt 


Ne WOT 
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TABLE 3 


EQUATIONS FOR NATURALLY ASPIRATED ENGINE CRUISE PERFORMANCE 


% 

c 


a 

c 




P c * 2TT IQcJ N_/ 33 , 000 

E APA L 


c 


V p 


E 


(3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 


SUMMARY 


1. COMPACT CRUISE PERFORMANCE MODEL 

• CORRECTED QUANTITIES 

• CORRECTED PERFORMANCE PLOTS 

• ALGEBRAIC EQUATIONS 

• MAXIMIZE R* WITH OR WITHOUT CONSTRAINTS 

• APPEARS SUITABLE FOR AIRBORNE MICROPROCESSOR 
IMPLEMENTATION 

2 . HARDWARE REQUIREMENTS 

• IGNITION TIMING REGULATOR 

• FUEL-AIR MASS RATIO CONTROLLER 

• MICROPROCESSOR 

• SENSORS AND DISPLAYS 

3. TYPICAL R* INCREASE = 20% TO 26% ABOVE CURRENTLY ACHIEVED 
R* VALUES. 


J 
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WIDE FIELD OF VIEW LASER BEACON SYSTEM 
FOR THREE DIMENSIONAL AIRCRAFT RANGE MEASUREMENTS 

E. Y. Wong 

Princeton University 


OBJECTIVE 


This is a project for NASA Ames Research Center to develop 
a system that can measure accurately the distance from an air- 
craft to a helicopter for rotor noise flight testing. This 
system will be able to measure the range and angles between 
these two aircraft using laser optics. This system could be 
applied in areas such as collision avoidance, robotics, and 
other measurement-critical tasks. 
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PRINCIPLES OF OPERATION 


A laser beacon is mounted on the front of a helicopter 
which has a pair of sweeping orthogonal fan-shaped beams of 
laser light. This light passes a set of detectors that are 
mounted on the back of another aircraft, which receives a 
set of pulses on each sweep. These pulses determine the vector 
orientation of each rotating beam; using this pair of vectors 
the vector pointing from the array to the beacon can be found. 
True range measurements can be calculated once the sweep rate 
and array geometry are known. This position data is then 
displayed to the test engineer in the aircraft, and error signals 
are sent to the helicopter pilot. 
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LASER BEACON SUBSYSTEM 


The laser beacon consists of two 5 mW He-Ne lasers, 
optics, and motor drive. The two lasers are placed under- 
neath the beacon and brought to the top by means of prisms. 
Beam expansion optics are used to expand the beams for the 
cylindrical mirrors. These mirrors are used to reflect the 
expanded beams into the fan-shaped light beacon. The motor 
drive unit powers two orthogonal shafts in which these mirrors 
are mounted. The drive will rotate the mirrors giving sweep 
rate from two to four revolutions per second. A belt drive 
and flywheel are added for increased noise reduction from the 
motor . 
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ROTATING 

MIRROR 
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DETECTOR AND ARRAY SUBSYSTEM 


The detectors are a pair of compound parabolic 
concentrators mounted end to end. This configuration results 
in a non- imaging wide-angle lens . A narrow-band interference 
filter is placed between the concentrators to reduce solar 
interference. A photodetector is placed at one end to 
collect the resulting light from a beam sweep. Further 
electronic filtering is done and is made compatible with our 
timing circuitry. The detector array that we are testing is 
a right tetrahedron of four detectors. This gives a simpler 
mathematical relation between the rotation axes of the beacon 
and the detector array. 
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MICROPROCESSOR SUBSYSTEMS 


A microprocessor is used to implement the vector 
relations and to do further signal conditioning on board 
the detector-carrying aircraft. This computer is also 
transmitting the error signals to the helicopter, driving 
the I/O interface with the test engineer, and outputting 
the measurements to a data recording device. This system 
consists of a 4 MHz Z-80 microprocessor, a 9511 math chip 
for high speed floating point arithmetic, an 8253 counter/ 
timer board for 0.5 microsecond timing resolution, and 
other peripheral drivers. A radio link is put into the 
system to transmit serial data from the aircraft to the 
helicopter. Another microprocessor is used in the helicopter 
for driving the displays. A glide slope indicator displays 
the Y-Z error measurement and a LED panel meter displays the 
X error measurement. This microprocessor will also have 
other warning lights and error checking functions. 
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TESTING 


Static tests were run to examine system performance 
without aircraft dynamics. These tests were used to determine 
the current noise levels in the complete system. The follow- 
ing test results have been obtained: 

1) A 1-foot standard deviation with approximately 
100 feet beacon to detector array distance 

2) A 1/2 degree standard deviation in angular 
measurements . 

In these tests there was no signal conditioning done by the 
microprocessor . 



FUTURE PLANS 


With the conclusion of the static tests, a digital 
filter can be introduced into the software to improve system 
performance. Flight tests then will be run at NASA Ames 
Research Center for actual measurement capabilities. Further 
modifications then can be made to the hardware and software to 
meet or exceed the design specifications. 
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